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DESIGN ANALYSIS OP A SUBSONIC NUCLEAR -POWERED LOGISTIC 

AIRPLANE WITH HELIUM-COOLED REACTOR* 

By R. H. Cavicchi, H. H. Ellerbrock, E. W. Hall, H. J. Heppler, 
J. N. B. Livingood, and F. C. Schwenk 




SUMMARY 


A detailed design analysis of a nuclear-powered turboprop logistic 
airplane with a helium-cooled reactor was made to determine how the 
nuclear and engine components might look, whether any unforeseen dif- 
ficulties might arise in the design, and what research is most needed 
to support such a system. A 400, 000-pound- gross -weight airplane designed 
for 0.72 flight Mach number at 30,000 feet was selected as the reference 
airplane. The resulting wingspread and fuselage length are 250 and 175 
feet, respectively; and reactor power of 98.5 megawatts is required. 
Chemical fuel is provided for thrust augmentation for takeoff, landing, 
and emergency use. For emergency, part of the neutron shield is assumed 
to be chemical fuel. Each of eight 4600-horsepower turboprop engines 
required is supplied with air at 1800° R turbine-inlet temperature from 
its separate heat exchanger. Mounted in the wings, the heat exchangers 
receive helium at 2250° R from the reactor. A helium pump, absorbing 
about 24 percent of the engine power, is required with each heat- 
exchanger unit. 

The airplane can be either logistic (carrying passengers and/or 
freight) or patrol type. A unit shield designed to limit the dose rate 
in the crew compartment to 0-025 rem per hour is provided. This limit 
allows approximately 40 hours of flight every 2 weeks. On the basis of 
the assumptions of the analysis, the following weight breakdown seems 
reasonable: Aircraft structure, 108,000; equipment, 22,000; powerplants, 

83,000; shield and reactor, 120,000; payload, 63,000; and chemical fuel 
for takeoff and landing, 4000 pounds. The heat-exchanger temperature 
limits the performance. Research on the reactor (especially the fuel 
elements), the heat exchangers, the pumps, and the helium leakage prob- 
lem is required. 


*Title, Unclassified. 
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INTRODUCTION 

The direct-air, liquid-metal, and circulating-fuel reactor cooling * 

systems for aircraft nuclear powerplants have been under investigation 
for several years. However, studies of the gas-cooled reactor in a 
closed cycle have, mainly been limited to cycle analyses. As a gaseous • 

coolant, helium is very attractive because it is inert and has excellent 
properties for heat transfer. Analysis of the performance of direct-air, 
liquid-metal, and helium-cooled systems for supersonic, turbo jet -powered, 
logistic aircraft indicates that helium compares very favorably with the 
other systems (ref. l) . However, no known detailed study was available 
on which to base many of the assumptions that were required in this anal- 
ysis. Therefore, it was decided to make a more detailed design study of 
a helium-cooled aircraft nuclear powerplant to determine how the compo- 
nents might look, whether any unforeseen difficulties might arise in the 
design, and what research is needed most to support a system of this 
kind. 

The choice of airplane type lay between a supersonic airplane and 
a subsonic airplane that may have less value but would be considerably 
easier to design and operate. The power requirements for a supersonic 
airplane are several times those of a subsonic airplane so that the 
weights and performance of the power system are much more critical in a 
supersonic design. The first nuclear-powered airplane, regardless of the 
type, will have many problems that cannot now be anticipated in the design. 
Accordingly, a conservative approach appeared to be justified, and a 
subsonic airplane was selected for the design analysis. 

For subsonic speeds, unless speeds near Mach 1 are required, the 
turboprop engine is much more efficient than the turbojet engine (see 
ref. 2) . Furthermore, the turboprop engine and propeller characteris- 
tics provide increasing thrust with decreasing speed, so that, if the 
engine lacks sufficient power to fly at the design speed, flight may 
still be possible at a lower speed. The characteristics of the turbo- 
jet engine, however, are such that it would be in more serious trouble 
under similar conditions. In addition, the turboprop has better take- 
off characteristics. Therefore, the turboprop engine was selected for 
the nuclear-powered airplane of this study. 

The present analysis is a design study of the airplane and the 
various components, rather than a cycle study of the effects of changing 
various parameters. In most instances the parameters are arbitrarily 
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chosen and in only a few cases are the effects of changing these param- 
eters investigated. The report is divided into discussions of the 
major airplane components, and the effects of changes in design or oper- 
ating conditions investigated are discussed in connection with each 
component . 

All methods used to calculate performance are reported in appen- 
dixes. Whenever existing methods are available, they are referenced. 

In many cases these methods are extended or revised to fit the condi- 
tions of this analysis. In some cases they are generalized to cover a 
wider range of conditions; these general analyses are included in the 
appendixes. The reactor components and assembly procedure are described 
in appendix I, by D. W. Drier. 

SYMBOLS 

The following symbols are used in the body of this report : 


Arp 

turbine frontal area, sq. ft 

AFD1,2,3,4 

annular flow dividers 

AR1,2,3,4 

annular reflectors 

AS 

annular support 


AS1,2,3,4,5,6 annular shields 


D 

airplane drag, lb 

d 

reactor core diameter, in. 

E 

energy, ev 

P 

net jet thrust, lb 


FS1,2,3,4,5 front support shields 

HP shaft horsepower 

ISP1,2,3,4,5 inner shield plates 


k eff 

effective multiplication factor 

LCFPP 

low cross-section fission products poison 

l 

reactor core length, in. 


OTPrmTPrn 
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Mq 

OSH., 2 
P 
Po 
p(t) 

Q 

RAS1,2,3,4,5 
RSI, 2, 3, 4 


T 

e 


T, 


He 


m 


W 


w_ 


w 


an 


v He 


flight Mach number 
outer shield plates 
helium pressure, lb/sq ft 
equilibrium poison (during operation) 
poisoning at time t (during operation) 
volumetric heat source, Btu/(sec)(cu in.) 
rear annular shields 
rear shields 

temperature of inner cladding surface of circular hole 
through fuel element (see fig. 25), °R 

temperature of cladding surface at inner radius of 
annular passage through fuel element (see fig. 25), 
°R 

temperature of moderator block surface at outer radius 
of annular passage through fuel element (see fig. 
25), °R 

temperature at outer boundary of moderator block (see 
fig. 25), °R 

temperature of helium, °R 

maximum temperature of fuel-element meat , °R 
time, sec 

uranium volume content per unit volume of compound 
weight , lb 

rate of airflow through engine, lb/sec 

rate of helium flow in annular passage through fuel 
element, lb/sec 

rate of helium flow through reactor, lb/ sec 
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w h 

* 

rate of helium flow in circular passages through fuel 
element, lb/sec 

Xe o 

xenon concentration at equilibrium, nuclei of 
xenon/ (sec ) (cu cm) 

r 

fission yield 

CO Sk 

ao 

CVl 

* 

reactivity 

decay constant, sec - '*' 

2 a 

macroscopic absorption cross section, cm - '*' 

2 f 

macroscopic fission cross section, cm - '*' 

a a 

microscopic absorption cross section, barns 

Of 

microscopic fission cross section, barns 

o 

9- 

equilibrium neutron flux, neutrons/ (sec) (cm ) 

9 r 

r 

relative neutron flux, neutrons/ (sec) (cm^) 

^1,2, 3, 4, 5, 6 

neutron flux for group 1, 2, 3, 4, 5, and 6, 
neutrons / ( s e c ) ( cm^ ) 

♦ 

turbine blade metal taper factor 


Subscripts : 


D 

decal in 

eq 

equipment 

G 

gross 

I 

iodine 135 

n 

neutron shield 

P 

payload 

Pm 

promethium 149 

PP 

powerplant 

R 

reactor 
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S 

st 

t 

th 

U 

Xe 

r, i 
r, o 

Superscript : 
FPP 




shield 
structure 
total 
thermal 
uranium 235 
xenon 135 

inner gamma shield 
outer gamma shield 

fission products poison 


GENERAL REQUIREMENTS 

The mission for which this nuclear-powered airplane is considered 
is logistic, for carrying heavy loads of either equipment or personnel 
over long distances. The increase in thrust with decrease in flight 
speed of the turboprop engine chosen for the airplane provides a per- 
formance margin that would be available if required. Thus, if the en- 
gine is designed for a fairly high subsonic flight Mach number, flight 
might still be possible at a reduced speed if less than design power is 
available from the engine or more than design power is required for the 
airplane. The design flight Mach number for this mission is therefore 
0.70 or higher. Prom considerations of desired Mach number, wing load- 
ing, and ambient air temperature, the design altitude was arbitrarily 
selected early in the investigation as 30,000 feet. Later, the effects 
of changing design altitude were investigated. The results are included 
in a later discussion. 

The characteristics of the nuclear-powered airplane and shield are 
such that the payload fraction of the gross weight increases rapidly 
with increasing gross weight for airplanes near the size and weight of 
existing ones. Nevertheless, in spite of the increased carrying effi- 
ciency with increased size, it is desirable to limit the size and weight 
to even less than those of some existing airplanes, provided that at 
least a reasonable payload can be carried. Larger gross weights would 
seriously l imi t the number of fields from which the airplane can operate; 
and, with the usual growth factors that accompany the design and use of 
an airplane, the gross weight would become prohibitive. Furthermore, in 
keeping with the conservative approach, more knowledge is available for 
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the design and construction of airplanes within the existing size class. 
For these reasons the design gross weight is to be kept under 400,000 
pounds if a payload of the order of 50,000 to 100,000 pounds can be 
carried. 

For similar reasons, the engine size is kept within the size of 
existing engines. Although this probably results in more engines on the 
airplane than is generally desirable, the reliability obtained from en- 
gines within the existing size class probably more than offsets the dis- 
advantages of multiengine installations. Engine size is based on the 
sea- level airflow. A value of 140 pounds per second at standard sea- 
level conditions was chosen as the max imum allowable airflow for an en- 
gine. For the assumed turbine-inlet temperature (1800° R), this airflow 
limit ensures that shaft power will not exceed 5900 horsepower. For a 
400,000-pound airplane with the speed and altitude considered, this lim- 
itation in size results in a requirement of eight engines. 

Except for the reactor, the heat exchangers are probably the most 
critical items in the airplane. The maximum turbine-inlet temperature 
of the engines is limited by the heat exchangers to considerably less 
than that employed in some existing engines. The engines then are con- 
servatively designed with respect to turbine-inlet temperature, which 
should contribute to their reliability. 

The payload-carrying ability of the airplane might also be increased 
by using a lighter-weight divided shield rather than a unit shield. How- 
ever, because of the increased radiation in the cargo compartment, its 
use would seriously restrict the usefulness of the airplane by limiting 
the payload to materiel or equipment and by complicating the transfer or 
loading operations and the airplane and engine maintenance and overhaul 
procedures. Since not enough is known about handling procedures, it was 
decided to concentrate on a unit shield. 

Some proposals for nuclear-powered airplanes are based on the use 
of two reactors to improve the safety. Providing the airplane with twice 
the power required to fly it imposes a serious penalty, since the two 
reactors must be independent; that is, each must be capable of critical 
operation without .the other, and this results in increased shield weight. 
When power in only one of the reactors is inadequate to fly the airplane, 
sustained flight would be possible only by using chemical fuel to provide 
the additional heat to the air. 

A more reliable system for the same weight, but which possibly pro- 
vides less emergency range, is obtained by eliminating the second reactor 
and shield and replacing them with an equivalent weight of chemical fuel. 
There are probably many instances, and perhaps even whole flights, where 
flight on chemical fuel alone is desirable (e.g., takeoff or landing, 
training, ferrying, etc.). Furthermore, this chemical fuel could replace 
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part of . the reactor shield, because, if the reactor were shut down, full 
reactor shielding would no longer be necessary. Because of the added 
complication of two-reactor design and operation, a single reactor is 
used in this airplane, and chemical fuel is used for emergency flight. 

The amount of emergency fuel is limited to that required in the reactor 
shield. Some additional chemical fuel is provided for takeoff and 
landing. 

It would, of course, be desirable to have to refuel the reactor 
as infrequently as possible. Refueling will Involve dismantling the 
reactor and shield and removing many radioactive components from the 
reactor. This operation should be avoided as much as possible. On the 
other hand, designing for a large fuel burnup not only makes the design 
more difficult, but also leads to more uneven power distributions and 
therefore more severe operating conditions. Since the reactor, like 
the other airplane components, should be inspected within reasonable 
intervals to forestall any failure, an operating time of about 1000 
hours at full power was selected as the burnup time for the fuel ele- 
ments. On a structural basis the reactor lifetime was taken as a min- 
imum of 10,000 hours. Although the reactor would probably not be oper- 
ated for this length of time (10 refuelings), this tends to provide some 
margin of safety from a structural standpoint. 

The missions to which this airplane Is best suited are exceptionally 
long range, perhaps beyond those normally possible with chemically fueled 
airplanes. If a total range of 16,000 nautical miles is assumed, the 
trip would require about 40 hours. One large unknown is the waiting and 
turnaround time; but, if it is further assumed that one such trip could 
be made every 9 days to 2 weeks, then about 30 to 50 weeks would elapse 
between refuelings of the reactor. Therefore, it is estimated that the 
reactor should be refueled once each year, after it has operated about 
1000 hours and has flown about 400,000 nautical miles. 

Although the calculation of crew dose rates is an uncertain proce- 
dure at the present time, an attempt was made to calculate the shield 
requirements, and a choice of allowable dose rates was therefore neces- 
sary. For an assumed dose rate of 600 millirems per week, which is 
double the allowable laboratory dose rate of 300 millirems per week, the 
previously assumed 40 hours of flight every 9 days to 2 weeks results in 
an allowable dose rate during flight of between 20 and 30 millirems per 
hour. For the present case a dose rate of 0.025 rem per hour was assumed. 


GENERAL ASSUMPTIONS AND TYPICAL AIRPLANE 

Throughout this study many arbitrary assumptions had to be made be- 
cause there was insufficient time to investigate all parameters and com- 
binations of parameters. Thus, even though the effects of varying many 
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conditions were studied, it was necessary to set standard or reference 
conditions fairly early in the investigation in order to establish a 
basis for comparison and to study in greater detail the design of the 
reactor and the heat exchanger. Considerably more effort was therefore 
spent on the design of this reference airplane than on any other; how- 
ever, this should not be taken to mean that the design and the operating 
conditions for this particular airplane are considered superior to 
others that may or may not have been investigated. In fact, later anal- 
yses showed improved performance for other values of several parameters. 

Because the design is considered to be fairly representative of the 
type of airplanes and conditions investigated, this reference airplane 
will hereafter be referred to as the typical airplane. The principal 
assumptions, specifications, and operating conditions for this typical 
airplane are given in table I. The weight breakdown of the various com- 
ponents is given in table II. An artist's sketch of the airplane is 
presented in figure 1, and a picture of a model showing the location and 
arrangement of equipment is presented in figure 2. The general over-all 
picture of the design and performance of the typical airplane is pre- 
sented here in order to provide orientation for later discussions of the 
various components. 

The airplane is a conventional subsonic type with control surfaces 
on the tail in the rear. The reactor and shield are housed in the rear- 
ward part of the fuselage, with helium lines running through the wings 
to eight heat exchangers located in the wings, each directly over an 
engine. Eight turboprop engines are located on the underside of the 
wings. The engine and heat-exchanger layout is shown in figure 3. Air 
is collected behind the engine compressor, ducted upward and through 
the heat exchanger in the wing, then through an afterburner- type com- 
bustor, and thence down and Is delivered to the turbine. No fuel is 
burned in the combustor during normal operation, but the combustor is 
provided to allow for emergency or off-design operation. An afterburner- 
type combustor is used so that very small pressure losses will be in- 
curred. Therefore, air can be passed through the combustor at all times 
and no bypass ducting or valving need be provided. 

By locating the heat exchangers in the wings away from the engines, 
the engines can be dismounted or changed without disrupting the helium 
system by disconnecting the air ducts and helium-pump drive shafts. The 
helium pump is located beside the heat exchanger and driven from the en- 
gine shaft. The helium lines between the heat exchanger and reactor are 
concentric, with the hot line coming from the reactor located inside the 
annular passage carrying the cooled helium back to the reactor. Heat 
transfer between the passages is reduced by using a multiple-walled duct 
with a corrugated inner wall between two smooth outer walls. A sketch 
of the inner pipe showing the corrugations between the two walls is 
shown in figure 4. 
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Each engine has a sea-level airflow of about 110 pounds per second 
(which is below the limit of 140 lb/ sec established previously) and a 
sea-level shaft power (including helium pump power) of 4600 horsepower 
at the rated turbine-inlet temperature of 1800° R. Six engines could be 
used if engines were available having a sea-level airflow of about 150 
pounds per second and a shaft power of over 6000 horsepower. Temper- 
atures higher than 1800° R could be used in the engine except for the 
added difficulty in the design of the heat exchanger. The engine is a 
one-spool type with the propeller; helium pump, and compressor all driven 
by one turbine with three stages. The compressor has 11 stages and a 
design-point pressure ratio of 11, which was selected on the basis that 
this is the highest value that can be obtained practically in a one- 
spool engine. 

The reactor is of two-pass construction, having a flow divider sur- 
rounding the core. An outline of the reactor is shown in figures 5 and 
6 (plan views) and 7 (which is section C-C of fig. 5). Helium enters 
the reactor through the shield at the rear face, passes through the re- 
flectors and thermal shields for cooling, then reverses and flows through 
the core containing the fuel elements and moderator, and thence back out 
the rear face of the reactor in the center of two concentric helium 
ducts. Figures 8 and 9 show closeup pictures of the reactor pressure 
shell in the airplane model with the external shield indicated. 

The reactor was originally designed on the basis of 100 megawatts 
of heat given up to the coolant. Later it developed that about 98.5 
megawatts were required at the design point. Helium enters the reactor 
at 1250° R (pump outlet temperature) , is heated to about 1300° R in 
cooling the reflector and shields, and is further heated to the final 
temperature of 2250° R in the core. The helium pressure entering the 
reactor is assumed to be 180,000 pounds per square foot. A total- 
pressure loss of 4 percent was assumed for the helium in passing through 
the reactor; of this loss, 2.6 percent occurs in the core. 

The fuel elements consist of tubes containing a total of 153.2 
pounds of uranium carbide UC (129.7 lb q n a molybdenum matrix with 

a molybdenum cladding both inside and outside. Helium flows both inside 
the tube and outside through the annulus between the tube and the moder- 
ator. The moderator consists of hexagonal blocks of beryllium oxide. 
Molybdenum is used as the supporting structure for the core. The side- 
reflector material is assumed to be beryllium metal, and the end reflec- 
tor beryllium oxide. The thermal shields and pressure shell are made of 
mild steel and type 347 stainless steel, respectively. The outside diam- 
eter of the pressure shell is about 64 inches. Further details of the 
reactor are given in a later section. 

There are 13 control rods of the following types located as shown 
in figure 7: two regulating rods, six shim rods, three shim-safety rods, 

and two scram rods. The six shim rods are placed in three pairs, each 
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pair operated so that the two rods move in opposite directions with a 
variable loading in each rod, thus tending to maintain an undistorted 
axial flux distribution. The six shim and two regulating rods are 
actuated by electric motors . The three shim-safety rods are hydrauli- 
cally operated, and the two scram rods are pneumatically operated. 

The biological shield (see figs. 8 and 9) surrounding the pressure 
shell consists of an inner gamma shield of iron or borated steel and an 
outer neutron shield consisting of lithium hydride and decalin (used 
for emergency chemical fuel) located as shown in the figures. Tanks of 
decalin are also used to reduce duct streaming. Chemical fuel for use 
in takeoff and landing is stored in tanks in the wings; the shielding 
decalin is used for emergency only. 

The payload is considered to consist of two parts: one part having 

a high density, such as materiel or equipment, and the other part hav- 
ing a low density, such as personnel or very bulky but lightweight ob- 
jects. The high-density part is located near the center of gravity of 
the airplane just ahead of the reactor and shield. Ahead of this is 
located the compartment that carries the low-density payload. Figure 
1 shows the airplane equipped with two decks of seats in the low-density 
compartment arranged for carrying personnel. 

For the assumed 400,000-pound airplane, the aircraft structure 
weighs about 108,000 pounds or 27 percent of the gross weight. In addi- 
tion, there are 22,000 pounds or about 5 percent for equipment weight, 
which includes the crew, instruments, and so forth. The poweiplant 
system, including the engines, heat exchangers, lines, and so forth, 
weighs about 83,000 pounds or 21 percent of the gross weight. Of this 
weight about one-third is attributed to the heat exchangers and lines. 
Four thousand pounds of chemical fuel for takeoff and landing are con- 
tained in tanks in the wings. The reactor and shield weigh 120,000 
pounds or 30 percent. The remaining weight is available for a payload 
of 63,000 pounds. 


AIRFLAME 

Only a cursory analysis was made before the design conditions for 
the typical airplane were selected. During the time that the typical 
airplane was being studied in detail, analyses were made to determine 
to what extent the typical airplane payload could be increased through 
the choice of different design conditions for the airplane and its com- 
ponents. In the case of certain of the design conditions, it was pos- 
sible to incorporate better choices into the typical airplane design 
without disturbing previously completed work; however, most of the de- 
sign conditions remained fixed at the initially assumed values. This 
section presents the results of the design- point studies that are related 
to the performance of the airplane. 
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Analysis 

The lift-drag ratios were computed by conventional methods . These 
are described in appendix A along with the methods for determining the 
structural weights. As mentioned in appendix A, structural weights de- 
pend on the distribution of the loads between the wing and fuselage. 
Consequently, the variation of structural weight with the load distribu- 
tion was calculated and plotted for each airplane geometry. These curves 
were used to determine the structural weight after calculations of the 
weights of the components attached to the wings (engines, nacelles, heat 
exchangers, lines, and wing-mounted landing gear) were completed. Meth- 
ods for calculating component weights are given in the appendixes. Pay- 
load is defined as the portion of the gross weight Wq. remaining after 

all component and structure weights are accounted for. 

The results of the analyses are presented as the variation of com- 
ponent weights and payload with one particular design variable, as in 
figure 10. All other conditions remain fixed at the values selected for 
the typical airplane . The airplane weights are divided into four main 
categories : (1) structure and equipment weight ( , including crew, 

instruments, and controls; (2) powerplant weight Wpp, including engines, 
nacelles, propellers, heat exchangers, lines, gearboxes, and helium 
pumps; (3) reactor and shield weight W R+ g; and (4) payload weight Wp. 

Part of the shield consists of a chemical fuel that would be avail- 
able for use in the combustors for emergency operation (reactor shut- 
down) . The weight of the fuel portion of the shield Wp is also indi- 
cated in figure 10, since this quantity determines the endurance under 
emergency conditions. The emergency flight duration is discussed later. 

A small amount of fuel (l percent of gross weight) to be used during 
takeoff and climb is allowed for and is shown in the figures (Wp^ imb) • 

This fuel is not part of the reactor shield; therefore it is accounted 
for apart from the emergency fuel supply. 


Effect of Design Flight Conditions and Landing Speed 

As mentioned previously, it was decided that the airplane should be 
designed for flight at a Mach number greater than 0.70 and at an arbi- 
trarily assigned altitude of 30,000 feet. Some preliminary calculations 
indicated that an airplane of 400,000-pound gross weight would require 
approximately 100 megawatts of reactor power for these flight conditions. 
Accordingly, reactor calculations and design were based on this power 
output. More refined airplane performance calculations showed that the 
required reactor power would be 98.5 megawatts at a flight Mach number 
of 0.72. 
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The effects on payload of variations in both the design Mach number 
and altitude (each independently) were studied. An intrinsic factor in 
these variations is the constraint placed on the design wing area by the 
allowed landing speed and lift coefficient at landing. For the typical 
airplane, the landing speed was set at 175 feet per second, and the lift 
coefficient at landing was 1.75. These two values remained constant 
during the studies of the effects of design altitude and Mach number. 

In fact, the landing-speed requirement instead of the flight condition 
established the wing area for most of the cases studied. The effects of 
varying landing speed are also presented. 

Design flight Mach number . - The effects of variations in design 
flight Mach number on the component and payload weights are shown in fig- 
ure 10. There is a decrease in payload from 93,000 to 8000 pounds as 
the Mach number is increased from 0.60 to 0.85. It appears that the 
greatest payload could be carried at a design Mach number somewhat less 
than 0.60. At the design Mach number for the typical airplane (0.72), 
the payload is about 63,000 pounds. 

The decrease in payload is due to an increase in powerplant and 
reactor and shield weights. These weights depend on the required power, 
which increases from 74 to 138 megawatts as the Mach number varies from 
0.60 to 0.85. The increase in required power can be traced to two 
causes. First, even with a constant drag, the power required varies 
directly with flight speed. Second, the airplane drag increases with 
increasing Mach number, because the wing area cannot be decreased. Wing 
area is specified by the landing conditions and not by flight Mach num- 
ber in these cases. Consequently, at 30,000 feet and in the Mach number 
range considered, the lift-drag ratio is never optimum. Lift-drag ratio 
decreases from 22.8 to 15.7 as the Mach number changes from 0.60 to 
0.85. At a higher altitude the situation would be modified, because an 
increase in altitude requires an increase in the wing area for optimum 
lift-drag ratio. 

Design flight altitude . - The effect of changing design altitude, 
at a flight Mach number of 0.72, on the weights of the airplane compo- 
nents is shown in figure 11. The large variation in powerplant weight 
with changes in altitude is a result of several effects. The maximum 
lift-drag ratio (minimum thrust and power) occurs at an altitude of 
approximately 40,000 feet. However, because of higher air density at 
lower altitude, the engine weight is a minimum at 30,000 feet. Minimum 
heat-exchanger weight occurs at 35,000 feet. As a result, the total 
weight of the powerplant (which consists of engines, nacelles, heat ex- 
changers, manifolds, pumps, and helium lines) is minimum at 35,000 feet. 
The lowest reactor and shield weight occurs at the altitude for minimum 
reactor power (40,000 ft) in spite of the reduction of air- scattered 
dose rate at higher altitudes. A maxi mum payload of nearly 70,000 pounds 
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occurs at an altitude of just above 35,000 feet compared with a payload 
of 63,000 pounds for the 30,000-foot design altitude selected for the 
typical airplane. 

The optimum altitude (maximum payload) depends, of course, on the 
design flight speed. Maximum payload would occur at a lower design 
altitude if a lower design speed were selected. Performance for simul- 
taneous variations in design speed and altitude, was not determined. 

Design landing speed . - Figure 12 shows the variation with design 
landing speed of the weights of the airplane components for a design 
altitude of 30,000 feet and a design Mach number of 0.72. The range of 
landing speeds presented varies from 140 to 250 feet per second. Design 
landing speed for the typical airplane was 175 feet per second. 

The payload increases from 25,000 pounds at a design landing speed 
of 140 feet per second to a maximum of 82,000 pounds at a design landing 
speed of 230 feet per second. For the range of landing speeds shown in 
figure 12, wing area depends on landing speed and not the flight condi- 
tions. Therefore, structure weight decreases with increasing landing 
speed. The maximum design lift-drag ratio occurs at a landing speed of 
220 feet per second; and, at the same speed, the weights of the power- 
plant, reactor, and shield are a minimum. 

For a design landing speed of 250 feet per second or more, the wing 
area is determined by the design flight conditions. Therefore, no fur- 
ther variations in the component weights with landing speed will occur 
for speeds above 250 feet per second. 


Effect of Airplane Gross Weight 

The effects of changing the airplane design gross weight on the 
component weights are shown in figure 13. Flight conditions were taken 
at the typical values of Mach number (0.72) and altitude (30,000 ft). 

As might be expected, the sum of the powerplant and structure weights 
is nearly a constant fraction of gross weight. At a gross weight of 
300,000 pounds, the structure, equipment, and powerplant account for 53.6 
percent of the gross weight. These items represent 55.2 percent of an 
800,000-pound airplane. Specifically, the powerplant weight decreases 
from 21 to 20 percent, the equipment weight decreases from 6.3 to 4.2 
percent, and the structure weight increases from 26 to 31 percent of 
gross weight as the gross weight increases from 300,000 to 800,000 pounds. 

As a consequence of these variations, the sum of the reactor, shield, 
and payload weights is also nearly a constant fraction (45 percent) of 
the gross weight. The weights of these items increase from 136,000 to 
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348.000 pounds for the range of airplane gross weights shown in figure 
13. Compared with this large variation, the reactor and shield weight 
varies only a small amount with gross weight. Of course, the reactor 
and shield weights depend on the power and on the crew-reactor separa- 
tion distance, which are related to gross weight. 

The methods used for calculating reactor and shield weights as a 
function of power and of crew-reactor separation distance are presented 
in detail in appendix J. Very briefly, the reactor core length was 
maintained constant, the core frontal area was assumed to he propor- 
tional to the power, and reflector and thermal- shield thicknesses were 
held constant. Based on these assumptions, the fuel-element - moderator- 
block geometry will not depend on reactor power, and the number of fuel 
elements will vary directly with power. In addition, it is expected 

o-zcr 

that the amount of required for each fuel element will not vary 

significantly with power; however, no calculations were made to deter- 
mine fuel loadings for the large-size reactors. 

An increase in gross weight from 300,000 to 800,000 pounds requires 
a change in reactor power from 86.5 to 182.4 megawatts. The associated 
increase in reactor and shield weight is from 110,000 to 146,000 pounds, 
which is a decrease in the fraction of gross weight from 37 to 18 per- 
cent. The net result is that the payload increases from 28,000 to over 

200.000 pounds for the gross-weight variation considered. However, the 
percentage of gross weight available for payload changes from 9.5 to 
over 25 percent. This result illustrates a fundamental characteristic 
of nuclear-powered aircraft, that a large increase in payload-carrying 
efficiency results from an increase in gross weight of the airplane. 

In the selection of design gross weight for a logistic airplane, 
some consideration should be given to the growth potential of the air- 
plane. The gross weight of many airplanes designed in the past has 
grown by as much as 30 percent of the initial weight. The growth in 
gross weight has usually been accomplished without increasing the struc- 
tural weight of the airplane, through redesign and refinement of the 
structure. Therefore, as the airplane structure is developed, an In- 
crease in gross weight results In a large increase in payload-carrying 
ability. Such improvements are greater than the variations of payload 
with design gross weight shown In figure 13. 

It is expected that a similar growth potential exists for nuclear- 
powered aircraft. Therefore, in order to take advantage of the improve- 
ments mentioned, the design gross weight of the typical airplane was 6et 
at 400,000 pounds. At this weight the payload-carrying efficiency is 
low compared with what could be achieved at the higher gross weights. 
However, with past experience as an indication, the 400,000-pound air- 
plane would grow in gross weight and payload-carrying ability. The 
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gross weight of the developed, airplane would not be excessive for cur- 
rently available manufacturing and runway facilities. 

The improvements due to structural refinements leading to increased 
gross weight and payload of particular airplanes in the past have led to 
an erroneous belief that the same effect could be achieved simply by 
increasing the design gross weight of the airplane. However, structural 
improvement s are available for any of the initial design gross weights . 

If the airplane gross weight is to be increased after design, ade- 
quate thrust must be provided for flight with the increased gross weight. 
Adequate thrust can be provided by designing oversized engines or by 
operation at other than design flight conditions. In the present typ- 
ical airplane, added thrust is available at flight speeds below the de- 
sign Mach number of 0.72. Since the design wing area of the typical 
airplane is larger than that required for optimum lift-drag ratio, 
flight at the design altitude of 30,000 feet k could be maintained with 
an increase in gross weight and some decrease in flight Mach number. 


Effect of Design Aspect Ratio 

The effects of design aspect ratio on the airplane component weights 
are shown in figure 14 for a range of aspect ratios from 6 to 15. Design 
aspect ratio for the typical airplane is 10. 

The payload maximizes at a design aspect ratio of 8 for this air- 
plane; however, only a small variation of payload with aspect ratio is 
indicated. The variation in payload is due to several factors. Struc- 
tural weight increases with increasing aspect ratio. The powerplant and 
reactor and shield weights decrease with increasing aspect ratio because 
an increase in aspect ratio improves the lift-drag ratio and effects a 
decrease in required power. 

Since aspect ratio has a relatively small effect on payload and 
since the heat exchangers can be located in the wings as shown in fig- 
ure 2, the final choice of aspect ratio may depend on the space required 
for the heat exchanger. 


Effect of Crew Dose Rate 

For the typical airplane the total weight of reactor, shield, and 
payload is 183,000 pounds. A decrease in crew dose rate results in a 
heavier shield and a correspondingly lighter payload. Figure 15 shows 
the change in shield and payload weights with a change in crew dose 
rate. The shield weights shown were calculated according to the methods 
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given in appendix J. No claim is made for the accuracy of these calcu- 
lations; the trends, however, are believed correct. A discussion of the 
calculations is given later. 

A tenfold change in crew dose rate from the nominal value of 0.025 
rem per hour causes a change in shield weight, and therefore payload 
weight, of about 35,000 pounds. A decrease in dose rate to the currently 
allowable laboratory rate of 0.3 rem per week (40-hr week) or 0.0075 rem 
per hour would require a decrease in payload of 17,000 pounds. 


ENGINE 

As indicated in the INTRODUCTION, turboprop engines were chosen for 
the airplane. This type of engine is well suited for driving the typ- 
ical airplane, which is designed for 0.72 Mach number at 30,000 feet. 

The section GENERAL REQUIREMENTS has explained the selection of eight 
engines to drive the 400,000-pound airplane. A shaft power of 5900 
horsepower per engine approximates that obtainable from turboprop en- 
gines in current production in the United States (e.g., Pratt & Whitney 
T34 and Allison T40) . 

In line with the conservative mission selected for the first nuclear 
aircraft, the engine design was also made conservative. As a conse- 
quence, it was felt unnecessary to make detailed designs of the engine 
components, since the magnitude of their problems should not be great. 
This section will, however, discuss some of the parameters of the com- 
pressor and turbine designs. 


Selection of Engine Design Parameters 

A range of theoretical turboprop engine designs is presented in ref- 
erence 3. For flight conditions approximating those of the typical air- 
plane, this reference indicates that a compressor pressure ratio of 11 
and a turbine pressure ratio of 9 would yield low specific fuel consump- 
tion in a chemically powered engine. Low specific fuel consumption is 
indicative of high propulsive efficiency. Therefore, it was reasonably 
expected that a combination of design parameters that result in low spe- 
cific fuel consumption would likewise result in a requirement of low 
reactor power. 

Preliminary design-point calculations were then made for a compres- 
sor pressure ratio of 11 and also for 8 and 15. For the latter two val- 
ues, turbine pressure ratios of 6.7 and 12 were used, respectively. A 
turbine-inlet temperature of 1800° R was specified. Figure 16 shows 
that, although net jet thrust per unit airflow F/w a decreases, airflow 
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per unit turbine frontal area vj and shaft horsepower per unit air- 
flow HP/w a both rise with increasing compressor pressure ratio. 

Higher propulsive efficiency is a further advantage of the higher com- 
pressor pressure ratios. On the other hand, pressure ratios greater 
than about 12 are not considered feasible for a one-spool compressor. 

On the basis of these considerations, a one-spool compressor with a 
pressure ratio of 11 was selected for design of the engine of the typ- 
ical airplane. Subsequent calculations indicate that a value of com- 
pressor pressure ratio of 11 yields max imum payload. 

Figure 17 was also used in the preliminary design in selecting the 
turbine-inlet temperature. This figure shows that F/w a and Wg/Aip 

decrease while HP/w a rises with increasing turbine-inlet temperature. 

A turbine-inlet temperature of 1800° R was chosen for the typical air- 
plane. Subsequent calculations, discussed later in this section, show 
that maximum payload is obtainable at 1900° R temperature. 

The performance of an experimental eight-stage compressor with a 
pressure ratio of 10 is presented in reference 4. This reference shows 
that such a high pressure ratio is obtainable in a one-spool machine 
that will operate with an adiabatic efficiency of 0.84. This value cor- 
responds to a polytropic efficiency of 0.88 at the pressure ratio of 10. 
With an equivalent tip speed of 1100 feet per second and an inlet hub- 
tip radius ratio of 0.48, this experimental compressor had a weight- flow 
capacity of about 30 pounds per second per square foot of frontal area. 
The corresponding values for the compressor design of the present report 
are 1000 feet per second equivalent tip speed, 0.5 hub-tip radius ratio, 
and 20.9 pounds per second per square foot equivalent weight flow. The 
compressor tip diameter is 29.4 inches. Eleven conservative compressor 
stages are required for the pressure ratio of 11. Compressor and tur- 
bine design parameters axe summarized in table III. 

The turbine is designed for 30,000-psi centrifugal stress at the 
rotor-exit hub radius. This is a conservative value for a turboprop 
engine because of the low stagnation temperature at the rotor exit 
(1103° R for the typical airplane). Reference 5 indicates that the life 
of the turbine blade metal should well exceed 1000 hours. 

Use of a design turbine-inlet relative Mach number of 0.6 at the 
hub radius, 0.6 exit hub-tip radius ratio, the required enthalpy drop of 
183.1 Btu per pound, and the tip speed corresponding to 30,000-psi stress 
result in a requirement of three turbine stages. This determination was 
obtained with the aid of reference 6. With such conservative design, 
there should be no difficulty in realizing a turbine polytropic effi- 
ciency of 0.85. Turbine tip diameter is 34.9 inches. 
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The methods of calculating design-point thrust, weight flow, horse- 
power, and engine component weights are given in appendix B. 


Effect of Engine Design Parameters on Component Weights 

Effect of turbine-inlet temperature . - The effect of turbine- inlet 
temperature on airplane component weights is presented in figure 18(a) . 
Design-point flight Mach number, altitude, and compressor pressure ratio 
were maintained at 0.72, 30,000 feet, and 11, respectively. Of prime 
interest on this plot is the increase in payload with temperature up to 
about 1900° R, above which the payload decreases. For the conditions of 
figure 18(a), reactor power decreases with rising turbine-inlet temper- 
ature because of improved cycle efficiency. As shown in figure 18(a), 
this results in decreasing reactor plus shield weight by about 10 per- 
cent over the range studied. 

Additional explanation for the payload variation with turbine-inlet 
temperature is provided by the variation in powerplant weight. Figure 
18(a) shows that powerplant weight decreases to a minimum, then increases 
as turbine-inlet temperature rises. A breakdown of powerplant component 
weight variation with temperature as presented in figure 18(b) is in- 
structive. Total powerplant weight, according to this figure, is a min- 
imum at about 1880° R. The weight of the engines decreases with rising 
turbine-inlet temperature because of the decrease in airflow require- 
ments. The increase in heat-exchanger weight with rising temperature 
has a major effect on powerplant weight. With constant temperature of 
the helium entering the heat exchanger, in order for turbine-inlet tem- 
perature to rise, additional surface area must be provided in the heat 
exchanger. Furthermore, the resulting higher heat-exchanger wall tem- 
peratures as turbine-inlet temperature rises pose a more serious stress 
problem. The walls must therefore be thicker. These two factors cause 
heat-exchanger weight to rise continuously with increasing turbine-inlet 
temperature . 

Figure 18(b) thus shows that, up to about 1880° R, the decrease in 
engine weight offsets rising heat-exchanger weight so that powerplant 
weight decreases. Above this temperature, the heat exchangers become so 
heavy that they cause powerplant weight to rise. 

Returning to figure 18(a), it will be noticed that structure and 
equipment weight increases to a maximum at approximately the temperature 
at which powerplant weight is a minimum. This is to be expected, for as 
powerplant weight is decreased, the distributed lift forces on the air- 
plane wings cause greater bending moments at the roots of the wings. 

The structure must be made stronger and thus heavier. 
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The net result of the variations in these airplane component weights 
with turbine- inlet temperature is that payload is a max imum at approx- 
imately 1900° R- At this temperature the payload is 6000 pounds greater 
than that at the 1800° R temperature used in the final design. 

Effect of compressor pressure ratio . - Figure 19(a) presents the 
effect of varying compressor pressure ratio on airplane component 
weights. Turbine-inlet temperature remains at 1800° R and flight con- 
ditions remain at their design values in this figure. Payload is essen- 
tially constant at compressor pressure ratios of 8 and higher. For the 
conditions of this figure, it can be shown that reactor power decreases 
as compressor pressure ratio rises because of increases in cycle effi- 
ciency. Figure 20 shows that this results in decreasing weight of re- 
actor plus shield. 

Figure 19(a) also shows that powerplant weight decreases to a min- 
imum and then rises as compressor pressure ratio is increased. Varia- 
tion in powerplant component weights with compressor pressure ratio is 
shown in figure 19(b). For the conditions of this figure, airflow first 
decreases and then increases with rising compressor pressure ratio. 

This accounts for the minimum in the curve of engine weight. 

In figure 19(b) heat-exchanger weight decreases with rising compres- 
sor pressure ratio up to a value of about 9.5. Further increase in com- 
pressor pressure ratio results in heavier heat exchangers. This effect 
is discussed in the section on heat exchangers. 

Reference to figure 19(a) shows that the weight of structure and 
equipment is essentially insensitive to change in compressor pressure 
ratio. The previously discussed variations in reactor plus shield and 
powerplant weights result in little variation in payload above a com- 
pressor pressure ratio of about 8. Maximum payload occurs at a pressure 
ratio of approximately 11. 


HEAT EXCHANGERS, PUMPS, AND LINES 

In order to deliver the energy developed in the reactor to the tur- 
boprop engines, a helium-to-air heat exchanger is required, as well as 
helium lines, air ducts, and pumps for the helium. Figure 3 shows the 
arrangement of these components. Each heat exchanger is located close 
to its engine so that the length of relatively large -diameter air ducts 
can be made as small as possible. The helium pumping power is extracted 
from the engine shafts, and there is one helium pump for each engine. 

The primary reason for separating the heat exchanger from the en- 
gine was to permit engine handling and maintenance without any disruption 
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of the helium system. In addition, calculations showed that the heat 
exchangers can he designed to he stowed in the wings as shown in 
figure 3. 


Heat Exchanger 

Materials considerations . - As mentioned previously, the turbine- 
inlet temperature and helium temperature at the reactor outlet for the 
typical airplane were set largely hy materials limitations in the heat 
exchanger. Strength and oxidation resistance at high temperatures are 
required of the material used in the heat exchanger. Some nickel-hase 
alloys are suitable for this application. Figure 21 shows 1000-hour 
stress -rupture data for a nickel-hase alloy. The stress values assumed 
for the various heat-exchanger designs are also shown in figure 21 at 
metal temperatures corresponding to the peak temperature in the heat 
exchanger. Corresponding turbine-inlet temperatures are also indicated 
(refer to discussion of effects of turbine-inlet temperature in ENGINE 
section) . 

Type of heat exchanger . - A shell-and-tube heat exchanger with the 
helium contained in the tubes was selected for this airplane. On the 
basis of size and weight considerations, it was also decided to use 
multiple passes of the helium tubes with counter crossflow of the air- 
stream. Calculations showed that a four-pass heat exchanger (including 
shells and structure) was somewhat lighter than a three-pass type and 
that the resulting size allowed stowage of the heat exchanger in the 
wings. Circular tubes were chosen for reasons of mechanical strength. 

There are many other items to consider in the design of the heat 
exchanger. Some of these are the size of the tubes, the spacing and 
arrangement of the tubes, the desirability of finning, the fin material, 
fin spacing and thickness, and allowable pressure drop in each fluid or 
gas. In addition, theoretical or empirical relations for the heat- 
transfer coefficients are required. The calculation of heat-transfer 
coefficients for flow inside tubes is well established (e.g., refs. 7 
and 8) ; however, the situation is much different for flows across banks 
of tubes. It appears that the dependence of heat-transfer coefficient 
and pressure losses on factors such as tube and fin geometry has not 
been established for all cases. Therefore, it was decided to limit the 
heat exchangers considered to those specifically given in the compila- 
tion of reference 9. Such a decision restricts the study of heat ex- 
changers and prohibits a complete optimization of all variables; however, 
through the use of empirical heat-transfer and pressure-loss data, the 
heat exchangers can be accurately designed. Design calculations were 
done on an IBM 653 computer, as outlined in appendix C. 
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Finned- and plain-tube heat exchangers were designed for identical 
conditions in order to determine the value of finning in the applica- 
tion. Figures 48 and 100 of reference 9 provided the geometry and heat- 
transfer characteristics for the plain- tube and the finned- tube heat 
exchangers, respectively. It was assumed that the fins and tubes were 
made from the same alloy. Although the sizes of these two types of heat 
exchanger were not equal, each type would fit into the wing. It was 
found that the finned heat exchanger was heavier than the plain-tube 
type; consequently, the plain- tube he at -exchanger was selected. The 
simplicity of fabricating a plain-tube heat-exchanger appeared to be an 
advantage, also. Of course, the use of fins might have proved worth- 
while with other fin-and-tube geometries or with fins constructed of a 
material having a high thermal conductivity such as copper. Copper fins 
would require oxidation protection for this application, however. 

The pressure losses of the air and helium are important in the heat- 
exchanger design, since the weights of many components are influenced by 
these quantities. Actually, the selection of optimum pressure losses 
for the heat-exchanger requires a complete analysis of the airplane sys- 
tem, because the helium and air pressure losses affect heat-exchanger 
weights, engine and nacelle weights, and reactor and shield weights. 
Analysis showed that very low pressure loss was required in the helium 
flow to keep pumping power to a reasonable level. For the typical air- 
plane, 834 horsepower per engine were required to pump the helium with 
a pressure ratio of 1.085. The corresponding helium pressure ratio 
across the heat exchanger (outlet to inlet) is 0.99, a value that is 
probably near the optimum for this system. Efforts were made to mini- 
mize the pressure losses on the air side of the heat exchanger. The 
analysis showed that the desirable value of outlet-to-inlet pressure 
ratio for the heat exchanger is approximately 0.92. The design ratio 
of turbine-inlet pressure to compressor-outlet pressure is 0.85 includ- 
ing pressure losses in the air ducts and burner. 

Heat exchanger for typical airplane . - The heat exchanger for the 
typical airplane was designed according to the geometry and data given 
in figure 48 of reference 9. The following are some of the features of 
the heat exchanger that meets the design specifications given in table 
I. Approximately 530 tubes of 0.32-inch inside diameter are required 
for each heat exchanger. For an allowed stress of 7000 psi, the outside 
diameter of the tubes must be 0.372 Inch. The peak metal temperature is 
estimated to be 1690° F with a turbine-inlet temperature of 1340° F 
(l800° R) • The total length of each tube is 30.3 feet, and the dimen- 
sions of each heat exchanger are 7.6 feet long, 1.4 feet high, and 3.3 
feet wide. (The width is measured in direction of airflow in fig. 3.) 

It is estimated that the total weight of each heat exchanger will be 
2000 pounds, 1500 pounds for the tubes and 500 pounds for shell, tube 
supports, structure, and headers. Figure 22 shows the arrangement of 
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the helium headers on the heat exchanger. The helium that leaves the 
pump is ducted around the high-temperature header mainly to reduce the 
pressure differential from the hot metal and thereby reduce greatly the 
header weights. 

In this design, pressure loads are carried by the cool outer jacket. 
Some additional discussion of the header design is given in appendix C. 

Effect of compressor pressure ratio . - The effect on airplane com- 
ponent weights of varying design compressor pressure ratio has been dis- 
cussed previously, but some additional information on the heat exchanger 
is required. Furthermore, a discussion of the effects on the heat ex- 
changer of varying pressure ratio will serve to illustrate the analyses 
of heat exchangers that were made for all design variations considered 
in the AIRPLANE and ENGINE sections of this report. 

The variations of some of the important he at- exchanger parameters 
with compressor pressure ratio are shown in figure 23 for a reference 
airflow rate of 100 pounds per second. Helium inlet and outlet temper- 
atures and turbine-inlet temperature were constant and equal to the val- 
ues assigned for the typical airplane. Constant percentage pressure 
losses were also assumed. Figure 23 shows that the tubing weight, which 
is directly related to the heat-transfer surface area, minimizes at a 
pressure ratio of 10.5. The variation in weight (surface area) results 
from a combination of effects. The increase in inlet air density with 
pressure ratio provides higher heat-transfer coefficients and lower 
weight. However, the increase in inlet air temperature associated with 
increasing pressure ratio requires increased effectiveness of the heat 
exchanger, which in turn means larger surface area and weight. The ob- 
served decrease in air inlet Mach number with increasing pressure ratio 
also increases the surface area and weight through a lowering of the 
heat-transfer coefficient. The inlet Mach number variation is required 
to maintain constant percentage pressure losses in view of the varia- 
tion in the airflow length (width of heat exchanger) . 

All heat-exchanger designs were calculated on the basis of 100 
pounds per second of airflow and then were properly scaled for the cor- 
rect engine airflow. The number of tubes, the weight of tubes, and the 
depth measured perpendicularly to the airflow all vary directly with 
airflow . 


Helium Pumps 

Since the helium pumping power is quite large (834 hp of the 3410 
hp developed by each engine at design flight condition) , pump efficiency 
is important in this cycle, because a decrease in pump efficiency will 
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require an increase in engine size and weight and reactor and shield 
weight. In view of the importance of pump efficiency, an axial-flow 
pump was chosen over the centrifugal type to take advantage of the poten- 
tially high efficiencies of axial- flow compressors. Furthermore, since 
the pumps are very small and low in weight, any weight advantages of 
either type of pump are unimportant. 

The design of the axial- flow pump is given in appendix D along with 
an estimation of the pump performance map. Some of the features of this 
pump are as follows. Three identical stages can produce a pressure ratio 
of 1.085 with an estimated punp adiatatic efficiency of 0.80. The design 
tip speed is 1500 feet per second, the tip diameter is 5.23 inches , and 
the rotational speed is 65,800 rpm. Inlet relative Mach numbers are less 
than 0.4. The aerodynamic design is straightforward; conventional tech- 
niques and "blading can "be used. However, the high rotational speeds and 
high temperature of the helium (800° F) will pose problems in bearing 
and shaft seal design. 

The bearing problem can be solved through the use of small, high- 
quality ball bearings, provided that they are well lubricated and in- 
stalled carefully. The shaft sealing problem requires more attention. 
Preliminary calculations showed that labyrinth seals would allow exces- 
sive leakage of helium. Face-type seals will keep leakage to a low 
value; however, the combined effects of high rubbing speed, temperature, 
and pressure differential make their design difficult. A possible solu- 
tion to the sealing and bearing problem is shown in figure 24. This 
proposal utilizes lubricating oil at a pressure of 600 pounds per square 
inch to reduce the pressure differential across the helium seal and to 
provide some cooling. The arrangement shown In figure 24 appeared rea- 
sonable to a manufacturer of face-type seals. 


Helium Lines 

Two problems are associated with the design of the helium lines. 
Sufficient mechanical strength must be provided (l) to hold the high- 
pressure helium and (2) to prevent destructive heating of the airplane 
structure. Both of these problems can be solved if concentric lines are 
used with the inner line containing the high- temperature (1790° F) 
helium. The cool helium (790° F) flows in the annulus and removes the 
high-pressure loads from the hot metal parts. As a result, the inner 
liner of the concentric ducts can be made thinner. The weight of a con- 
centric line will be approximately one-fifth the weight of two single 
lines . 

Heat transfer between the two helium streams in the concentric 
lines was investigated, and it was found that the temperature of the 
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helium decreased approximately 50° F in transit from the reactor to the 
outboard engines if the inner pipe was uninsulated. However, the heat 
loss from the high-temperature helium can he reduced significantly by 
using a lightweight, three-layer construction for the walls of the hot 
pipe. Typically, the hot pipe can be built from a tube having corrugated 
walls placed between two smooth-walled tubes as shown in figure 4, The 
insulating effect of this construction reduces the heat transfer to the 
extent that the temperature of the helium decreases only 2° P in flowing 
from the reactor to the outboard engine. 


REACTOR 

Considerations Leading to Selection of Reactor 

The reactor designed for the typical airplane has already been 
described briefly. The reactor layout is shown in figures 5 to 7. Many 
considerations enter into the selection of the reactor type to use, some 
of which will be discussed in this section. 

As mentioned previously, this study was made using a helium-cooled 
reactor. Helium has several advantages: 

(1) Has good heat-transfer characteristics 

(2) Does not become radioactive by neutron bombardment, obviating 
an intermediate heat exchanger and shielding of external 
circuitry 

(3) Is chemically inert; consequently the materials it contacts 
will suffer no corrosion nor mass-transfer difficulties 

Reference 10 points out that "the only contamination in the helium sys- 
tem would be due to the leakage of radioactive materials from the reactor 
into the helium stream." Helium has the disadvantage that it diffuses 
readily through imperfections in the containment system; hence a great 
deal of care in fabrication is required. 

The use of the inert gas, helium, leads to the use of high- 
temperature refractory materials such as molybdenum, columbium, and so 
forth, that cannot be used with air without cladding because of oxida- 
tion and corrosion effects. The strengths of these materials are gener- 
ally adequate for reactor design up to temperatures of approximately 
3500° R. A solid-core reactor was chosen for the present study to take 
advantage of the characteristics of these materials. The limiting tem- 
perature in the reactor - he at -exchanger system with helium as the re- 
actor coolant will be the heat-exchanger wall temperature. Materials 
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must be used in the heat exchanger that are compatible with the air from 
the engine, which is on one side of the heat-exchanger walls. 

A helium temperature of 2250° R out of the reactor was set as the 
design value from considerations of the he at -exchanger limitations. A 
temperature of 1250° R into the reactor was determined from considera- 
tion of heat-exchanger size and the transfer of the heat to the engine 
air. The helium is required to be at high pressure to reduce the void 
in the reactor as well as to increase the heat-transfer coefficient. A 
high pressure also makes the heat exchanger and helium lines more amen- 
able to fabrication with respect to size. A value of 1250 pounds per 
square inch into the reactor was used. 

The high pressure of the gas in the reactor leads to the problem of 
containing the gas. A pressure shell was chosen for the following 
reasons : 

(1) It avoids use of headers for the fuel tubes with large pressure 
differentials across them. 

(2) It removes pressure loads from high -temperature (fuel element) 
region. 

(3) A high-density gamma shield is needed anyway. 

Thermal shields, located between the reactor core and the pressure shell, 
were required to reduce the rate of heat generation and consequently the 
thermal stresses in the pressure shell. A cylindrical geometry was 
chosen for the core shape. One benefit of the cylindrical geometry is 
that all fuel-element tubes are the same length. 

A small reactor size permits the design of a low shield weight, an 
item that can make or break nuclear airplane feasibility. To obtain 
small reactors, the fuel must be highly enriched. The performance gain 
due to the relatively small reactor will be of enough advantage to off- 
set the extra cost of the fuel. For the present reactor the fuel chosen 

was 92.5 percent and 7.5 percent U 


Design Procedure 

The main steps in the design of the reactor are as follows : From 

the known values of reactor-coolant inlet and outlet temperatures and 
reactor power distribution and assumed values of fuel-element length, 
moderator-to-reactor volume ratio, and heat generation, 

(a) Determine the fuel-element and moderator-block dimensions to 
meet specified temperature and coolant pressure-drop limitations. From 
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these calculations the ratio of moderator cross-sectional area to frontal 
area and power per fuel element are obtained. 

(b) From the power of the reactor and the dimensions determined in 
the first step, determine the number of fuel elements and get the reactor 
diameter. 

(c) For an estimated excess reactivity, from nuclear calculations 
determine the fuel quantity required for the reactor obtained in steps 
(a) and (b) . Excess reactivity is built in because of the effects of 
poisons and burnup. 

The first calculations made along the foregoing lines are of a pre- 
liminary nature, because the structure and control rods are omitted for 
simplicity. The calculations are repeated with other fuel-element and 
moderator dimensions, as required, until a satisfactory design is 
achieved. Then the nuclear calculations are continued with these final 
dimensions and with the items that were omitted. 

After the core and reflector aspects of the reactor are finalized, 
the thermal shields, pressure shell, supporting structure, and control- 
rod mechanisms are designed. The process requires much iteration, as 
can be inferred from the preceding discussion. 


Fuel Element 

A solid-core reactor can be constructed from parts made in numerous 
ways and of differing shapes, sizes, and materials. The choices made 
and some reasons for the choices of the parts used in the reactor are 
given in this section on the fuel element and following sections of the 
report . 

Materials . - The use of helium leads to possibilities of using 
molybdenum or other refractory metals in the fuel elements. To minimize 
radiation damage in fuel elements, the fuel is contained in a matrix of 
nonfissile material, and superposed on the matrix is cladding consisting 
of only the nonfissile material. Matrix materials with melting points 
above 4000° F are about the only ones of interest in an application such 
as the one considered herein. Even though the outlet helium temperature 
is only 1800° F, the maximum fuel-element temperature will be much higher. 
If the refractory materials that appear to be scarce and expensive are 
omitted, the metals that are left for consideration are tungsten, tanta- 
lum, molybdenum, and columbium. Considerations of the strength, ductil- 
ity, ease of fabrication, thermal conductivity, neutron absorption cross 
section, and state of development of these four materials led to the 
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choice of molybdenum as the best material to use for the fuel elements. 
This choice was mainly based on the state of development of molybdenum 
compared with the other materials. The choice does not mean that molyb- 
denum is superior in all the properties mentioned. Consideration of the 
materials at a later date might lead to another choice, but this would 
not greatly affect the results . 

The dispersant fuel phase (or the uranium compound) of a fuel ele- 
ment should be chemically stable with regard to its environment at the 
operating temperature. It should contain the highest possible density 
of uranium per unit volume. Among possible fuel dispersants, uranium 
dioxide U0 2 and uranium carbide UC or UCg are the most promising. Some 
properties are as follows : 


Compound 

(a) 

Melting 

point, 

°F 

U 

1.0 

2072 

UC 

.69 

4132 

uc 2 

.56 

4352 

U0 2 

.53 

4532 


Cranium content per unit 
volume of dispersed phase 
relative to density of 
uranium. 


The table shows that the several uranium compounds permit much higher 
operating temperatures than does natural uranium. Uranium dioxide stands 
out because its technology and handling are well understood, although all 
the compounds listed have relatively high density of uranium, UC being 
better than U0 2 in this respect. It was decided to use UO2 in the reac- 
tor fuel elements, and calculations were made using this material. Later 
considerations caused a change to UC. These considerations will be taken 
up later in the report. 

Shape . - In choosing the fuel-element shape, considering the most 
general types used (i.e., solid cylindrical rods, hollow cylinders, and 
plate type), the hollow cylinder or tube type was finally chosen. The 
basis of this choice is that it is strong and that it can be cooled by 
forced convection by passing coolant through the hole in the element as 
well as through an annulus formed by the outer radius of the fuel element 
and the Inner radius of a moderator block. In order to get a large 
surface- to-volume ratio for good heat transfer, the element must be small. 
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Analytical methods used in design . - The analytical methods used 
for determining the temperatures and stresses in the fuel elements and 
the pressure drops in the coolant passages are given in detail in appen- 
dixes E and F. In the temperature analysis the fuel element was assumed 
to be composed entirely of molybdenum, because the amount of uranium in 
each element was unknown at the start of the calculations, and to include 
the thermal conductivity of UOg in the analysis would have complicated it 
unduly. The conductivity of UOg is much less than that of molybdenum, so 
that the results of the temperature analysis used herein would be opti- 
mistic for the fuel elements where large amounts of UC >2 are used. The 
power distribution assumed for the heat-transfer calculations at first 
was constant in a radial direction and "chopped" 2/3-cosine law in the 
axial direction (see appendix E). Consequently, the fuel elements near 
the reactor core periphery contain the most fuel, and the calculated tem- 
peratures would be more in error for these elements than for the elements 
near the center of the core. It was assumed that 90 percent of the re- 
actor power was generated in the fuel elements . 

The assumption was also made in the stress analysis that the fuel 
element was made of one material only, or the cladding and matrix can be 
represented by one ring. If the element had been divided into two clad- 
ding rings and one "meat" ring, each would have to be treated separately. 
Such a stress theory development was considered beyond the scope of this 
report. Consequently, stresses were calculated using the analysis and 
assuming that the fuel element was made entirely of molybdenum. In this 
way some idea of the upper limit of the strength of the element could be 
obtained. 

Tube design . - The fuel element used in the reactor is shown in fig- 
ure 25. It consists of a sandwich of three concentric tubes, the tubes 
being continuous through the core and the end reflectors as shown in fig- 
ure 6. The middle tube contains molybdenum plus the uranium compound 
through the core length; it is only molybdenum through the end reflector 
lengths. The two outer tubes form the cladding. The cross section shown 
in figure 25 is section D-D of figure 5. Also shown in this figure, be- 
sides the fuel element and moderator block, are a molybdenum spring 
spacer for positioning the fuel element within the moderator block and a 
molybdenum bushing on which the spacer bears rather than on the fairly 
brittle BeO moderator block. A number of these spacers and bushings are 
placed along the fuel-element length as shown in figure 6. The dimen- 
sions of the fuel element, moderator block, and other parts are included 
in figure 25. The inner coolant passage is 0.25 inch, the cladding 
thicknesses 0.007 inch, and the "meat" thickness about 0.030 inch. The 
outer coolant passage, the annulus, has a spacing of 0.071 inch. 

The reactor was designed on the basis of 100 megawatts (95 Mw of 
heat in core, 5 Mw of heat in side reflector and thermal shields), as 
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stated previously. Because the coolant picks up heat from the thermal 
shields and side reflectors and then passes through the core, the total 
heat picked up By the coolant is 100 megawatts. On the basis of some 
simple preliminary calculations, it was decided to cool the fuel ele- 
ment on both sides with helium as shown in figure 25. The temperatures 
of the fuel element were calculated with an allowable pressure drop of 
approximately 2.5 percent of the inlet pressure (which was 1250 lb/sq 
in.), helium temperature in and out of the core of 1300° and 2250° R, 
respectively, core length of 24 inches, and end reflector thickness of 
4 inches. (There is a rise of 50° F in helium temperature through the 
thermal shields and side reflectors, which is the reason for the 1300° R 
temperature into the core. ) The helium flow rate to the reactor had been 
determined to be 76 pounds per second. The number of fuel elements fi- 
nally used was 640. Eighteen positions where fuel elements would be are 
occupied by through-bolts. Some of the coolant flows along the control 
rods and through-bolts. This flow was estimated, on the basis of obtain- 
ing the same temperature rise of the coolant or 950° P, to be roughly 2 
percent. The flow through each fuel-element - moderator-block assembly 
is then about 0.1143 pound per second. 

Axial temperature variations of coolant, meat, and cladding are 
given in figure 26. The pressure drop through the core is about 2.6 
percent, a value considered satisfactory with regard to the pressure 
drop allowable for the whole helium system, with flows of 0.0459 and 
0.0684 pound per second through the hole and annulus, respectively, of 
a fuel-element - moderator-block assembly. The max imum fuel-element 
temperature was 3235° R (fig. 26). The temperatures and pressure drop 
being considered satisfactory, the next step was to determine the 
stresses for this design. 

The temperatures given in figure 26 were then used to evaluate the 
properties in the stress equations for the case when the fuel element 
is considered to be entirely molybdenum. Built into the stress equation 
are the effects of the thermal gradients resulting from using a partic- 
ular material. These effects are separate from the property terms appear- 
ing in the equations. The sources of data for the properties of Mo re- 
quired in the calculation of the stresses are given in appendix F. Since 
relatively long life is required of the Mo, the stress-rupture strength 
was considered the best criterion to use. The allowable strength data 
for Mo at 1000-hour rupture life are shown in figure 27 and were extrap- 
olated from data for Mo plus 0.45 percent titanium in reference 11. 

These allowable data were used because adding the small amount of tita- 
nium increased the strength of Mo appreciably. 

The results of the stress calculations with only Mo in the fuel 
element are shown in figure 28, where the allowable and actual stresses 
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are plotted against the axial distance along the fuel element. For 
quite some distance along the fuel element it is apparent that the cal- 
culated actual stresses are greater than the allowable stresses. It 
should be pointed out that the theory in appendix F is for the elastic 
case and plastic flow is not included. It might be expected that plas- 
tic flow would relieve the stresses and the fuel element would be safe 
even in the regions indicated as unsafe in figure 28. It would be much 
better to have a fuel element that is safe along the entire length using 
elastic stress theory; then plastic flow would provide a further margin 
of safety. 

A means of changing both the allowable and calculated axial stress 
distribution is to vary the axial volumetric heat-source distribution, 
and thus the axial power distribution, in some manner different from the 
cosine distribution used to obtain the results shown thus far. Such 
a different heat-source distribution will result, of course, in a dif- 
ferent temperature distribution from that presented in figure 26. Con- 
sequently, calculations were made of the fuel-element temperatures, heat- 
source distribution, and stresses for a case that would result in the 
actual stress curve paralleling the allowable stress curve. The calcu- 
lations were obtained using methods given in appendixes E and F and from 
a simplified analysis not included herein. 

The axial variations in the fuel element and helium temperatures are 
shown in figure 29 for the case described previously. The axial varia- 
tion of the volumetric heat-source distribution for this case compared 
with the distribution for the 2/3 "chopped cosine" law is shown in fig- 
ure 30. The resulting stresses are shown in figure 31. The conditions 
used result in a "safe" fuel element from consideration of the calculated 
thermal stresses. The fact that the curves in figure 31 are not parallel 
may be attributed to the simplifying assumptions in the analysis. The 
distribution of fuel required to obtain the variation of volumetric heat 
source resulting in the stresses of figure 31 will be taken up in a later 
section of the report with the nuclear characteristics. 

No calculations were made considering the fuel element made up en- 
tirely of the uranium compound, because of the results described pre- 
viously. All efforts should be made to keep the volume of compound in 
an element as small as possible. Therefore, it was decided in the final 
design to change from UOg to UC on the basis of the table previously pre- 
sented. If some fuel elements do have large amounts of UC in them, the 
calculations presented will be rather meaningless with regard to the 
stresses, for reasons explained before. In any event, experiment alone 
■will determine whether the fuel elements are satisfactory. 
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Core and Reflector 

Preliminary nuclear characteristics . - For the application of nu- 
clear power to airplanes, it is important to design the reactor as small 
as possible so that the shield will not be overly heavy. Consequently, 
preliminary calculations of the uranium investment, fission spectrum, 
and flux for a series of reactors of varying core diameters and side- 
reflector thicknesses were made to determine how small a reactor seemed 
feasible. The smallest core diameter chosen was 17 inches. The control 
rods and structure were not included in the calculations. The analysis 
used for these nuclear calculations, and those to be presented later, 
are given in appendix G, with the exception of a few details that will 
be discussed in connection with some characteristics as they arise. The 
length-diameter ratio Z/d of the reactor core was chosen to be a little 
less than 1.0 for all cases because, on the basis of equations in refer- 
ence 12, it can be shown that this results in minimum amount of fuel 
needed for a required effective multiplication factor . 

The results showed that to keep the uranium investment from getting 
too high, the reactor (including side reflector) should be at least 38 
inches in diameter, as shown in figure 32. No values are shown on the 
figure ordinate on the amount of uranium because of the factors that 
were omitted. The trend of the curve is probably valid for reactors 
that include all the factors omitted. Another result was that the small 
reactors were in the high epithermal range, and in keeping with the idea 
of conservatism for the whole project this was not desirable. In addi- 
tion, there are doubts about the accuracy of the nuclear theory for such 
reactors, and the resulting calculations could be doubtful. Thus, even 
though smaller diameters would be desirable with respect to shield 
weight, these preliminary calculations indicated that the reactor should 
be greater than 38 inches in diameter (including reflector) on the basis 
of designing a more conservative reactor. 

Core -s true ture nuclear effects . - A few nuclear calculations were 
made after those just described in order to find out the effects of using 
different materials in the structure of the reactor. A calculation was 
made in which it was assumed that the structure would increase the 
frontal area of the reactor by 20 percent, and the multiplication factor 
was determined first assuming ZrC as the structural material and then 
Mo. The conditions used in the calculations are as follows: 
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Reactor power = 100 Mw 

Fuel amount = 70 lb uranium (92.5 percent 

Number of fuel elements = 528 

Pressure drop of helium through reactor = 4 percent 

Diameter of reactor * 32 in. 

Length of reactor = 24 in. 

No end reflector 

Side reflector = 4 in. 

Structure = 20 percent of total volume (increases frontal area) 

With ZrC structure, kg.^ = 1.06 

With Mo structure, kg-f-f = 0.74 

The result showed that ZrC would give a much higher multiplication factor, 
1.06 as compared with 0.74 for Mo, for the conditions used. Again from 
the standpoint of conservatism it was decided not to use a cermet for the 
structure, even though it gave a much better kg^f, because of possible 
brittleness. It was decided to use molybdenum for the same reasons given 
previously in connection with the choice of material for the fuel element. 
Columbium could probably be used equally as well at some future time when 
its development is advanced to the point where molybdenum is now. This 
would probably make considerable changes in reactor characteristics com- 
pared with molybdenum, but the over- all airplane results would probably 
remain unchanged. 

It should be pointed out that the percent structure in the final 
core was much less than that used in these calculations, so that the dif- 
ference between using a cermet like ZrC and a refractory metal will not 
be as great as that determined in the calculations. 

Final nuclear characteristics . - The foregoing calculations led to 
a final reactor with the characteristics presented in table IV and the 
general design features given in figure 6. The core radius was 17.36 
inches, and the core length 24 inches. To meet specified coolant 
pressure-drop limitations of the helium, 640 fuel elements were used. 

The result was a 2.6-percent pressure drop discussed previously in the 
fuel-element design. Both mechanical and control-rod structure (the 
sizes of the latter were roughly estimated) were included in the reactor 
to obtain the nuclear characteristics. In table TV, UC is the fuel. 
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From the calculations made without the structure, it was estimated 
that 130 pounds of uranium would be required to obtain the estimated re- 
quired k e ff and that there would be 10 pounds of fuel bumup in the 

estimated 1000 hours of operation proposed for the reactor. In making 
the calculations, it was estimated that the moderator and core structure 
would be operating at about 1400° F; and, since this structure is about 
90 percent of the reactor, the reactor core would not be operating at 
a temperature much higher than this. Since temperatures up to 1800° F, 
as compared with 1400° F, did not affect the cross sections appreciably, 
it was decided to use 1800° F in determining the thermal energy. Calcu- 
lations were made of k eff for the startup both cold and hot, and at 
1000 hours for a "clean" reactor and for one including poisons effects. 
The latter will hereafter be called "dirty" reactor. The calculation 
for the "clean" reactor at 1000 hours was made to determine the effect 
of fuel burnup alone on k e ^ . For the cold startup calculation the 
thermal macroscopic cross sections were determined at of 0.025 

electron volt. During the operation of the reactor there will be a 
buildup of xenon, samarium, and other fission product poisons which must 
be accounted for. They will reduce the kg^. To make the calculation 
in a rigorous manner, the nuclear calculations would be done over the 
entire operating time in steps corresponding to the "cycle" of operation 
from the start of a flight to the start of another flight. This would 
be very tedious and was thought to be unnecessary for this investigation. 
It was assumed that the reactor would be operating without shutdown over 
the entire period of operation of 1000 hours; this would result in a 
pessimistic value of kgff for the dirty reactor. In the calculation, 
the thermal macroscopic absorption cross sections of the fission product 

■npp Sni 

poisons 2 , , , of the xenon Z AC ., , and of the samarium Z D +v , were 

* a,th a,th a ^n 

added to the thermal macroscopic absorption cross section used for the 

clean reactor at 1000 hours. The microscopic fission product poisons 

cross section , used to get was obtained from reference 

a,tlr a,th' 

13. The cross sections Z^ e ., and 2?°+. were obtained by methods 

a,th a,hh 

given in reference 12 (pp. 332 and 338, respectively). The microscopic 
cross section o^ e ., used in Z? e + . was obtained from reference 14. 

a^uLL 

The resulting values of k e ff from these calculations are shown in 
figure 33. At startup, ^- e ff of 1.18 was obtained for the cold reactor 
and 1.12 for the hot reactor. After 1000 hours, the value of k e ff for 
the hot clean reactor was 1.11 and for the dirty reactor was 1.08. A 
dashed line is drawn between the hot condition at startup and the dirty 
condition at 1000 hours because, as stated before, unless the cycle of 
operation for a flight is decided, this curve cannot be drawn. These 
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results will be discussed in relation to their effect on the control 
problem in a later section on reactor control. The values of neutron 
flux obtained for the six energy groups used in the calculations for hot 
startup are shown in figures 34 and 35 for flux along a radius at the 
midplane of the reactor and for that along the axial centerline, respec- 
tively. In order to plot the fluxes of all six groups on one figure, 
some of the flux values were multiplied by a factor shown on the curves. 
For instance, the fluxes of group 1, cpj_, were multiplied by 20. The 20 
cp]_ were plotted on the ordinate. The resulting numbers are labeled 
relative flux cp r . 

The percentage fissions in the fourth, fifth, and thermal groups, 
the average flux of each group, and the energy range of each group are 
shown in table V for the hot startup and "clean" reactor at 1000 hours. 

The methods of obtaining these data and their use are discussed in 
appendix G. The code discussed in this appendix permits fissions to 
occur only in the fourth, fifth, and thermal groups. The resulting 
k e ff for these conditions is pessimistic. The values of the flux at 
thermal energy are useful in determining time to start up the reactor 
after a shutdown and will be discussed more fully in the reactor control 
section. The number of fissions is generally well divided among the 
lower three energy groups as shown in the table, which indicates an 
epithermal reactor. To obtain a thermal reactor, an increase in size 
would probably be necessary. Such an increase would be detrimental to 
airplane performance. 

The weights of the constituents in the core, the core weight, and 
the core constituent atomic densities are given in table VI. The moder- 
ator amounts to a little over 70 percent of the weight of the core, the 
latter being about 2500 pounds. 

Power and uranium distributions . - The nuclear calculations were 
made assuming a uniform fuel distribution, as discussed previously. The 
fuel-element and moderator temperatures were based on uniform radial 
power distribution, which results from a nonuniform radial fuel distri- 
bution. In the axial direction, these temperatures were based on a 
"chopped" 2/3-cosine-law power distribution, which results from an approx- 
imately uniform axial fuel distribution. Calculations were made to de- 
termine the radial fuel distribution for the uniform radial power dis- 
tribution assumed in the first temperature calculations with UOg and 
then with UC as the fuel. The reactor core was divided into 10 radial 
segments, the thickness of each being the same. The required volume of 
U0 2 (or of UC) in each segment for the startup case (140.5 lb of uranium) 
was calculated and compared with the available volume in the "meat" of 
the fuel elements in each segment. The details of the calculations are 
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not included herein except to note that the assumption was made that the 
radial flux distribution for constant radial power was the same as that 
for uniform fuel distribution. The flux for the latter was part of the 
nuclear results obtained with the computing program. 

The results of these calculations are shown in table VII. In the 
center of the reactor core (reactor sector l) there is no volume avail- 
able because no fuel elements are in this part. The two central control 
rods and through-bolts occupy the space near the core center. Compari- 
son of available and required volume for the fuel shows that except at 
reactor sector 1 there is enough volume to handle the UC. This is not 
the case for the UOg. It was for this reason that it was finally decided 

to use UC in the reactor. Because less volume would be required by the 
UC than by the UOg, a stronger fuel element would result. The conduc- 
tivity of UC is much higher than that of UOg; hence the use of UC will 
improve the dissipation of heat. Reference to table VII shows that the 
UC required in sector 1 could be put into sector 2 in addition to that 
already in sector 2. 

Comparison of the core composition using UC with that using UO^ 

(the nuclear calculations were made using the latter fuel) showed little 
difference. Consequently the nuclear calculations are considered ade- 
quate for UC. It is also thought that these nuclear calculations , from 
which & e ff vas obtained and which were made for a uniform fuel distri- 
bution, are a fair approximation with regard to uranium required for 
uniform radial power distribution. 

In addition to determining the fuel distribution for the conditions 
used for table VII, the fuel distribution using UC only was determined 
for a uniform radial power distribution and for the axial distribution 
resulting in satisfactory fuel-element stresses. Because no nuclear 
calculations were made with these power distributions, the radial and 
axial flux distributions for this case were again assumed to be the same 
as those for uniform fuel distribution. The results are shown in figure 
36. The sketch at the top shows an end view of half the reactor core 
and the available volume for UC in each of the 10 radial segments. The 
reactor was divided into 20 axial segments. Thus each outer semicircular 
segment of 1.2-inch length has a volume of 2.255 cubic inches available 
for UC. The volumes required in each segment are shown at the bottom of 
the figure. The numbers apply to the entire semicircle and for this 
reason were not repeated in the right half of the bottom sketch. The 
segments for which the volumes available are not adequate to hold the UC 
required are cross-hatched in the figure. The available volume at the 
center of the reactor is zero for the same reason as that discussed in 
connection with table VII. It appears from the figure that many seg- 
ments are inadequate, and some of the fuel in these segments must be 
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shifted to near segments. Some of the available volumes have no appar- 
ent pattern radially because of control rods occupying positions in 
which fuel elements might be placed. 

No nuclear calculations were made with this fuel distribution. The 
total amount of uranium in figure 36 is 140-5 pounds, the amount deter- 
mined as required for the startup with uniform fuel distribution. It 
would be expected that, for the power distributions used to obtain fig- 
ure 36, the effective multiplication factor would not be appreciably 
different from the values shown in figure 33. These values were deter- 
mined for uniform fuel distribution. It is expected, however, that the 
total fuel quantity required would be much greater than the 140.5-pound 
total of figure 36. Further nuclear analysis of this reactor will be 
required if future work is to be done on this project. 

Moderator characteristics . - The moderator blocks in the core were 
made of BeO, and the shape and size are given in figure 25. Temperatures 
and stresses of the blocks were calculated using the methods in appen- 
dixes E and F. The outside wall temperature Tf and inside wall tem- 
perature T e of the moderator blocks and the helium temperatures Tg e 
axially through the core are given in figure 37 for the power distribu- 
tion that resulted in the satisfactory fuel-element stresses of figure 
31. The temperatures continue to rise axially through the core, the 
highest temperature being about 2370° R. These are "safe" absolute tem- 
peratures, being far below the melting point. 

Figure 38 shows the tensile strength of BeO as a function of temper- 
ature as obtained from reference 15 (p. 843). This was used with the 
temperatures Tp of figure 37 to determine the allowable axial stress 
variation through the core shown in figure 39. The actual thermal 
stresses calculated using the temperature results of figure 37 are also 
shown in figure 39. The results show that there is a large margin of 
safety between the calculated allowable and actual moderator stresses. 

Reflector characteristics . - The side reflector was made of Be for 
strength reasons, as mentioned previously. The methods used for deter- 
mining the temperatures in the reflector pieces axe given In appendix H. 
The side reflectors can be considered shields for the pressure shell, 
because the radiation is attenuated through them. Five percent of the 
reactor power was considered to be dissipated outside the core. From 
this value, the heat flux could be obtained as a starting value from 
which the heat source in each reflector piece could be obtained as shown 
in appendix H. It was assumed that the coolant flows through each pas- 
sage, including thermal shield passages, will be metered so that the 
coolant temperature in each passage at a given axial position from the 
inlet will be equal. This will result in a fixed ratio of coolant weight 
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flows in adjacent passages throughout the length of passage, and the max- 
imum temperature in a part will remain at the same radial position for 
all axial positions. 

The general method of design was to assume the temperature of the 
reflector pieces at the coolant exit (hecause this is where the highest 
temperatures will occur) and to calculate the thicknesses of the slabs 
from a thermal stress consideration. Only thermal stresses were con- 
sidered, because these members carry small loads and can be designed 
with no restraints. After the thicknesses were determined, the temper- 
atures were calculated to see if the temperatures were satisfactory 
(i.e., below the 1000° F assumed for the stress calculations). The 
properties, thermal conductivity, gamma radiation absorption coefficient, 
modulus of elasticity, coefficient of thermal expansion, and Poisson's 
ratio, of the Be reflector pieces used in the calculations were obtained 
from references 15 and 16. 

The side reflector must be divided into a series of concentric cyl- 
inders to provide for reflector cooling and to limit the thermal stresses. 
In addition, the inner portion of the side reflector is used also as a 
framework for the core structure (fig. 6). Therefore, the inner re- 
flector must not distort in service and should be designed for low 
thermal stresses. A design stress of 1000 psi was selected for the In- 
ner reflector piece. The corresponding creep rate is much less than 
0.001 percent per hour (p. 76, ref. 16) at a temperature of 1000° F, 
which is the highest temperature expected in the reflector. Based on 
the simplified calculations of heat generation given in appendix H and 
flat-plate thermal-stress theory (ref. 15, p. 703), the allowed thick- 
ness of the inner portion of the side reflector is 0.50 inch. The re- 
mainder of the side reflector is made from cylinders having wall thick- 
nesses of 0.75, 1.00, and 1.50 inches. The largest thermal stress (5000 
psi) occurs in the 1.50-Inch piece. This stress value will be relieved 
by plastic deformation at the operating temperature. 

A sketch of a plan view of the four side-reflector pieces showing 
the channels between them is given in figure 40. The coolant-flow rates 
in each passage, the helium in and out temperatures, and the reflector 
temperatures at the coolant exit are shown. The latter are shown an 
curves on each piece with the temperature scale to the left. The high- 
est temperature shown is about 1380° R (920° F), which is below the 
1000° F used for the properties calculations in the determination of the 
stresses. 


Thermal Shields and Pressure-Shell Characteristics 

The pressure shell must be designed to contain the high-pressure 
helium and to withstand thermal stresses due to neutron and gamma 
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heating. Consequently, the pressure shell and internal thermal shield 
must he considered as a system in the design. Furthermore, since the 
internal thermal shield and pressure shell will serve as a part of the 
total biological shield, the weights of these components are not primary 
design criterions (except for handling problems) in a unit shield 
configuration. 

Since the pressure shell is the basic structural member of the core 
assembly, it must be constructed of a material with good creep proper- 
ties at operating temperatures of about 1000° F. An austenitic stain- 
less steel (type 347) will creep at a rate of approximately 1 percent in 
10,000 hours at a stress of 25,000 psi and 1000° F temperature (ref. 16). 
This material wan chosen for the pressure shell. Because of the chrom- 
ium in this steel, it will be radioactive when the reactor is shut down. 
An alloy such as Hastelloy alloy B could probably be used with less 
radioactive effect. Creep data on this metal, however, are scarce. 

Mild steel was selected for thermal shielding. 

On the basis of the simplified calculations of neutron and gamma 
heating rates described in appendix H, it was determined that 6.5 inches 
of steel thermal shielding would limit the total stress in the pressure 
shell to 25,000 psi if the pressure- shell thickness was 1.75 inches. 

(The properties of type 347 stainless steel were taken from ref. 16. ) 

The same limiting stress in the pressure shell could be obtained with 
less thermal shielding and a thinner pressure shell. However, since 
iron is also required for the biological shield exterior to the pres- 
sure shell, no detailed analysis to determine the minimum stress or min- 
imum thicknesses of the thermal shields and pressure shell was made. 

Plain carbon steel (SAE 1040) can serve as material for most of the 
thermal shielding, because corrosion is not a problem, and adequate creep 
strength can be achieved at the expected peak operating temperature of 
1000° F. Molybdenum is used as thermal shielding in the hot gas stream 
at the outlet of the reactor. The cylindrical thermal shield, which is 
6.5 inches in thickness, must be divided into six concentric cylinders 
(fig. 7) to allow for cooling and to reduce thermal stresses. The thick- 
nesses of the individual cylinders are 0.5, 0.75, 1.00, 1.25, 1.50, and 
1.50 inches. Thermal stresses are less than 8000 psi, and the resulting 
creep rates are less than 1 percent per 1000 hours for SAE 1040 steel 
operated at 1000° F. Cooling of the thermal shielding located at the 
ends of the core is provided by the holes shown in the shields in 
figure 6. 

The results of calculations of the cylindrical thermal- shield tem- 
peratures are shown in figure 41 for the final reactor design conditions . 
This figure is a sketch similar to that for the side reflectors in fig- 
ure 40. The shield temperatures shown are those at the helium outlet, 
because this is where the highest values occur. The methods of calcula- 
tion are given in appendix H. The helium flow rate w He is shown in 
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each passage, the criterions of design being the same as that for the 
side reflectors discussed before. The highest temperature, which 
occurred in the second thermal shield, was about 1390° R (930° F), and 
then there was a general reduction in temperature to an average of about 
880° P in the shield next to the pressure shell. Since the stresses 
given previously were based on properties at 1000° F and were satis- 
factory, these temperatures will provide a greater margin of safety. 


Control 

The reactor control primarily regulates the source energy required 
by the turboprop powerplants. That is, the nuclear components should 
provide performance and stability consistent with the over-all system 
demands. Other important considerations in the nuclear control design 
are reliability and safety. Normal operation results in components oper- 
ating continuously over extended periods of time. During operation the 
safety of the personnel and the vehicle demands a minimum of dangerous 
excursions or reactor scrams. 

A potential for large excursions exists in this reactor design, 
since the fuel loading results in a large excess reactivity. The excess 
reactivity to be controlled varies substantially with operating condi- 
tions. Table V, discussed in the analysis of the nuclear characteris- 
tics, shows a change of reactivity of 6 percent from a cold clean to a 
hot clean conf iguration. This variation in reactivity, added to lesser 
amounts from fission product poisoning and fuel depletion, results in 
the need for continuous shim control of the reactor. 

A block diagram showing the more important parameters that affect 
reactivity is shown in figure 42. The forward path has a transfer func- 
tion composed of the basic reactor kinetics with power or neutron flux 
as an output and reactivity as an input. Power changes, in turn, affect 
reactivity and are shown by the feedback paths that close the loop. The 
feedback paths according to their environment are classed in two groups : 

(1) the internal feedback path that is inherent in the core design and 

(2) the external feedback paths composed of the control rods and the 
reactor load. A more complete discussion of nuclear control loops is 
given in references 12 and 17. 

Poisoning effects . - Calculations were made for the individual 
worth of poisoning resulting from equilibrium and peak xenon-135, equi- 
librium and peak samarium-149, and the other fission products. These 
calculations supplement the steady-state limped values obtained in the 
nuclear characteristics section. 

The methods of calculating xenon concentration as described in ref- 
erence 12 were applied to this high-temperature epithermal reactor. 
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Equilibrium xenon poisoning, a balance between the rates of forma- 
tion and loss, can be calculated from the following equation: 


Po = 


Xe o g a,Xe 

2 a,U 


a a,Xe^ T I + T Xe 




a,U 




( 1 ) 


The development of the equilibrium poison equation can be found in ref- 
erence 12. Several of the parameters of equation (l) vary with energy, 
namely the xenon cross section, the ratio of the fuel fission to total 
cross sections, and the neutron flux. The xenon cross section varies 
greatly with energy in the region of interest, as shown in reference 18. 
Thermal energy associated with the high moderator temperature at rated 
power is 0.102 electron volt. At this energy, the xenon cross section 
is at a threshold of sharp reduction and has a Maxwellian average of 

— 

1.32x10 barns. The next energy level, group 5, covers the energy range 
from 0.411 to 50.6 electron volts, with a corresponding xenon cross- 
section variation from 8.5x10^ to 28.5 barns. The average xenon cross 
section in the fifth group is several orders of magnitude below that of 
the thermal group. This results in negligible poison contributions from 
the epithermal group. 


Equilibrium xenon poisoning was calculated for two conditions: 

(l) a freshly loaded core and (2) a partially depleted core with neutron 
flux shown in table V. The values of poisoning equivalent to reactivity 
of -0.0107 and -0.0112, respectively, were obtained. These values are 
adjusted to consider that approximately one-third of the fissions occur 
at thermal energy. 


Xenon peak poisoning results when the xenon formation rate exceeds 
its loss rate because of rapid decreases in the neutron flux. Neutron 
capture by xenon becomes ineffective, causing xenon to build up through 
a maximum as a function of time. The time variation of xenon from equi- 
librium conditions for a step decrease in neutron flux as derived in 
reference 12 is as follows : 


p(t) 


ff / "V 

_( A Xe " 2 a,U V 



( 2 ) 


The xenon buildup affects the thermal region only, as in the equi- 
librium xenon calculations. Since the thermal flux is relatively low, 
with approximately one-third of the fissions in this region, the peak 
xenon is limited to small values . Under the most stringent conditions 
(i.e., a step decrease in neutron flux from rated power), the negative 
reactivity remains below 2 percent. Figure 43 shows the xenon time 
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variations for these conditions with a freshly loaded core and a par- 
tially depleted core. The maximum values are reached in about hours. 

The over- all xenon poisoning is greatly reduced in this high- 
temperature epithermal reactor over what would be expected from a low- 
temperature thermal reactor. A slight additional fuel inventory in the 
core provides the ability to override xenon, resulting in a safer and 
more flexible system. Since xenon "buildup" is limited to small values, 
the related problem of "burnout" upon restart at peak values is also of 
a minor nature. 

The effect of samarium poisoning is calculated in a manner similar 
to xenon poisoning. The equilibrium poison equation is simplified, 
since there is no direct samarium yield from fission and no decay of 
samarium. Therefore, the equilibrium equation reduces to the following 
equation, which is independent of neutron flux and samarium cross 
section: 

Zf 

P ° ' 

The equilibrium poisoning results in a negative reactivity of 0.0114 
and reaches equilibrium in approximately 35 days . Since the time to 
reach equilibrium is long, the effect on criticality is fairly negligible. 

Samarium poisoning as a function of time for a step decrease in flux 
is very small and can be neglected. 

Temperature effects . - The temperature variations within the reactor 
have an important effect on reactivity. These temperature changes are 
caused by local variations in the fission energy and variations in the 
coolant parameters. The helium coolant is a non-neutron absorber and 
has a relatively small moderating value. Therefore, the coolant does 
not directly influence reactivity. Changes in coolant-flow parameters 
and in the heat-exchanger load do affect reactivity by varying reactor 
temperature, but in a more gradual manner. 

A change in reactor temperature alters the multiplication factor in 
two ways. The change in the mean energy of the thermal neutrons affects 
the thermal cross sections, and density changes affect buckling and leak- 
age. The first effect, which is called the nuclear temperature coeffi- 
cient (ref. 12), is the stronger and quicker and therefore the more im- 
portant. Table V shows a 6-percent drop in k e ff from the cold clean 

to the hot clean configuration. This amounts to an over-all temperature 
coefficient of -3.76x10“^ 5k/°F. Although there are other factors af- 
fecting the total -temperature coefficient, this value should be repre- 
sentative and contribute toward the inherent stability of the reactor. 
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Fuel bumup . - The calculations presented in the nuclear character- 
istics section based on an estimated 1000-hour load life resulted in ef- 
fective multiplication factors as shown in figure 33. The following 
analysis shows the rate of fuel bumup and the expected operating time 
before reloading. 


The following equation gives the fuel bumup rate per megawatt of 
power: 


T t 235 Bumup rate _ .. ... p a,U grams 
Unit energy ~ * Of Mw-day 


= 2.29X10 -3 


p a,U lb 
Of Mw-day 


(4) 


The ratio CT a ,u/ p f Yar i eB somewhat with energy, and the average value 

at thermal energy 1.228 was used in the calculations. The rate of 
bumup per megawatt adjusted for nonfission capture is 2. 81*10“ ^ lb/Mw- 
day. At this rate it takes 37.4 days (approximately 900 hr) of contin- 
uous operation at rated power to burn up 10.5 pounds of fuel. In this 
time the reduction of k e ff from fuel bumup is 1.0 percent, as shown 
in figure 44. Loss of reactivity from poisoning reduces k e ff another 
3 percent to a value of 1.08. Considering that peak xenon contributes 
less than 1-percent negative reactivity at this point, there is a margin 
of reactivity of better than 7 percent. A rough extrapolation indicates 
that reloading of the reactor will be required after 195 days of opera- 
tion at full power. This extrapolation does not take into account the 
fission spectrum shift towards a more thermal reactor and the associated 
increase in xenon and samarium poisoning. The actual flight plan does 
not call for continuous operation and will affect the bumup rate also. 
Therefore, the reloading time is an approximate number to be adjusted 
when more exact calculations are made. 


Control rods . - The maximum excess reactivity that occurs with the 
cold clean configuration is 18 percent. This excess reactivity consid- 
ered with a shutdown reactivity of 10 percent results in a total worth 
of 28 percent. Therefore, the design of the control rods has a worth of 
28 percent distributed between the shim and safety systems. 

The tentative control-rod system is composed of two scram rods, 
three combination shim-scram rods, three pairs of shim rods, and two 
regulating rods. Figure 7, which is the end view of the reactor, shows 
the radial position of the various rod types. The actuators used to 
drive the control rods are powered by electric motors, hydraulic servos, 
or pneumatic servos. Table VIII lists the drive mechanisms and worth of 
each rod when all the rods are in the core. 
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Two regulating rods with a worth of 0.005 each are located at the 
center of the reactor. One of the regulating rods will normally be in 
a standby position fully withdrawn from the core. In the event of a 
malfunction, the standby rod is used with an effective worth of 0.005 
or less, depending on the shadowing of the primary rod. The regulating 
rods have individual drive mechanisms located within external thimbles, 
as shown in figure 5. Linear reluctance motors similar to those de- 
scribed in reference 17 provide the drive motion directly without the 
need of gearing. 

Shim control during operation at power level is attained by three 
pairs of rods located 120° apart and half the radial distance from the 
center to the reflector, as shown in figure 7. Each pair of rods has a 
worth of 0.04 and is motivated by a common synchronous reluctance motor. 
The rotary motion of the motor is transformed by gearing, screw, and a 
movable nut to provide linear motion of opposite directions to the indi- 
vidual rods. The boron loading of each rod varies in a linear manner 
over its length but in opposite directions, as shown in figure 5. This 
results in a worth of the pair of rods that is equal at all times along 
the length of the reactor core, therefore minimizing axial power dis- 
tortion normally associated with a rod partially inserted in the core. 
Variation in total worth is accomplished as the rod pairs are simultane- 
ously inserted and extracted. 

Dual-purpose rods used for both shim control and scram are located 
equidistant between the shim- rod pairs. These three rods have individ- 
ual worths of 0.04 each and are used as shim rods for starting. Upon 
reaching power level, the rods are almost completely extracted from the 
core and can be used as scram rods exclusively. The drive mechanism can 
be hydraulic or pneumatic with slow speed operation in or out in conjunc- 
tion with a special dump port to provide fast insertion for scram. The 
pneumatic system seems more advantageous because it can make use of the 
same inert gas, helium, used as the coolant. The high temperatures, 
radioactive field, and possibility of contaminating the coolant by leaks 
make the hydraulic system less favorable. Individual rods have their 
own drive that can be externally connected to work in unison or singly. 

Two scram rods of relatively small worth, 0.02 each, are located 
180° apart at a radius larger than the shim and shim-scram rods, as shown 
in figure 7. These rods are pneumatically operated with slow extraction 
and very fast insertion speeds. The position of these rods from start 
to shutdown will be fully extracted. 

All the control rods are cooled by helium bled off the main supply 
and passed through an external air heat exchanger. This cooled helium 
at 950° R enters the control-rod sleeves near the actuators, passes 
through the sleeve around the control rod, and mixes with the primary 
coolant to flow through the reactor core. 
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Automatic control in the power range is obtained by varying the 
regulating rod, with followup shim- rod control when the regulating rod 
reaches a limit. 


Reactor Structure 

Core structure . - The geometry of the 640 fuel elements and the 
moderator blocks has been described, and figure 45 shows a section (E-E 
of fig. 5) of them through the core. Figures 5 and 7 show plan and end 
views of their assembly in the reactor, and figure 25 shows cross sec- 
tion D-D of figure 5. Figure 46 shows a complete assembly of a fuel 
element and moderator block. The fuel element is continuous, but the 
moderator block is divided into six lengths, as shown in figure 6. The 
lengths between the outer moderator block and fuel-element support 
plates 0SP1 and 0SP2 (see fig. 6) and the inner moderator block and 
fuel-element support plates ISP1 and ISP5 form the end reflectors. Five 
of these inner support plates are used in the core, as shown in figure 
6. The detail of these plates is shown in figure 47. The outer support 
plates are similar, except that larger holes (l-in. diam. ) are in them 
to give more area for the flow of greater volume of hot helium at the 
reactor exit. 

The principal problem encountered in the pressure-vessel-enclosed 
reactor design was the differential radial expansion between the core, 
reflector, thermal shield, and pressure vessel. This incompatible dif- 
ferential expansion results from the varying expansion rates of the dis- 
similar materials used in these structures and the high negative radial 
temperature gradient from the core to the pressure vessel, especially 
during startup. 

To eliminate this radial expansion problem and the stresses result- 
ing from expansion interferences, the supporting structure of the core, 
reflector, and shields was designed to allow these parts to expand freely 
radially. This was accomplished without leaving these parts free to 
cause impact loading between each other during periods of acceleration 
or deceleration. It was especially important to protect the brittle BeO 
moderator blocks of the core from any type of impact loading. 

In the core the BeO moderator blocks are individually supported by 
the molybdenum plates shown in figure 47 and similar ones at the outlet, 
so that clearances for cooling are provided between the blocks. The Mo 
plates are each supported by 24 radial pins around their circumference 
as shown in the figure. These radial pins fit into radial holes in the 
plates and the first Be side reflector (AR1 of fig. 6), which is also 
used as a support housing for the core. These radial, pins allow free 
radial expansion between the core structure and the side reflector, and 
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also the horizontal and near-horizontal pins allow the first reflector 
to help support the core weight. The pins finally fit into slots in 
the second reflector (see fig. 6, reflector AR2), which also helps to 
support the core weight. Similarly, in the core the moderator blocks 
and support plates are held together axially by low-stressed through- 
bolts, making it possible to use thin plates (l/8 in. thick) to support 
the weight of the moderator blocks and fuel elements. The moderator 
blocks held together in this manner act as stiffeners for the support 
plates . 

The control rods are supported independently from the core with 
necessary clearances to isolate them from the core. Further details of 
this structure will be given later. 

The fuel elements are positioned in the core in the following man- 
ner: Each element is attached to an end positioning block by a Mo 

spring spacer as shown in figure 48. The spacer is brazed to both the 
positioning block and to the fuel element. The positioning blocks fit 
into holes in support plate 0SP2 as shown in figure 6. The remaining 
spring spacers that support the fuel element (fig. 25) are brazed to 
the Mo support bushings shown in figure 46. The five inner support 
bushings are positioned axially by the moderator blocks, while the end 
support bushing at support plate 0SP1 is positioned axially by a pin 
through the moderator block. The fuel-tube - positioning-block assembly 
is slid through the spring spacers and is positioned axially by the 
block, which has a protuberance as shown in figure 48. The protuberance 
bears on the shield plate RSI. The positioning block is held between 
the moderator block and the plate RSI. 

The core assembly is surrounded by a thin molybdenum cylindrical 
flow divider as shown in figure 6. This divider separates the helium 
flowing through the side reflectors and thermal shields in one direction 
from the helium flowing through the core in the other direction. 

Reflector structure . - The structural features of the first side- 
reflector cylinder AR1 have been discussed in the preceding section. 

The weight of the core i6 further supported through the Be reflector 
cylinders AR2, AR3, and AR4 by means of spacers between them, as shown 
in figure 6. The four reflector cylinders are centered and supported 
by four radial vanes, in both the front and rear of the reactor, over 
which the cylinders slide (see fig. 6). The vanes at the rear of the 
reactor are shown in figure 49. The four front radial vanes are 
attached to the front support shield FS5, and the four rear vanes to the 
rear support shield RS4 (see figs. 5 and 49). FS5 is centered and sup- 
ported by the centering pin in the front, and RS4 by the centering col- 
lar in the rear. 
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Thermal shields and pressure-shell structure . - The internal shield- 
ing consists of the six side cylindrical shields AS1, 2, 3, 4, 5 , and 6 
for which the stresses and temperatures were previously given, four cir- 
cular front shields FS1, 2, 3, and 4 (fig. 6), front shield FS5 (fig. 6) 
which has been described, three circular rear shields RSI, 2, and 3, (fig. 
6), five rear annular shields RAS1, 2, 3, 4, and 5 (fig. 6), and rear 
shield RS4 (figs. 6 and 49) which has also been described. The front 
shield FS1 and rear annular shield RAS1 are shown in figures 50 and 51, 
respectively. The other shields are similar to the two types shown in 
these figures except that they vary in thickness. Included between RS3 
and RS4 are a plenum and ducts for the helium leaving the reactor. The 
thin cylindrical flow separator around the core butts against this ple- 
num chamber. A cylindrical annular support AS (fig. 6) is placed be- 
tween the rear circular and annular shields. 

The radial vanes described before (the rear ones are shown in fig. 

49) also center and support the iron shields. These supports are such 
that they will give positive location at assembly and allow relative 
expansions during startup and shutdown. The shield structure is axially 
positioned by four clamping bars that are axially fixed in the pressure- 
vessel head. Some details of the shields, ducts, vanes, and clamping 
bars are shown in figures 7 and 52, which are sections C-C and B-B, re- 
spectively, of figure 5. 

It is expected that the shield structure and pressure vessel will 
creep because of the weight and pressure loads at the design operating 
temperatures . For this reason spacers have been placed on the annular 
shields to maintain minimum cooling -pass age openings and control-rod 
housing clearances. 

Axial motion is not considered serious, since the axial accelera- 
tions are not expected to be as high as the vertical and turning accel- 
erations. Provision to prevent axial impact loads due to creep loosen- 
ing have therefore not been made; they could be incorporated, however, 
if necessary. 

Further details of the way the foregoing structure is supported and 
assembled in the pressure shell can best be obtained from appendix I, 
which gives a complete listing of the reactor components and a step-by- 
step assembly procedure. Following the procedure will clarify the struc- 
tural picture of figure 6. 

Control-rod structure . - The control rods are loaded with (Mo^Bg 
and have a cladding of molybdenum around this material, as shown in fig- 
ure 53. Each three-layer rod is attached to a molybdenum cylinder by 
means of molybdenum spring spacers as shown in figure 53. The latter 
cylinder then slides in a stationary molybdenum sleeve shown in figures 
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5 and 45. As mentioned, low -temperature helium flows through a hole in 
each rod as well as around the rod in the annulus formed by the rod 
outer cladding and the cylinder to which it is attached. Higher- 
temperature helium also flows through an annulus formed by the support 
sleeve and the moderator blocks- Details of the helium flow can be ob- 
tained from figure 5. 

The control rods are independently supported from the core with 
necessary clearances to isolate them from the core. The sleeves are 
supported in front by knife-edge supports in FS5 (see fig. 5) and in the 
rear by similar supports on split positioning rings, shown in figure 5, 
in the pressure-shell head. The positioning rings are held in place by 
housings in which the control-rod mechanisms are placed. This type of 
support for the control rods ensures that the control -rod housings will 
not be loaded by the core expansions and movements and cause the control 
rods to jam. 


Reactor Assembly and Stresses 

The reactor structural features have been discussed in some detail 
in a previous section of the report. The method of assembly of the re- 
actor is given in appendix I. This appendix augments the discussion 
given previously so that a clear picture can be obtained of the reactor 
structural aspects from it and the material presented herein. A detail 
worth mentioning is that the fuel-element tubes are designed for removal 
from the rear end of the core. Replacement is possible by removing the 
pressure- vessel head and the rear reactor shield assembly. After these 
are removed, the fuel tubes can be slid in or out for replacement or 
checking. The fuel tubes are held in place in the reactor by position- 
ing blocks attached to the fuel tubes and clamped between 0SP2 and RSI 
(see fig. 5). This positioning arrangement permits the tubes to expand 
freely and also facilitates the replacement of the tubes. 

Calculations of some reactor stresses were made, and the results 
are shown in table IX. They do not include internal stresses due to 
high internal thermal gradients. The major stresses tabulated for the 
pressure-vessel wall and flanges are due to the internal helium pressure. 
Gussets were added to flanges to reduce bending stresses and to reduce 
maximum combined stresses to 25,000 psi or less for low creep rates. 

The primary stresses tabulated in table IX for the reactor core, 
control rod, and shield components are for weight loads, the assembly 
of these components being such that no thermal-interference stresses 
occur between mating parts. The weight loads were based on the opera- 
tion of a logistic airplane with low maneuvering and landing accelera- 
tions during normal operations. The assumed normal maximum acceleration 
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in the vertical direction was 3 g's. Side and longitudinal accelera- 
tions are expected to be less. The normal maximum angular acceleration 
about the longitudinal axis was assumed to be 6 radians per second^. 

Six radians per second^ is equivalent to a vertical acceleration of 
3 g's about the longitudinal axis of a point on the wing 16 feet from 
the fuselage centerline. 


The maximum deflection for the control-rod housing tube is given 
for a 3-g acceleration load for a shim-scram rod in the most unfavorable 
position. With this low deflection, these rods should be able to move 
freely during normal maximum expected accelerations. 


BIOLOGICAL SHIELD 
Primary Biological Shield 

The determination of the size and weight of a unit biological 
shield was accomplished by use of the procedure discussed in appendix J. 

7 For the typical airplane considered herein, flying at 0.72 Mach 

^ number at an altitude of 30,000 feet, calculations revealed a required 

power of 98.5 megawatts. The dimensions of the reactor, cylindrical 
in shape, were found to be 24 inches in length and 34.66 inches in diam- 
eter. A 4- inch BeO reflector surrounded the reactor. Thermal shields, 
with provision for cooling, and a pressure shell of maximum 4-inch thick- 
ness surrounded the reactor and reflector. The thermal shields and 
pressure shell were considered part of the gamma shield. Additional 
required gamma and neutron shielding, exterior to the pressure shell, 
was determined so as to yield a crew-compartment dose rate of 0.025 rem 
per hour at a distance of 90 feet from the reactor-core center. 

o 1° 

The configuration was divided into six 15° and four 2 2^ sectors , 

as shown in figure 54. The shield external to the pressure shell was 
shaped (the thicknesses of gamma and neutron shields varied for each 
sector) so that for the specified total dose rate of 0.025 rem per hour 
at the crew compartment (for air scattering and direct radiation) , the 
shield weight is minimum. Initial calculations were made for 100 mega- 
watts and a 100-foot separation distance; the results could then he ad- 
justed for the 98.5-megawatt and 90-foot case. 

The assumptions and approximations made in the shield evaluation 
are discussed in appendix J. Any attempt at an analytical evaluation of 
a shaped shield requires many simplifying assumptions, and at best the 
shield weights obtained are approximate. The shield method used is 
based on core radiations only, and the very important source of gammas 
from captures in the shield is not taken into account. It is hoped 
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that borated materials used in the shield will reduce these capture 
gammas to the extent that neglecting them does not cause considerable 
error in the shield-weight results . 

The first set of calculations was made using various combinations 
of gamma and neutron shielding materials. The results of these calcu- 
lations, for 100 megawatts and 100 feet, presented in table X in increas- 
ing order of total shield weight, reveal considerable differences in 
weight for the different combinations of materials. The combination 
Fe-Fe-decalin-LiH was selected for the present reactor shield. Reasons 
for this choice follow. The typical airplane is of a conservative de- 
sign, and one which could be built almost immediately. Iron (or borated 
steel) was selected for use in the inner gamma shield, pressure shell, 
and outer gamma shield because this material can withstand the temper- 
atures and pressures imposed, can be fabricated in the sizes necessarily 
required for such an application, is a reasonably good gamma shield, and 
is relatively inexpensive. Some other of the gamma materials considered 
would result in a lighter shield, but fabrication with these materials 
to the required sizes might require a long development program and their 
costs might be prohibitive. A second requirement imposed in the current 
design was the use of a chemical fuel as part of the neutron shield. 

This fuel could be burned in cane of emergency. Decalin was chosen for 
this fuel, and lithium hydride for the remaining neutron shield. The 
total shield weight including the thermal shields and pressure shell for 
the chosen materials for the 100-megawatt 100-foot case was 103,000 
pounds (see table X) . Of this total, 23,000 pounds were decalin; this 
weight amounts to about 3100 gallons. Table XI presents the calculated 
exterior gamma- and neutron- shield thicknesses for each sector used in 
this weight determination. If depleted uranium were used for the gamma 
shield throughout and lithium hydride for the neutron shield, the weight 
could be reduced to 73,000 pounds (see table X) . 

After selection of the shielding materials (Fe-Fe-decalin-LiH) , cal- 
culations were made for different reactor powers and different separa- 
tion distances. The thermal-shield and pressure- shell thicknesses, 
reactor-core length, crew dose rate, and reflector thickness were all 
held fixed, but the core diameter was varied by the relation 

Core diameter = 3. 466-^Reactor power in Mw in. 

The diameter is proportional to the square root of the power, since only 
the number of tubes and the flow through them varied. Reactor powers of 
50, 100, 150, and 200 megawatts and separation distances of 50, 100, and 
150 feet were considered. The results of the calculations are presented 
in table XII. 

The variation of shield weight with altitude was calculated for the 
100-megawatt, 100-foot, 0.025-rem-per-hour case. The scattered neutron 
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and gamma dose rates were altered by consideration of the variations in 
air density. The results are shown in figure 55. 

The variation in shield weight with dose rate for the so-called 
standard condition (100 Mw, 100 ft, 50,000 ft) and for the chosen shield 
materials is given in figure 56. 

Calculations were also made for the standard case with dimensions 
exterior to the reflector held constant, but with different reflector 
thicknesses and the corresponding variations in core size. For example, 
the following cases were considered: 


Reflector 

Core 

Core 

thickness. 

length. 

diameter, 

in. 

in. 

in. 

2 

28 

38.66 

4 

24 

34.66 

6 

20 

30.66 


The total shield weights for these three cases differed by only 700 
pounds . 

Other calculations were made for the standard case with the reactor- 
core size fixed, but with reflector thicknesses of 2, 4, and 6 inches. 

The thicknesses of the thermal shields and pressure shell were held 
fixed, but the distance from the reactor-core center to the pressure 
shell varied as the reflector thickness was changed. The total shield 
weights obtained were as follows : 


Reflector 

Shield 

thickness. 

weight. 

in. 

lb 

2 

95,000 

4 

103,000 

6 

113,000 


The final set of calculations was made, again for the standard con- 
ditions, but with the decalin replaced by either JP fuel, water, or void. 
For JP fuel or for H 2 0, the dose rate remained at about 0.025 rem per 
hour. However, with the void, the dose rate increased to 295 rems per 
hour. Hence, once the reactor is off and the decalin is used for fuel, 
the reactor must not be restarted until the decalin tanks are refilled. 
However, if the plane happens to be in a location where decalin is not 
available, water can be used as a substitute. This, however, will remove 
the safety factor originally leading to the use of decalin and will also 
increase the shield weight, but it may make possible the return of the 
plane to a site where decalin can again be used. 
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Additional Shielding Required for Streaming 
Through Ducts and for Voids 

In addition to the primary biological shield, additional shielding 
is required to compensate for the presence of ducts and streaming through 
the ducts. As stated in appendix K, if considerations for these effects 
are made for neutrons, the resulting additional shielding will also be 
satisfactory for gamma shielding. 

The geometrical considerations of the present problem prevented the 
use of optimum-sized ducts, and streaming calculations were made by the 
procedure described in appendix K. It was found that an additional 
shield weight of 3500 pounds was required to compensate for streaming. 

An additional 3000 pounds of shielding was also required to compensate 
for the void in the primary shield caused by the presence of the ducts . 
This combined additional shielding was added in the form of hemispherical 
tanks of decalinj about 900 additional gallons of decalin became avail- 
able for emergency use. 


Additional Weight Required to Provide for Shield Cooling 

As yet, no discussion has been made regarding provision for cooling 
the biological shield exterior to the pressure shell. To allow for such 
cooling, the calculated diameter of the biological shield was increased 
8 inches. The additional weight required to compensate for this arbi- 
trarily selected void was calculated to be about 5000 pounds. 


Adjustment of Primary Shield Weight for the Typical Airplane 

It was previously stated that the so-called standard case (100 Mw, 
100 ft, etc.) could be adjusted to the 98.5 megawatts and 90 feet neces- 
sary for the typical airplane. This was done, and the resulting weight 
of the primary biological shield increased from 103,000 to 104,000 pounds. 


Total Shield and Reactor Weight 

The total weight of the shield plus reactor for the typical airplane 
can now be found as follows: 

Primary biological shield, lb 104,000 

Weight for ducts and voids, lb 6,500 

Weight for shield cooling voids, lb 5,000 

Reactor and reflector weight (given previously in report) , lb . . 4, 500 

Total, lb 120,000 
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Structural Features 

The previous discussion of the shield has indicated the amount of 
gamma and neutron shielding required. Part of the gamma shield (thermal 
shields) is inside the pressure shell, which surrounds the reactor. A 
suggested method of constructing the shields exterior to the pressure 
shell is indicated in figure 9. This figure shows the reactor and shield 
mounted in the fuselage. Plan, elevation, and end views are presented. 
The plan view shows that the shield and reactor will occupy a large part 
of the fuselage width. 

Construction . - The gamma shield on the exterior of the pressure 
shell is a part of the permanent airframe and is a circular yoke, as 
shown in the end view. This yoke is made of iron containing a small 
percentage of boron. A method of constructing the neutron shield, con- 
sisting of lithium hydride and decalin, is indicated in figure 9. Deca- 
lin was selected as the neutron shield forward of the reactor, and lith- 
ium hydride in the rear. The additional shielding required for streaming 
and void effects was chosen as decalin and located in tanks exterior to 
the primary shield and in direct line with the various duct legs. Fig- 
ures 6 and 9 show the proposed ducting. The lithium hydride was placed 
in containers that are stacked, like bricks, in such a way that the 
joints of successive rows are staggered. The decalin was installed in 
a series of tanks, as shown. The lithium hydride containers and the 
decalin storage tanks around and forward of the iron-boron shield are 
installed permanently. The lithium hydride and the two decalin storage 
tanks at the rear of the reactor surrounding the flange, ducting, and 
control-rod actuators must be removable for reactor installation or 
replacement. 

Mounting in fuselage . - The circular-yoke gamma shield, on which 
the reactor is designed to be mounted, is a permanent part of the air- 
frame. It is attached to the airframe structure by means of the girders 
indicated in figure 9. The reactor is first attached to a lift outside 
the airplane. The rear flange of the reactor pressure vessel has a 
circumferential groove into which the hydraulic lift fixture is secured 
by hydraulically actuated pins in the lift fixture. The lift mechanism, 
supported from a railroad track alined with the airplane fuselage, lifts 
the reactor to its vertical position in the fuselage through bay doors 
that swing open from the fuselage bottom. The front of the reactor is 
then slid into the iron-boron shield yoke, with sufficient clearance 
allowed for cooling. The reactor is then positioned in front by a 
centering pin and around the flange by a circular I-beam. This I-beam 
is a part of the permanent airframe structure. An electrical drive 
mechanism then opens and closes the split ring that positions the reac- 
tor axially by fitting in a circumferential groove in the reactor center- 
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ing pin. The ducts are welded into place after the reactor is secured 
in the fuselage structure. The lithium hydride containers and the two 
decalin storage tanks at the rear of the reactor are then installed. 


AIRPLANE OPERATION 
Propulsion-System Control 

The propulsion system is composed of eight turboprop engines with 
their individual heat exchangers and helium pumps, powered by a common 
reactor. Since the control scheme is the same for all engines, this 
section considers the operation of an individual engine. 

Basically, the control system can be divided into three Individual 
component control loops. They are the reactor nuclear control loop, the 
reactor coolant-flow system, and the turboprop engine control, as shown 
in the block diagram of figure 57. The nuclear and coolant-flow loops 
are coupled at the reactor, while the coolant-flow loop and engine are 
coupled at the heat exchanger and at the engine-driven coolant pump. 
Interaction between the loops is possible, even to the extent of an en- 
gine disturbance being reflected in a nuclear loop response. 

The philosophy of the integrated control loop design is to minimize 
interaction between components without penalizing the component controls 
of the system; that is, to arrive at a complete system that has both 
stability and adequate dynamic response time. 

The primary change from a conventional turboprop system is the re- 
placement of the conventional fuel-metering valve and burner with a re- 
actor heat source, coolant gas, and heat exchanger. A chemical fuel- 
burning engine normally uses engine parameters such as temperature and 
engine speed to vary propeller pitch and fuel flow for control. An anal- 
ogous system could be used for the nuclear-powered engine. This system 
could control the heat source by nucleonic means, or use a controllable 
coolant bypass about the heat exchanger. Some disadvantages of these 
modes of operation are as follows: (l) Changes in engine power demand 

place the burden of control on the nuclear loop, increasing the possi- 
bility of discontinuous operation; and (2) the transient response to an 
engine power demand is dependent on the nuclear period, thermal lags in 
the reactor and heat exchanger, and transport time of the coolant lines. 

A more desirable control scheme is shown in figure 58. The engine 
control consists of two loops, propeller blade pitch controlled by en- 
gine speed and power modulation from a turbine bleed control. Normal 
operation after startup is at a fixed engine speed and fixed reactor 
power. Fixed engine speed results in a constant airflow through the 
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compressor and heat exchanger. Since the coolant pump is driven hy the 
engine, it also will operate at constant speed. Therefore, the load on 
the reactor is essentially constant, reducing the number of nuclear 
transients. Engine power variations are obtained by regulating the 
amount of heated airflow through the turbine and the turbine bypass. 

With the bypass valve fully closed, all the air flows through the tur- 
bine, producing rated engine thrust. As more air bypasses the turbine, 
the rotor speed tends to drop, but the automatic speed control varies 
the propeller pitch to maintain the desired speed at a loss of propeller 
power. Therefore, power settings are changed with a minimum of inter- 
action between engine and reactor. Further detailed investigations are 
needed to study transient effects, loop stability, and automatic turbine 
bypass operation. 

An additional fuel input is located between the heat exchanger and 
turbine bypass control where decalin can be metered with an afterburner 
type flameholder. This chemical fuel system can be used as the primary 
system when the reactor is inoperative or in conjunction with the re- 
actor to increase turbine temperature up to its limiting value. The con- 
trol of the chemical fuel, whether automatic or manual, depends on 
turbine-inlet temperature . 

The heat exchanger operates near its critical temperature limit, 
and in certain phases of the flight plan it may overheat. Therefore, a 
temperature-l imi ting control is necessary to restrict reactor power to 
the temperature l imi t of the heat exchangers. 

Since the integrated control system is a complex network with a 
reactor and eight turboprop engines, it is reasonable to assume that a 
crew including both a nuclear engineer and a powerplant engineer is 
needed to assist the pilots in operation. 


Startup and Shutdown Procedure 

The startup procedure is a matter of programming the various oper- 
ations in the reactor, the coolant loop, and the engine to provide a 
safe, efficient operation. Until the power range of the reactor is 
reached, the coolant temperature is unchanged, essentially uncoupling 
the reactor from the engines. Therefore, initially the reactor flux can 
be increased to the power level and the engines started on chemical fuel 
independently. For reasons of safety, one or two of the engines should 
be operative during reactor startup. This provides circulation of the 
coolant to absorb an accidental overshoot into the power range. 

Increasing reactor power through the power range to rated power is 
a more critical operation. Several reactor limitations restrict the 
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rate of increasing power. They are the reactor period, reactor temper- 
ature, and thermal shock to the reactor components. 

The reactor period is governed by control-rod manipulations and 
is limited by the nucleonics of the system. During the positive reactor 
period, the reactor temperature will increase and must be maintained 
within limiting values . The coolant flow must increase with reactor 
power to maintain temperature below the temperature limit of the reactor. 
Since the coolant pumps are engine -driven, the engine speed loop control 
regulates the coolant flow through the reactor. Therefore, coolant flow 
or engine speed must maintain reactor temperature below its maximum 
value. Coordination is necessary between the nuclear variations and 
coolant-flow control to minimize sharp temperature changes. These quick, 
large temperature variations may cause reactor components to exceed 
their thermal stress limits. 

As the usable power of the reactor becomes available at the heat 
exchangers, the chemical power can be decreased accordingly. When the 
reactor is delivering rated power, the engines are at rated speed, with 
the excess power bypassing the turbine. To take off, the turbine bypass 
is closed, increasing the propulsion power. Chemical power is used as 
needed to augment the nuclear power. 

The shutdown process is comparatively simple. Insertion of the 
reactor control rods shuts down the reactor with a stable period of 
approximately 80 seconds. Chemical fuel is used as needed to maintain 
an adequate coolant pumping power to remove reactor afterheat. 


Off -Design Performance 

The design-point performance represents an important phase of air- 
plane operation; however, the characteristics at off-design conditions 
must be satisfactory in a worthwhile airplane system. Furthermore, some 
of the demands on the control system can be uncovered by examining the 
off-design performance. Another important point to be considered is the 
emergency range supplied by the fuel used as shielding should the reactor 
become inoperative. 

Information on the performance of the typical airplane at off- 
design flight conditions was determined by calculations based on the 
following assumptions. The calculations assumed constant shaft speed 
(engine and pump) type of operation, as mentioned in the discussion of 
the Propulsion-System Control. The temperature of the helium leaving 
the reactor was not allowed to exceed the design value of 2250° R, and 
the temperature of the air leaving the heat exchanger was limited to 
1800° R. These temperature restrictions were imposed largely because 
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of properties of the heat-exchanger material. The compressor perform- 
ance was estimated from the data of reference 4, and the pump perform- 
ance map is given in appendix D. 

Performance at various altitudes . - Calculated variations of thrust 
and drag with Mach number are shown in figure 59 for altitudes of 25,000, 
30,000, and 35,000 feet. The thrust curves were calculated for purely 
nuclear operation. The curves show that the maximum flight Mach numbers 
at 25,000 and 35,000 feet are less than at the design altitude. This re- 
sult is due to a combination of factors. At high Mach numbers, the drag 
of the typical airplane at 35,000 feet is less than the drag at 30,000 
feet because of the large wing area, as discussed previously. On this 
basis alone, an increase in Mach number with increasing altitude would 
be expected. However, the reduction in engine airflow rate with increas- 
ing altitude requires that the. helium temperature leaving the reactor be 
reduced below the design value to prevent the air temperature at the 
heat-exchanger outlet from exceeding the 13-mi ting value of 1800° R. As 
a result, the heat-exchanger material limitations do not allow the power 
potentialities of the reactor to be utilized at altitudes above design. 

A similar effect occurs at the design altitude (30,000 ft) for flight 
Mach numbers below the design value of 0.72. In this case, the engine 
airflow rate decreases with flight Mach number (changes in ram density 
ratio), and a decrease in reactor outlet temperature is again required. 
For example, at a flight Mach number of 0.46 at 30,000 feet, the reactor- 
outlet temperature must be 2170° R. 

At lower than design altitude, the maximum Mach number decreases 
largely because of the reduction in lift-drag ratio (increase in drag). 
There is, however, an additional effect caused by temperature limita- 
tions. At 25,000 feet and with the maximum reactor-outlet temperature 
of 2250° R, the turbine- inlet temperature is lower than the design value 
of 1800° R because of an increased engine airflow rate at the lower 
altitude . 

Takeoff and climb . - Chemical augmentation is required for takeoff 
and some phases of climb. Calculations showed that, with no burning, 
the power produced by the turbine was almost entirely absorbed by the 
compressor and helium pump at sea-level static (takeoff) conditions. 

At these conditions, the temperature of the air leaving the heat ex- 
changer is only 1460° R because of the high engine airflow rate. If 
enough fuel is burned to raise the turbine-inlet temperature to 1800° R, 
an adequate takeoff thrust (86,000 lb) is achieved. 

An estimate of the amount of fuel consumed during transition from 
a climb-out condition to the design flight condition was calculated. 

The climb was assumed to start from sea level at a Mach number of 0.22 
and to progress at constant dynamic head up to 30,000 feet and a Mach 
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number of 0.41. Subsequently, the airplane is accelerated to a Mach 
number of 0.72. During this maneuver, the turbine-inlet temperature is 
maintained at 1800° R by burning fuel at altitudes below 30,000 feet. 

At 30,000 feet, no fuel is used; in fact, the reactor-outlet temperature 
is, necessarily, reduced below the design value to meet the temperature 
limitations discussed previously. 

The entire climb maneuver was divided into seven steps for calcu- 
lation purposes, and the midpoint conditions of each step were used to 
compute rate of climb and fuel consumption. According to the approxi- 
mate calculations, 6 minutes are required for the climb and 200 pounds 
of fuel are consumed. These results show that the 4000 pounds of chem- 
ical fuel allowed for takeoff and climb are sufficient. 

Emergency range . - Calculations were made to determine the range 
provided by the fuel used as part of the neutron shield should the re- 
actor become inoperative. The weight of fuel is about 30,500 pounds in 
the typical airplane. Time did not allow a complete survey of flight 
and operating conditions for the purpose of finding maximum emergency 
range, and only two flight conditions (a Mach number of 0.5 at altitudes 
of 25,000 and 30,000 ft) were studied. The Mach number of 0.5 was chosen 
because, at this condition, the airplane operates close to maximum lift- 
drag ratio or peak aerodynamic efficiency. A turbine-inlet temperature 
of 1546° R was selected for the emergency range calculations. It is at 
this inlet temperature that net thrust equals airplane drag. 

The fuel used as shielding provides 1.73 hours of flying time and 
a range of 590 statute miles at 30,000 feet. At 25,000 feet, the time 
is 1.55 hours and the distance is 540 miles. Other flight or operating 
conditions may provide greater range; however, more fuel may be required 
in the interest of safety than is provided in the typical airplane. 

The possibility of burning part of the neutron shield as a means of 
providing emergency propulsion raises several points for discussion. It 
might be supposed that the after- shutdown radiation may become Intoler- 
able because of removal of the decalin. Calculations have shown, how- 
ever, that the decay of neutron flux is much more rapid than the decalin 
removal . 

Another item of interest is the penalty in payload if the entire 
neutron shield consists of decalin. It can be shown that the entire 
shield would weigh an additional 10,000 pounds, with a corresponding de- 
crease in payload. Nevertheless, the weight of decalin would then be- 
come 48,000 pounds. The extra amount of decalin would permit the air- 
plane to fly an additional 61 minutes at 0.5 Mach number at 30,000 feet. 
The emergency range would then become 935 statute miles . It should be 
mentioned here that the typical airplane already suffers a 5400-pound 
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penalty In payload compared with that available if the entire neutron 
shield is lithium hydride. 

As mentioned in the section BIOLOGICAL SHIELD, if all the decalin 
has been burned the reactor must not be restarted. It would be neces- 
sary to land the airplane and fill the decalin tanks with water (if no 
fuel is available) before starting up the reactor. 


FABRICATION AND DESIGN FEATURES 

In the designs described in this report it has been stated that 
certain items would be fabricated, without indicating the means of fab- 
rication. In addition, several items were not designed; it was merely 
stated that they were provided. The most important among the latter 
were several features of the control-rod design. No calculations were 
made of coolant-flow rates or control-rod temperatures nor was the radi- 
ator designed for cooling the low-temperature helium passing through the 
control rods. The amount of (Mo)£B 5 required in each rod also was not 
calculated. Statements were also made that cooling is provided around 
the pressure shell and through the biological shield, but no detailed 
calculations were made of these temperatures. None of these problems 
are considered serious enough to prevent the attainment of the reactor 
design proposed, and consequently for this type of study they are not 
considered in detail. 

In connection with the fuel elements, it was stated that they would 
be made of molybdenum and UC. Experience is lacking for mixtures of 
this nature. The following discussion is known to apply to UOg and Mo 
fuel elements, and it is expected that similar procedures would apply to 
the fuel elements used herein. The UOg and Mo, both in powder form, are 
mixed in proper proportions, blended, and cold-pressed in a die to form 
annular disks. The disks are then assembled to form tubes. The tubes 
are made in this way so that the fuel can be nonuniforml y distributed 
both axially and radially, since the disks vary in fuel quantity. In 
the mixture of UOg and Mo, there is a tendency for the Mo to react with 

the oxygen. The pressed disk is heated to about 3500° F in a Mo furnace 
with hydrogen atmosphere. The cladding may be put onto the meat by the 
use of molybdenum hexacarbonyl Mo(C0)g. The carbonyl is volatile at a 

low temperature. When passed over the fuel element, it decomposes, de- 
positing Mo on the meat and releasing CO. 

The stresses calculated for the fuel elements did not include pres- 
sures in the material caused by fission-product gases buildup. Uranium 
and its alloys tend to swell during irradiation because of such buildup 
(e.g., xenon and krypton). This tendency increases at high tenperatures . 
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Fission-product gases produce excessive pressures in metal already weak- 
ened by high temperature. 

UOg is free from the swelling and corrosion problems associated 
with uranium metal. When compacted to a high density, sintered, and 
irradiated at moderate temperatures, UOg can retain fission-product 
gases with no swelling at much higher burnups than metallic uranium 
elements . It is believed that extra space in the crystal lattice pro- 
vides for the fission products. It appears that UC combines the best 
characteristics of U and UOg- Furthermore, UC has excellent thermal 
conductivity and good irradiation stability. In the fabrication of the 
reactor, the fuel element proposed would first have to be developed to 
prove some of the expectations discussed. 

Lithium hydride was proposed as a shielding material, and some of 
its physical properties will be reviewed in this section to indicate 
its probable feasibility as a shielding material. Lithium hydride LiH 
comes as a crystalline solid, but is also available as a powder. It is 
in the latter form that it is proposed herein to use it. 

The finely divided powder is quite hazardous, since it is highly 
flammable if contacted with a small amount of water. In fact, powdered 
LiH ignites spontaneously upon exposure to air on a very humid day. A.t 
elevated temperatures, LiH reacts violently with oxygen and the halogens. 
In addition, LiH dust is irritating to nose, throat, and skin. 

The container proposed for the powdered LiH must be impervious to 
hydrogen and capable of withstanding relatively high pressure. The con- 
tainer must also resist attack of LiH at high temperatures, have high 
tensile strength, and be sealed hermetically. The tensile strength 
should be about 50,000 psi at 1200° F. Low-carbon stainless steels are 
suitable and easily available container materials. Satisfactory canning 
techniques for preparing hydrides for shielding purposes are still in 
the development stage. Glass, quartz, and enameled containers should 
not be used, since LiH has the ability to crack these materials at low 
temperatures . 

The fact that in the present reactor the LiH shield is out far from 
the core results in an environment in which heating and flux are not too 
high. In the present design, air-cooling has been provided for further 
reducing the heat to which the LiH is subjected. 


CONCLUDING REMARKS 

This study was made to determine the payload capacity and the de- 
sign feasibility of a helium-cooled nuclear-powered turboprop subsonic 
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airplane. A gross weight of 400,000 pounds, comparable to some airplanes 
in current production, was selected; and conservative design principles 
were used whenever possible. 

A 63,000-pound payload was obtained. This is about 16 percent of 
the gross weight and compares favorably with payloads of current conven- 
tional large airplanes. 

This nuclear airplane appears feasible, subject to experimental con- 
firmation of several components . The nuclear reactor will require the 
most research and development. Fuel-element fabrication and testing and 
control- rod fabrication, testing, and operation are of primary importance. 
Experimental data on the adequacy of the shielding are also required. 

In addition to reactor research and development, the heat exchangers and 
helium pumps require experimentation. The heat exchanger restricts the 
turbine-inlet temperature to 1800° R because of material limitations. 

Heat exchangers must be developed even for this low a value of turbine- 
inlet temperature. If the turbine-inlet temperature could be increased 
100° R, a 6000-pound increase in payload would result. Although the 
aerodynamic design of the helium pumps is conservative, the high rota- 
tional speed (65,800 rpm) will pose problems in bearing and shaft seal 
design. 

No problems are anticipated in the compressors and turbines because 
of their conservative designs. In addition, the weight and size of the 
airplane create no runway problems. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 10, 1959 
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APPENDIX A 


AIRPLANE AERODYNAMIC PERFORMANCE AND STRUCTURAL WEIGHT 

The equations and methods that were used to calculate lift-drag 
ratios and structural weights of the airplanes are presented herein. In 
calculating the performance and weight it is assumed that the following 
items are either known or specified: flight conditions (Mach number and 

altitude), gross weight of airplane, fuselage size and geometry, wing 
cross section and taper ratio (but not sweep or area), and maximum land- 
ing speed. A photograph of a model of the typical airplane is shown in 
figure 2. The lift-drag ratio of the airplane is calculated first. In 
this calculation the geometry and size of all the parts are completely 
specified. This information is then used to calculate the structural 
weights . 


A 


At 

si 

B 



Symbols 

The following symbols are used in appendix A: 

surface area factor; multiplier to define skin weight as a 
function of surface area 

horizontal plus vertical plan-form tail area, sq ft 
wing plan-form area, sq ft 
aspect ratio, b 2 /A w 

bending structural factor; multiplier to define bending 
weight as a function of load transmitted 

aerodynamic span, ft 

structural span (distance from wing tip to tip along mid- 
chord) , ft 

total airplane drag coefficient, based on wing plan-form area 
friction drag coefficient, based on surface area 
induced drag coefficient, based on wing plan-form area 
wing lift coefficient, based on wing plan-form area 
total airplane drag, qA^Cp, lb 
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<% 

equivalent fuselage diameter, ft 

f 3' f 4 

functions defined by equations (A23), (A24) , (A25 ) , and 
(A26), respectively 

K 

ratio of design to optimum lift coefficient 

K eq 

ratio of equipment weight to gross weight 

i K lg 

ratio of landing-gear weight to gross weight 

L 

wing lift, lb 

l F 

fuselage length, ft 

M 

flight Mach number 

M div 

drag-divergence Mach number (eq. (A5)) 

4 m 

sweep efficiency factor, induced drag without sweep divided 
by induced drag with sweep 

N 

1.5 

" 

n 

normal load factor, number of design g's divided by level- 
flight g of 1 

<1 

*1 g 

dynamic pressure, -g pv , lb/sq ft 

Re 

Reynolds number 

Sp 

fuselage surface area, sq ft 

V 

velocity, ft/sec 

% 

equipment weight that is independent of gross weight, lb 

W F 

fuselage structural weight, lb 

^F , avail 

available load in fuselage, lb 

W F ,1 

total weight of fuselage and contents and tail, lb 

W G 

airplane gross weight, lb 

W R,w 

total weight on wing but excluding wing weight, lb 
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weight of tail or empennage surfaces, lb 
W wing weight, lb 

W v ava -^p available load on wing, lb 

X ratio of specific heats for air 

A wing leading-edge sweep 

A taper ratio, ratio of tip chord to root chord 

p density, slugs/cu ft 

t wing thickness ratio 

Subscripts : 


F 

fuselage 

max 

maximum 

opt 

optimum 

w 

wing 


Lift-Drag Ratio 

The lift- drag ratio was calculated for the airplane without the 
powerplant or nacelles. The drag of these components is accounted for 
by subtracting their drag from the engine thrust. Thus the airplane 
drag, as hereafter referred to, does not include the engine or nacelle 
drag. 


The airplane lift-drag ratio is taken to be the lift of the wing 
divided by the sum of all the drags of the airplane (except engine and 
nacelle) . 

For subsonic flight the pressure or wave drag is assumed to be zero 
in all cases. For the wing, the sweep is set to assure that this is 
true. The airplane drag, therefore, is composed only of skin-friction 
drag and the induced drag (drag due to lift) of the wing. 

Skin- friction drag . - The skin- friction drag coefficient Cp is 

calculated from two different equations depending on whether the flow is 
laminar or turbulent. In either case the value 0.0005 is added to the 
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usual equations to account for skin roughness. When the Reynolds number 

Re is less than 2 x]_o® , the skin- friction drag coefficient is defined hy 
the following laminar- flow equation (including the term for roughness): 


C D,fr = + 0.0005 (Al) 

This equation, without the roughness correction term 0.0005, is given 
in reference 19. If Re is greater than 6x10®, Cp is defined by 
the following turbulent-flow equation: 


“T>,fr 


0.0306 


(Re) 1 /^ 


1 + h 2 W 1 


+ 0.0005 


(A2) 


This equation, minus the roughness correction, is taken from reference 
20 and is a flat-plate equation in which the temperature is the arith- 
metic average of the wall and free-stream temperatures. The Reynolds 
number is based on the fuselage length for the fuselage and on the root 
chord for the wing. For Reynolds numbers between 2xio® and 6x10®, com- 
bined laminar and turbulent boundary layers are considered. That portion 
or surface of the fuselage or wing from the forward part to the distance 
that results in a Reynolds number of 2x10® is taken as laminar with Reyn- 
olds number of 2x10®, and the remaining surface is assumed to have tur- 
bulent flow with Reynolds number based on the full length. With tapered 
wings (i.e., tip chord less than root chord), it was arbitrarily assumed 
that the same fraction of the wing surface was laminar for the entire 
wing as at the wing root, even though the root chord is longer than any 
other chord for tapered wings. 


Induced drag. - The induced drag or drag due to lift is found from 
the parameter Th e theoretical value of Cp j/c£ for subsonic 

airplanes is l/rt sf for elliptical wings and elliptical pressure dis- 
tribution. In practice, the induced drag is higher than this theoretical 
value. With appropriate values for the airplane type under consideration 
and some allowance for inefficiencies, the following equation was used: 





m \0. 95 its/ 


+ 


0.006 


') 


(A3) 


where the value of the sweep efficiency factor m is a function of the 
sweep angle A and is given by 


m = -y/cos A 


(A4) 


Declassified 1-24-1963 

M— WJLUJJ. 




Declassified 1-24-1963 


66 


&EUKEX 


The induced drag coefficient is, of course, found by multiplying 

c d,l/ c l by the square of the lift coefficient. 

Sweep angle . - The sweep angle of the wing is set so that the drag- 
divergence Mach number M^ y of the wing is equal to or higher than the 
flight Mach number. The drag-divergence Mach number is defined by 

M div = [ (0 ' ?5 " 0 ’ 1333 Cl) ( 1 + " ° ,0:L ] (A5) 

where 

N = 1 - ^ (A6) 

In order to provide some margin between M^ iv and the design flight 

Mach number M, the following equation was used to find the sweep 
angle A : 

(cos A - 0.02)" = r t| ■ - 1 d ; 02 ' ) [(0.75 - 0.1533 C L ) (l + - O.OlJ 

(A7) 


Tail drag . - The tail surfaces are assumed to have no lift, so that 
there is no induced drag of the tail. The skin- friction drag coefficient 
of the tail surfaces is assumed to be the same as that for the main wing, 
regardless of the Reynolds number. Thus, the friction drag of the main 

( 

wing is increased by the factor ^1 + — j. The tail area A^ was assumed 
to be 0.40 times the wing area A^. 


Lift-drag ratio . - The lift-drag ratio is given by the ratio of lift 
coefficient to the total drag coefficient: 


L _ £l 
E - Cj, 


(A8) 


where 

°D = C D,L + 2 ^- + A^ C D,fr,w + ^ C D,fr,F ( A9 ) 

and both and Cp are based on the wing plan-form area. For sim- 

plicity, the wing and tail surface areas are taken to be twice the 
plan-form area. The surface of the fuselage Sp is, for simplicity, 

taken to be the surface of a cylinder (excluding the ends) having values 
for length and diameter equal to those of the fuselage. This is almost 
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exactly true for a fuselage having a main cylindrical section and ellip- 
soidal ends each with length-diameter ratio of 2.5. 

Optimum wing area . - The wing lift is given hy 


The dynamic pressure q is determined hy the flight conditions (Mach 
number and altitude), and the lift L is equal to the airplane gross 
weight Wq. Thus, the product is specified by the gross weight 

and flight conditions. The specification of the wing area or alter- 
natively the lift coefficient is difficult; it is usually set to 

provide the best compromise of lift-drag ratio and structural weight of 
the airplane. The lift coefficient C^ is generally made larger than 
that for maximum lift-drag ratio, because this results in a smaller and 
therefore lighter wing. The optimum value of for best over- all per- 

formance may vary considerably, depending on the flight conditions and 
many other factors. At the same time the lift coefficient for maximum 
lift-drag ratio also varies, and therefore it was found convenient to 
specify the ratio of lift coefficient to the lift coefficient for 

maximum lift-drag ratio CL,opt • The value of .j. can be calcu- 

lated from known fuselage drag and flight conditions. Therefore, by 
specifying the ratio K, where 


the wing area can be determined. Various values of K can be tried 

to determine which value gives the best over-all performance. 

The wing area for a given airplane also determines the landing 
speed. Thus, if the landing speed is not to exceed a certain value, the 
wing area must be greater than a certain minimum value for a given air- 
plane gross weight. Consequently, the wing area is taken to be the 
larger of the values determined (l) by the maximum landing speed or (2) 
from the lift coefficient as determined by Q p^- and the specified 
value of K. ' 

The wing area determined from this lift coefficient is found from 
the following equations. The wing area is given by 


The lift coefficient is determined from equation (All) , where K is set 
arbitrarily by experience to give the best airplane performance. The 


L = qCjA^ 


(A10) 



(All) 
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optimum lift coefficient C L is defined as that value which gives 

the maximum lift-drag ratio. From equations (A3) , (A8) , and (A9) , 


c L 


D ■ c 2 


(A13) 


L L / 1 \ ( At\ Sp 

m" V0.95 73 + 0 -°° 6 ) + 2 \1 + A^) c D,fr,w + C D,fr,1 


Setting the derivative of L/D with respect to C^ equal to zero and 
solving for C^ Q p^. result in 


^L,opt 




(A14) 


s (o^h^ + °' 006 ) 

Combining equations (All), (A12), and (A14) and solving for \ yield 


Aw = 


' (§) I * 0 -°° 6 ) 

2 ( x + ^)°D,fr,w 


(A15) 


An iterative procedure is sometimes required, since m depends on 
C L or A w . In other cases, however, where the wing area is determined 

hy landing speed, this iteration is avoided. From equation (A12), the 
minimum wing area for a given weight is determined hy the maximum value 
of q (set hy landing speed at sea level) and the maximum value of 
(1.75). 


Structural Weights 

The total airplane gross weight is made up of (l) several known or 
assumed fixed weights, such as fixed equipment; (2) equipment weight that 
is proportional to or a function of the gross weight; (3) some weights 
not presently known hut to he determined later, such as powerplant or 
shield weight, which can only he determined after the drag or required 
thrust is known; (4) structural weights; and (5) payload. Methods for 
calculating the structural weights will he discussed in this section. 

The structural weights are, to some extent, dependent on the distri- 
bution of the loads in the airplane. Principally, removing loads from 
the fuselage and locating them on the wing, where the lift occurs, de- 
creases the bending moments in both the wing and the fuselage. Since the 
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load distribution is not always known initially, the structural weights 
for each airplane were first determined for a range of load 
distributions . 


The structural weights of both the wing and fuselage are assumed to 
consist of two parts : One part is a function of the surface area, and 

the other part is a function of the bending moments . The wing weight is 
given by 


4 


W w =1.5^+ 0.03 ~ + 2(10" 6 ) 


bgn(l + A) 


W, 


F ,1 


(A16) 


For the fuselage, 


w p = 3d| ^2 + 2.72(lO-°)na F ^j W Pjl 


-5> 


:V 


(A17) 


In both equations (A16) and (A17) the last term is a function of the 
bending moments. The term Wp ^ is the total load causing bending 
moments at the wing root and is equal to the total load in the fuselage 
plus the weight of the fuselage and tail. 

Equations (A16) and (A17) can be rewritten as 

K = fiOO + W y ,jf 2 (B„) (A18) 

and 

% = f s( A F) + W P,l f 4(®F) (A 19 ) 

where all the factors except are determined by the geometry, size, 

and load factor n. 


The weight of the empennage or tail is given by 


3/2 


W + = 7 — f-i 

t Aw 1 


K> + ^ (iq) 


(A20) 


The factor 


n - 1 
n 


accounts for the fact that in normal flight the 


tail was assumed to carry no load, so that during turns or climb when 
the angle of attack is above normal the lift coefficient on the hori- 

n - 1 


zontal tail surface is 


n 


times that of the main wing. 


The structural-weight equations presented are not available in any 
reference. There are very few sources that present equations that are 
useful in a study of this type. In general, they are either much too 


Declassified 1-24-1963 


ajse±oiT 


Declassified 1-24-1963 


70 


O-EiUIUSX 


detailed to be used in a study that must consider the complete airplane 
or so simple that they do not include the trends or variations that are 
desired when considering airplanes of different sizes, types, geometry, 
and so forth. Those equations that have been given in the literature 
and are suitable for a study of this type are generally in agreement as 
to trends but seldom in agreement regarding the magnitude of the trends. 
In the equations presented herein, an attempt was made to include any 
significant trend that might affect the weight and to obtain values that 
are in general agreement with other available sources. 

The structural weights and allowable loads can now be confuted in 
the following manner. All required values are either known from the 
calculation of lift-drag ratio or selected. The aerodynamic span b 
is given by 

b = vv? (A21) 


The structural span (distance from wing tip to tip along the midchord) 
b g is given by 



. _ b 

s cos A 

(A22) 

Also, 

4 

flW) = 1-5 A w + 0.03 3 

(A23> 


- b^n(l + A) 
f 2 (B w ) - 2(10- 6 ) ^ 

(A24) 


- «| ^ 

(A25) 

and 

f 4 (%) = 2.72(l0- 5 )nA F (^) 

(A26) 

Various values 

are then assumed for W R ^ w , the load that 

is placed on 


the wing but excluding the weight of the wing. 


Then the wing weight is given by 


fl(Aw) + (W 0 - W R Jf^) 
1 + f 2 (Bw) 
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and the total fuselage load Wf? , l 83 previously defined is 

«F,I = W Q - W E,w - M w ( A2e > 

Then the fuselage structural weight is 

W F = f 3 (A p ) + W y#1 f 4 (Bp) (A29) 

and the tail weight W^. (including the vertical fin) is given by equa- 
tion (A20) . 

All equipment is assumed to be carried in the fuselage so that the 
remaining load available in the fuselage for payload, reactor, shielding, 
or powerplant is 

avail - V F,J * »F " »eq - Vd - hg v F,l " *t < A3 °> 

It is assumed that the landing-gear weight (which is KjgWg) is 

divided between the fuselage and the wing in proportion to the gross 
weight of these components. Thus, the landing-gear weight in the fuse- 
lage is KjgWji and that on the wing is Kj g (W q - 

The load available on the wing, then, for powerplant or externally 
mounted reactor and shielding is 

H w, avail - «R,v ' K Jg< M G ' W F,J> < A31 > 

Values of W-p nv H .-n can be plotted as a function of W w ava il 88 
Wj^ w is varied. When the load on the wing av ail is known, the 

weight available in the fuselage W p ava j_i_ can te determined. 

Assumed Airplane Design Variables 

Where not otherwise noted, the following values were used in the 
airplane analysis: 

Fuselage fineness ratio, Z p / d p 10.0 

Aspect ratio, si 10.0 

Wing thickness ratio, t 0.15 

Taper ratio, A 0.30 

Maximum lift coefficient at landing, C-^ max 1.75 

Maximum landing speed, v-| ar ^ , ft/sec 175.00 

Normal load factor, n 2.0 

Landing-gear- to gross-weight ratio, K^g 0.075 

Equipment weight (independent of gross weight), W e( ^, lb .... 10,000 

Equipment weight factor, K g( ^ 0.030 
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In all cases investigated the wing area was determined by the land- 
ing speed. 
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ENGINE DESIGN AND WEIGHT 

The performance of the turboprop-engine type selected for this mis- 
sion was calculated by conventional methods with the aid of references 
6 and 21. Both engine weight and performance were calculated by methods 
which, as for the airplane, show the effects of changes in the design 
or operating conditions and give results that are in general agreement . 
with the existing technology. As with the airplane, rather conservative 
performance and adequately heavy components were assumed, thus covering 
cases where changes have to be made to existing equipment or new engines 
must be built with a minimum of development time and effort. 

Calculations were first made to determine the shaft power and net 
jet thrust for the assigned engine operating conditions, turbine-inlet 
temperature , compressor pressure ratio, and pressure losses. Sizes of 
the various engine components were next calculated for a range of engine 
airflows by assigning values to flow parameters of the compressor and 
turbine. Finally, engine weight and drag were calculated for the same 
range of engine sizes. The equations for calculating the engine weight 
are presented later. As described in appendix A, the engine drag was 
subtracted from the engine thrust instead of added to the airplane drag. 
In this way the engine size is not required when calculating aiiplane 
performance. The drag of the engine was calculated by the same methods 
as were used for the airplane fuselage. 

As with the airplane, most calculations were made for a standard 
set of conditions . Variations were then made in some of these condi- 
tions, varying them one at a time, to determine their effect on over-all 
airplane performance. Conditions typical of those that were varied are 
compressor pressure ratio, turbine-inlet temperature, and he at- exchanger 
and line pressure losses. Other design conditions, such as efficiencies 
and flow parameters, were held constant throughout the analysis, and 
these values are listed in table III. 


Symbols 

The following symbols are used in appendix B: 
d tip diameter, in. 

AH change in total enthalpy, Btu/lb 

l length, in. 
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^ 2^1 compressor pressure ratio 
HP shaft horsepower 

W weight, lb 

Subscripts : 

A accessories 

C compressor 

E engine 

g reduction gear 

m miscellaneous 

N exhaust nozzle 

nac nacelle 

st structure 

T turbine 


Equations Used 

In reference 22, empirical formulas are given for calculating the 
component weights of gas-turbine engines. These equations were modified 
to give slightly heavier components, and the resulting equations used in 
this analysis are given herein. 

Compressor: 


W c = 0.133 

where 



Turbine : 

W T = 0.0052(AH) T d^ 


(Bl) 

(B2) 

(B3) 
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Structure : 

W gt * 3 [50 + 0.1(W C + W T )J (B4) 

Accessories : 

W A = 125 + 0.02(W C + W T ) (B5) 

Exhaust nozzle: 

W R * d|^>.05 + 0.01 (B6) 

Reduction gear: 

W g = 0.1 HP (B7) 

Miscellaneous : 

W m = 0.2(W C + W T + W gt + W A ) (B8) 

The total engine weight W E is equal to the sum of these component 
weights. In addition, a nacelle weight was added, taken to be 

W nac = 0 ' 25W E (B9) 

The nacelle dimensions were assumed to be defined by the following 
equations : 


^nac ® 1*2 dp 

and 

^nac “ ^nac 


(BIO) 


(fill) 
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APPENDIX C 

HEAT EXCHANGER, HEADERS, AND LINES 
Heat-Exchanger Core Design 

A program was devised for an IBM 653 computer for designing heat 
exchangers of a variety of configurations encompassing parallel-flow and 
crossflow geometries. References 8, 9, and 23 served as the Bases for 
the he at -exchanger calculations. 

Symbols 

The following symbols are used in appendix C : 

A,B empirical constants 

c specific heat at constant pressure, Btu/ (lb) (9R) 

ir 

d diameter, ft 

d^ hydraulic diameter, ft 

f friction factor 

g acceleration due to gravity, ft/sec^ 

h heat- transfer coefficient, Btu/(sq ft)(sec)(°R) 

k thermal conductivity, Btu/ (sec) (ft) (°R) 

1 length, ft 

n number of passes 

P pressure, lb/sq ft 

Pr Prandtl number, Cp|j./k 

Re Reynolds number, pvd/p. 

S wetted surface area, sq ft 

St Stanton number, h/ pvc 

T temperature, °R 
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Tu number of heat-transfer units 

U over-all heat-transfer coefficient, Btu/(sq ft)(sec)(°R) 


V 

velocity, ft/sec 

w 

weight flow, lb/sec 

Z 

ratio of heat capacities 


over-all heat-exchanger effectiveness 


defined by eq. (C6) 

^,1 

defined by eq. (C7) 

M- 

absolute viscosity, lb/(ft)(sec) 

P 

density, lb/cu ft 

Subscripts : 

A 

first fluid 

B 

second fluid 

cir 

circular 

i 

inner 

max 

maximum 

min 

minimum 

n 

number of passes 

0 

outer 

1,2 

stations 

Superscripts : 

a,b,c 

empirical constants 
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Calculation Procedure 

The relation employed in calculating heat-transfer coefficients was 

(Re) (St) = a i? - e I a (Cl) 

(Pr)* 

Friction factors for pressure-drop calculations were determined from 

t?d 

f = B(Re) c (C2) 

In these two equations, A, B, a, b, and c are empirical constants 
applicable to a particular geometry. The quantities Re, Pr, and St 
are evaluated for average bulk or film properties depending on the type 
of data correlation used as source material for the empirical constants. 

For example, the data on compact heat exchangers in reference 9 are pre- 
sented in terms of the bulk properties of the heat-exchanger fluid. 

General method . - For this project, interest centered on multipass 
crossflow heat exchangers with the heating fluid, helium, contained in 
tubes. Air flows perpendicularly across the tubes, which are arranged 
to give multiple helium passes. Values of the empirical constants in 
equations (Cl) and (C2) were taken from reference 8 for the helium side 
and from reference 9 for the air side of the heat exchanger. 

Input to the heat-exchanger design calculations consisted of the 
following: 

(1) Inlet and outlet total temperatures of both gases 

(2) Inlet total pressures of both gases 

(3) Weight-flow rate of air 

(4) Allowable total -pres sure losses for both gases 

(5) Allowable tube wall stress 

(6) Tube and fin metal thermal conductivities 

(7) Geometrical data such as the spacings and diameter of the tubes 
and the pitch and thickness of fins 

The design calculations then provided the weights of the tubes and 
the fins, the number of tubes required, and the dimensions of the heat 
exchanger. Header and shell weights were determined in a separate 
calculation. 
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Calculation details * - Reference 9 presents a partly graphical 
method for calculating multipass crossflow heat exchangers. For use in 
the automatic computer, some of the design curves in reference 9 were 
approximated by polynomials, and these details may be of interest. 


The surface area of the exchanger is computed from 


(wc ) . 
n v p'min 

U 

The quantity (wc is the smaller of the two products of weight- 

flow rate and specific heat, and (Tu) n is the number of heat-transfer 

units for an n-pass exchanger. This latter quantity is a function of 
the over- all heat-exchanger effectiveness t)^., the number of passes n, 
and a ratio of heat capacities Z: 




( wc P^min 

( vc p) max 


Z < 1 


(C4) 


The procedure for evaluating (Tu) n is as follows, 
heat- exchanger effectiveness is given by 


The over-all 


^t * 


( wc p) A ( T A,l " T A,2) 

(wc) (T-n -i - T a t 
P min B > L A ^ 1 


T 


(C5) 


The calculation procedure then required that the over-all effectiveness 
be converted to an equivalent effectiveness of a single-pass heat 

exchanger : 



(C6) 


A relation between the effectiveness \ ± of a 
changer for which Z - 1 and the effectiveness 
heat exchanger for which Z f 1 


DLS 


single-pass heat ex- 
of a single-pass 


1.37 - 0.37 Z 


(C7) 
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Equation (C7) is an approximation determined from figure 5 in reference 
9. The desired value for the number of heat-transfer units (Tu) n was 
obtained from 


(Tu) n = n(l.4 t^-l + 14 T)i,i) (C8) 

This equation was also obtained from figure 5 in reference 9. The re- 
sulting value of (Tu) n permitted the calculation of heat-exchanger sur- 
face area by means of equation (C3). 

Helium and air properties . - Helium properties were calculated from 
the following relations based on reference 24: 

Absolute viscosity \i: 

p = 2 . 3072KL0 -7 (T)° * 647 lb/ (ft) (sec) (C9) 

Thermal conductivity k: 

k = 4. 085X10 -7 (I) 0 ' 650 Btu/ (sec) (ft) (°R) (CIO) 

Specific heat of helium was taken constant at 1.25 Btu/ (lb) (°R) . 


Air properties were based on reference 25 and were calculated from 
the following expressions: 

Absolute viscosity p: 


p = 1.672xl0" 7 (T) 0 - 684 lb/ (sec) (ft) 

(Cll) 

Thermal conductivity k: 


k = 2.055 x 10 _8 (T)°" 846 Btu/ (sec) (ft) (°R) 

(C12) 

Specific heat Cp : 


c = 0. 117 (T) 0,112 Btu/ (lb) (OR) 

or 

(C13) 

The temperature T in all the property equations is 
degrees Rankine. 

measured in 


Helium Lines 

As mentioned in the body of this report, concentric helium lines 
are used for transporting helium from the reactor to the engines and 
back again. The high- temperature helium flows through the inner lines, 
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and the low-temperature helium flows in an annulus surrounding the inner 
line. The method used for determining diameters of the lines will now 
he discussed. 


CO 

00 

w 

I 

w 


An expression is derived for the required diameters of the helium 
lines for specified conditions of length l , pressure drop AP, and 
weight-flow rate w. The derivation assumes an incompressible, constant- 
temperature flow, which is satisfactory for this situation, since the 
allowed pressure drop is small and the heat flow from the lines can be 
reduced to relatively low values. 

Circular cross section . - The equations necessary for determining 
the required inside diameter d of a circular pipe are as follows: 


i 

3 


AP = 4f - Efl 

a 2 g 

Jt ^2 

w = - d pv 


4f = 0.316(Re) 


^ ^ (eg. (131 ) , p. 81, ref. 7) 
Re = 

M- 


(C14) 

(C15) 

(C16) 

(C17) 


A combination of these equations yields for d: 

r i 

or 


,19 


■( 


4 ,*» 

0.16 \/ 4 f l 



d = 0.0115 




p AP 
4/19 


JJW 


1/19 W 7/19 


(CIS) 


(C19) 


Equation (C19) gives the diameter of a circular pipe required as a func- 
tion of weight-flow rate w, length l , density of fluid p, allowable 
pressure drop AP, and viscosity p. 

Annulus formed by two concentric circles . - In this case, the inner 
diameter of the annulus is t alien as a known quantity, and the outer 

diameter d_ is determined fr 


w 


the following equations : 


4f — p — 
d h 2 g 

(C20) 

i&l - a?)pv 

(C21) 
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4f 


= 0.316(Re) -1 / 4 


(C22) 



(C23) 


where d^, the hydraulic diameter for an annulus, is given by 

^ = ~ d i 

Equations (C20) to (C24) yield 


(C24) 



(C25) 


It is noted that the right side of equation (C25) is identical to the 
expression derived for the diameter of a circular pipe in equation (C18) . 
Thus, the hydraulic diameter d^ of an annular pipe can be related to 

the diameter of an equivalent circular pipe d c ^ r by the following 
equation: 


This relation is plotted in figure 60 for convenience in determining 
annulus dimensions. Given fluid properties, length of line, flow rate, 
and allowable pressure loss, the diameter of an equivalent circular 
pipe d c j- r is found from equation (C18) . Then for a specified inside 
diameter clj_, the required outside diameter cLq of the annulus can be 
determined with the aid of figure 60. 


Headers or manifolds are required to distribute helium to the heat- 
exchanger tubes and to collect the cooled helium into a single line for 
pumping. Both the inlet and outlet headers must have leak-proof connec- 
tions with 540 tubes and must be stressed to withstand a helium pressure 
of nearly 1200 psi. The greater design problem is, of course, associated 



(C26) 


In terms of outside diameter <3^, equation (C26) becomes 



(C27) 


Headers 
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with the inlet header, since the helium temperature at the inlet is 
1790° F compared with the helium outlet temperature of 790° F. The in- 
let header described here is one solution to this design problem. 

Figure 22 is a schematic sketch of the arrangement of the headers 
for the heat exchanger. The outlet header consists of three cylinders 
arranged for series flow. Helium flows to the pump and then to the in- 
let header, which is designed to allow the cool helium to flow around 
the inlet header. 

One purpose for this type of design is to allow the high-temperature 
parts of the inlet header to he nearly free of pressure loading. Pres- 
sure loads are transferred to the outer pipe, which operates at rela- 
tively low temperature • Actually, the pressure of the low-temperature 
helium is approximately 100 psi greater than the pressure of the high- 
temperature helium. 

In addition, a leak-proof seal is not required where the tubes are 
joined to the high-temperature inner shell, because leakage flow at these 
points would not cause any loss of helium from the system. Consequently, 
the design of a leak-proof system is eased, since only the low-temperature 
outer shell is involved. However, should a sizable leak develop in the 
inner shell, the mixing of the low- and high- temperature streams would 
reduce the cycle efficiency and power output. 


Declassified 1-24-1963 


BJSCtUSl' 


84 


Declassified 1-24-1963 

LJJiUXUil 


APPENDIX D 

AXIAL-FLOW HELIUM PUMP 

As mentioned in the tody of this report, either a centrifugal-flow 
or an axial-flow pump could he employed in the helium system; however, 
the axial-flow type was investigated for the typical airplane described 
herein. It was decided to use one helium pump for each engine. This 
appendix describes some features of the aerodynamic design of the pump. 

Symbols 

The following symbols axe used in appendix D: 
a sonic velocity, ft/sec 

c specific heat at constant pressure, Btu/(lb)(°R) 

ir 

D diffusion factor 

g acceleration due to gravity, ft/sec 2 

i incidence angle, deg 

J mechanical equivalent of heat, ft- lb/ (lb) (°R) 

P stagnation pressure, lb/sq ft 

r radius , ft 

T stagnation temperature, °R 

U blade speed, ft/ sec 

v flow velocity, ft/sec 

w weight flow, lb/sec 

(3 flow direction measured from axis of rotation, deg 

y ratio of specific heats 

8 p i/ p i,a 

■n polytropic efficiency 
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s *l/%i 

a solidity, ratio of chord to spacing 

cp flow coefficient 

>1/qi dimensionless temperature-rise parameter 

Subscripts : 
d design point 

h hub 

m midspan 

R rotor 

S stator 

t tip 

z axial component 

1 inlet 

2 outlet 
Superscript: 

1 relative to blade row 


Aerodynamic Design 

Specifications . - The helium pump was designed for the following 
specifications : 


Weight-flow rate, w, lb/sec 9.5 

Inlet stagnation pressure,, lb/sq ft 166,000 

Inlet stagnation temperature, T^, °R 1200 

Over-all stagnation-pressure ratio, Pg/Pq 1.085 

Expected polytropic efficiency, .... 0.80 

Design stagnation-temperature rise, AT, °F 50 


The low flow rate and high pressure require that the pump diameter 
be rather small. The sonic velocity in helium is 5000 feet per second 
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at 1200° R , and therefore the relative Mach numbers in the pump will be 
quite low for any practical blade speed and velocity diagram. Further- 
more, the flow in the pump will be effectively incompressible because of 
the low Mach numbers and pressure ratio. As a result, the aerodynamic 
design is straightforward, and conventional techniques and blade shapes 
can be utilized. 

Blade speed and pump size . - High blade speeds are required if few 
pump stages are desired. In this design, the blade speed is limited 
primarily by the allowable stresses . Other conditions that influence 
the choice of blade speed are the hub-tip radius ratio and the inlet 
velocity diagram of the pump. By compromising aerodynamic and mechan- 
ical considerations, the following parameters were selected: 

Blade tip speed, U^- , ft/sec 1500 

Hub-tip radius ratio, r^/r^ 0-8 

Inlet flow direction measured from axis of rotation, |3- 1_, deg .... 0 

Ratio of inlet flow velocity to inlet sonic velocity, v-^/a-L .... 0.1 

Inlet velocity, v^, ft/sec 500 

For these conditions, the centrifugal stresses at the roots of un- 
tapered blades are estimated at 40,000 psi. Since the temperature of 
the pump blades will not exceed 800° F, no problems are anticipated in 
the selection of blade materials. Low-creep-rate materials would be 
required for the pump blading to prevent excessive radial growth and 
rubbing of the blades on the pump casing. Type 347 stainless steel 
would.be suitable as a blade material. 

The pump tip diameter is calculated as 5.23 inches for the pre- 
scribed conditions, and the required rotational speed is 65,800 rpm. 

The angle of the inlet relative velocity is 71.55° at the rotor tip. 

One of the problems associated with small compressors for air is 
the occurrence of low Reynolds numbers and the resulting penalty in 
efficiency. The blade-chord Reynolds number for the helium pump consid- 
ered herein will be on the order of 10®, which is above the critical 
Reynolds number of 2*10® observed in compressor tests. Therefore, scale 
effects due to low Reynolds numbers are not expected in this helium pump. 

Velocity-diagram calculations . - Velocity diagrams are calculated 
for the mean radius only, because of the short blade height. The blades 
have a slight twist; the mean line essentially represents the complete 
blade. The following assumptions are made: 

(l) Equal temperature rise in each stage (ATqqi) 
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(2) Constant annulus area; therefore, equal axial velocities at 
exit from each blade row (incompressible flow) 

(3) Constant hub and tip radii 

(4) Ho inlet guide vanes ; axial flow leaving each set of stator 
blades 


These assumptions simplify the design to the extent that all stages are 
identical. There does not seem to be any reason to complicate the de- 
sign by using unlike stages other than to alter the range characteristics 
of the pump if this proves to be necessary. 

A dimensionless temperature-rise parameter typ is defined as 


_ gJc^ (di) 

U m 

Then, from equation (B7) of reference 26, tp can be related to the rela- 
tive flow direction at the rotor exit and the flow coefficient cp as 


and 


= 1 - <p tan |3£ 


cp = 




U, 


m 


tl 


m 


(D2) 


(D3) 


It is necessary to design the pump so that the blade loadings do 
not exceed the values that experience has shown to produce good effi- 
ciency. The diffusion factor D of reference 27 is used as a measure 
of the blade loading and can be written, for this simplified design 
situation, for the rotors as 


% 


= x _ V(i - » T > 2 + 


9 




'1 + cp‘ 


2a, 


R 


+ 9 2 


and for the stators as 



(D4) 


(D5) 


where Op and Og are the rotor and stator solidities for the midspan 
of the blade rows. 


Declassified 1-24-1963 




Declassified 1-24-1963 


88 


DJ!iUEU!U 


Typical airplane . - For the three-stage pump used in the typical 
airplane, the following values are taken: 

U t = 1500 ft/sec 

r h /r t =0.8 

v l/ a l ~ 0.1, a^ = 5000 ft/sec 
cp = 0.37 

AT = 50° R, AT st = 16.67° R 

ip = 0.286 

3£ = cot -1 cp = 69.7° 

1 - 

3 ’ = tan -1 * 62.6° 

Ci (p 

Turning angle in rotor, A3' = 7.1° 

Solidity, a = 1 
= 0.382 


-1 ^T o 

Flow direction entering stator 62 = tan — = 37.7 

(equivalent to turning angle in stators for this 
design) 


Dg = 0.515 

The values of diffusion factor for the rotors and stators are within 
the range for which high efficiency might he expected. They are not so 
low as to allow for a two-stage pump. 

Although little experimental data on small pumps are available, the 
assumed efficiency of 0.80 for the pump design is probably reasonable. 


Estimated Off-Design Performance of Helium Pump 

It was necessary to construct a performance map for the helium pump 
to allow for the calculation of powerplant performance at off-design 
conditions. Since the engine controls sure devised to hold constant 
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CM 

fH 

I 


d 


mechanical speed for the engine shafts the pump map was estimated for a 
constant actual blade speed. 


In order to simplify the task of estimating the off-design perform- 
ance of the helium pump, it was first hypothesized that only a small 
operating range (variation of incidence angles) would be required. Later 
calculations proved this hypothesis to be reasonable. Since the pump 
operates with low relative Mach numbers (<0.4), the limit on the oper- 
ating range allowed the assumption of constant polytropic efficiency. 

This is practically equivalent to an assumption of constant blade -element 
losses, as is the case in low-speed cascades for a limited range of in- 
cidence angles greater than and less than the optimum incidence angle. 

It was also assumed that the flow direction leaving each blade was con- 
stant, a valid assumption for a limited operating range. 

In order to match calculations for other portions of the powerplant, 
the pump map was computed using the following variables : 

(1) Ratio of flow coefficients — » — ^ 

°Pa w a 5 

T i 

(2) Inlet temperature ratio & = — 

x i,a 

(3) Inlet pressure ratio 5 = P^/P-^ ^ 

( 4 ) Over- all pressure ratio Pg/P-^ 

The following equations are require! for computing the pump map: 


an! 




(D6) 


(D7) 


(D8) 


The mean blade speed U m and the rotor-outlet relative flow angle Pg 
are constants and are equal to the design values. 
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The estimated pump performance map is shown in figure 61. The fig- 
ure also gives high and low limits of flow-coefficient ratio cp/cp^, 

which indicate the operating range for which the off-design calculation 
is expected to apply. These limits were determined from an estimation 
of the low-loss incidence- angle range (ref. 28) for the rotor blade 
elements . 
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APPENDIX E 


ANALYSIS OF TEMPERATURE DISTRIBUTION AND PRESSURE 

DROP IN FUEL ELEMENTS AND MODERATOR BLOCKS 

The reactor fuel element is circular and consists of a fuel "meat" 
g with inner and outer surfaces clad as shown in figure 25. The moderator 

w block has a circular inner surface and a hexagonal outer surface. The 

w coolant flows in the central hole of the fuel element and in the annulus 

formed by the inner surface of the moderator block and the outer surface 
of the fuel element. 

For this fuel element, heat will flow from the "meat" to channel 1 
and to channel 2 (see fig. 62). Because of flow differences in the 
channels, T^ and T a will be different and there will be some maximum 
temperature at radius in the "meat." If Tg and T-^, the 

coolant temperatures in the channels, are different, will vary axi- 
ally. For one set of conditions, call them on-design, radius e can 
be determined when radii a and d are known, so that T2 and T^ 
will be the same and Rjjj will remain fixed axially. The coolant is 
also considered to be flowing between the moderator blocks but at very 
slow speed. Its heat-transfer effect is neglected, and dT/dr at the 
equivalent radius f is assumed equal to zero. The equivalent radius 
f is the radius of a circle whose area is equivalent to that of the 
hexagon in a moderator block. 

The "on-design" method of analysis involves, as mentioned, the de- 
termination of the radius e so that the coolant temperatures are equal 
in both channels at a given axial station for a given total coolant flow. 
An equation for Rjjj is required, and then the coolant temperatures and 
fuel-element temperatures can be obtained. The design must also con- 
sider the coolant pressure drops allowable through the reactor and the 
amount of heat that must be picked up by the coolant from the reactor 
to provide the power for the engine. 

The "off-design" method of analysis takes a given fuel-element de- 
sign and determines the coolant conditions for a given total coolant 
flow through the fuel element. In this case the coolant temperatures in 
the channels at a given axial station will be practically the same, and 
the radius R m will also remain practically unchanged axially. 

The following analysis develops equations for use in both design 
methods, and further details on how to use them will be given after the 
equations are developed. 
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Symbols 

The following symbols are used in appendix E: 

A values given in table XIII to be applied to eqs. (E49) and 

(E50) , (cu ft) (sec) (°R)/Btu 

a radius of hole in fuel element, ft 

B values given in table XIII to be applied to eqs. (E49) to 

(E5l), (cu ft) (sec) (°R)/Btu 

b inner radius of fuel element, ft 

constants, Btu/(sec) (ft 5 ) , Btu/(sec) (ft 4 ) , and Btu/(sec) (cu 
ft), respectively 

c inner radius of outer cladding of fuel element, ft 

Cp specific heat of coolant, Btu/(lb)(°R) 

D values given in table XIII to be applied to eqs. (E49) and 

(E50), °R 

d outer radius of fuel element, ft 

hydraulic diameter in eqs. (E57), (E6l), and (E63), ft 

E values given in table XIII to be applied to eqs. (E49) and 

(E50), °R 

e inner radius of moderator block, ft 

F values given in table XIII to be applied to eq. (E51) , °R 

f equivalent radius for outer hexagonal surface of moderator 

block, ft (jtf^ = area of hexagonal block) 

f friction factor in eqs. (E6l) and (E62) 

G ratio of coolant mass-flow rate, w, to cross-sectional area 

of channel, lb/(sec)(sq ft) 

g gravitational constant, 32.2 ft/sec^ 
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h 

K en 

k 

Z 

N 

P 

<» 

^VV 

p 

Q 

Q 

<1 

Re 

r 

S 

T 

t 

U 

w 

X 

x 

Y 


convective heat-transfer coefficient, Btu/(sec)(sq ft)(°R) 

coefficient for entrance pressure loss into reactor passages, 
see eg. (E65) 

thermal conductivity, Btu/(sec) (ft) (°R) 
length , ft 

number of fuel elements 

total pressure of coolant, lb/sq ft 

power or heat picked up by coolant in reactor core and re- 
quired by powerplants, Mw 

static pressure, lb/sq ft 

volumetric heat source at any axial position x, 

Btu/(sec)(cu ft) 

average of axial heat sources, Btu/(sec)(cu ft) 
heat transferred radially, Btu/sec 

radius to point of maximum temperature in fuel-element meat, 
ft 

Reynolds number 

any radius in fuel-element - moderator-block assembly, ft 
surface area, sq ft 
total temperature, °R 
thickness , ft 

over-all heat-transfer coefficient, Btu/(sec)(sq ft)(°R) 

coolant-flow rate, lb/sec 

defined by eq. (E26), Btu/(sec)(sq ft) 

axial distance from core entrance, ft 

defined by eq. (E27), Btu/(sec)(sq ft) 
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a defined "by eq.. (E66) 

H viscosity of coolant, lb/(ft)(sec) 

p density of coolant, Ib/cu ft 

Subscripts : 
a radius a 

av average 


B 

b 

c 

cl 

d 

e 

eff 


bulk 
radius b 
radius c 
cladding 
radius d 
radius e 
effective 


F film 

f radius f 

in into reactor upstream end reflector 

m radius 

mod moderator 

mt meat 

out out of reactor downstream end reflector 

r reflector 

w wall 

1 hole 

2 annulus 
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Required Data and Assumptions 

The following data are assumed known before making a fuel-element 
design: 

(1) Total pressure of coolant flowing into reactor, P^ n 

(2) Total pressure of coolant flowing out of reactor, P QU -t 

(3) Total temperature of coolant flowing into reactor, T^ n 

(4) Total temperature of coolant flowing out of reactor, T ou ^. 

(5) Total coolant flow through reactor, w 

(6) Power required by powerplant including pumping power or power 
required from reactor to heat coolant (does not include power 
dissipation from reactor core to reflector, thermal shields, and 
pressure shell), ^ 

Irlr 

(7) Dimensions of fuel element (radii a, b, c, and d) 

(8) Ratio of moderator area to frontal area 

(9) Reactor-core length, l 

(10) Reflector thickness, t r 

The following general assumptions are mader 

(1) There is constant radial heat distribution in the core. There- 
fore, each fuel element picks up the same amount of heat, and only one 
fuel element need by analyzed. 

(2) There is no heat transfer between moderator blocks. Therefore, 
as mentioned before, dT/dr at radius f is zero, and all heat from the 
moderator flows to the coolant in the annulus. 

(3) 90 Percent of total reactor heat is released "locally” in the 
fuel element. This heat originates from fission fragments and beta par- 
ticles (see ref. 15, p. 638). 

(4) 5 Percent of total reactor heat is released "nonlocally" in 
core or in moderator. This heat derives from gammas and neutrons (see 
ref. 15). 

(5) 5 Percent of total reactor heat is released "nonlocally" in re- 
flector, thermal shields, and pressure shell. This heat originates 
mostly from gamma radiation (see ref. 15). These values apply to 
"thermal" reactors but are assumed for the present reactor. 

Declassified 1-24-1963 


96 


Declassified 1-24-1963 

SBUKJST 


(6) No internal sources of heat from nuclear reactions exist in 
the coolant. 


Heat-Transfer Analysis 


The development of the heat- transfer equations is similar to that 
in two special problems of reference 29 (pp. 130-138). Some repetition 
of that work is included herein in order to obtain continuity in the 
various analyses of the present project. The development could also 
have proceeded from heat-transfer equations given in reference 15 (pp. 
657 and 662) . 

"Meat" temperatures . - For r > (see fig. 62), the heat gener- 
ated between 'R m and r in distance dx is 

dq = J* Q ftit dx 2rtr dr (El) 

^m 

where is the volumetric heat source in the "meat" at station x 

from the entrance. 

For the assumption that is constant radially, equation (El) 

integrates to 

dq = *Qmt( r2 " ( R 2) 

By Fourier's law, 

d( i = -Vt 2ltr % ( E3 ) 


Equating equations (E2) and (E3) and integrating from radius 
to radius c give 

At the outer "boundary (from eq. (E2)), 

dq c = rtQ mt^ c2 ’ B to) dx ( E5 ) 

When r < R, a corresponding development will result in the following 
equations : 


%t / o R m R 4 - V \ 

- % - $ - " - 3 ) 

(E6) 

d % = " fc 2 ) 3 * 

(B7) 

OT?n"mi?n 
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Cladding temperatures . - The cladding Is very thin, and the assump- 
tion is made that heat generation in it can be neglected. Then the heat 
from the meat passes through the cladding by a pure conduction process. 
Also, because of the thinness of the cladding, the assumption is made 
that the conduction of heat through it can be represented by the 
equation 

q = -kS ^ (E8) 


and for S the mean area will be used. Then, for the outer cladding 
the difference between the temperatures at radii c and d becomes, 
using equation (E8) , 


- T* = 


(d - c)(dq c /dx) 

2 ”(° + ^ T ^) 1 


(E9) 


2 I cl 

The difference between the temperatures at radii a and b becomes 

(t - a)(dq b /dx) 

x b “ x a ~ 


2n| b 




(E10) 


2 }\l 


Surface-to-coolant temperature differences . - The differences in 
temperature between the fuel-element surfaces and the coolant are 


dq c /dx 


T d “ T 2 " 2rtd he 


(Ell) 


and 


d ^b/ dx 


T - Tt = 
x a X 1 2rta h-, 


(E12) 


where h^ and are the convective heat-transfer coefficients of 

the channels. 


Temperature change from maximum meat temperature to coolant . - The 
temperature difference from radius Rm to the coolants in channels 1 
and 2 can be obtained by combining equations (E6), (E7), (E10), and 
(E12) and equations (E4) , (E5), (E9), and (Ell), respectively. The re- 
sults are as follows : 


T m - 


T 2 = 


2k mt \ 2 



Vt^ 2 - R m> 

2U 2 d 


(E13) 


QTTP'DTjrn 
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Qmt 


T m - T i = 2kmt 




^mt^m " ^ ) 
2U x a 


(E14) 



(E15) 


(El 6) 


Determination of radius R m . - By subtracting equation (E14) from 
equation (E13) and solving the resulting equation for R m , the following 
equation is obtained: 



(E17) 


Because the cladding is so thin, c « d and b ** a; therefore, 



(E18) 


Moderator temperatures . - On the basis of the assumption that all 
the heat in the moderator flows to the coolant in the annulus, and there- 
fore dT/dr =0 at radius f, the temperature difference between radii 
f and e can be represented by equation (E6) with changes appropriate 
to the moderator block: 


■ m = 
f “ e 2k mod 




(E19) 
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In a manner analogous to getting dq-^, dq^ becomes 

d % = «w ( f2 - e2 ) dx ( E2 °) 

The temperature drop from the surface at radius e to the coolant 

dq^/dx 


T 1 _ m — 

e 2 2ne h. 


Combining these equations gives 


T m ^od 


( 2 


2f 2 in X + e 2 - f 2 ) + 9m ° d ^ e ~ 6 1 


(E21) 


(E22) 


Axial temperature variations . - Prom equations (Ell), (E12), (E13), 
or (E14) , (E21), and (E22), the maximum and surface temperatures of the 
fuel element or the moderator block, respectively, can be obtained for 
any axial position x. Before solving the equations, the coolant tem- 
peratures at the various axial positions must be determined. For any 
position x the coolant temperatures depend on the equations 

w 2 c p| T 2 ( x ) " T 2 ,inl = f d < 3 c + / ( E23 ) 

L J O 0 


WiC p [ti(x) - Ti^ n J - J 


d % 


(E24) 


Using equations (Ell), (E12), (E13), (E2l) , (E22) , (E5), (E7), 
(E20) , (E23) , and (E24) , it can be shown that 


, v *(c 2 - Rm) 

T Jx) = To - + — — X 

2 ' # 2,m w 2 c p 


, «(f 2 - e 2 ) y 

w 2 c p 


(E25) 


where 


and 


x = f 0 
Y - f 

o 


(E26) 

(E27) 


e ramnan 
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Also, 


, . *(r£ - 1 2 ) 

*iW-®L,ln+ WlCp ^ 

(E28) 

m , , „ , i . W*X K 1 - 

T a (x) = Tl (x) + ^ 

(E29) 

T d (x) = T 2 (x) + ^ ^ " 

(E30) 

• e2) 

(E31) 


[4 


T m( x ) = T 2< x ) + W x ) 

T f (x) = To(x) + Qm °^ . X) 


'^ ln £ 

mt \ d 


(c 2 - bS) 


2Ugd 


1 


|_W 


2 f e 2 - f 2 -p2 _2 

r In - + — + 1 e 


^^mod 


h 2 e 


3 


(E32) 

(E33) 


The next step in the solution of the equations is to determine the 
axial heat-source variation in the reactor. 


Axial distribution of heat source . - For a uniform axial fuel dis- 
tribution, reference 15 (p. 643) shows that the heat distribution in a 
cylindrical unreflected reactor follows a cosine law. Because the re- 
actor has end reflectors, a modified power distribution is used. This is 


W*) = «mt(|) c °=(|f - |) 


(E34) 


where Q^(Z / 2) is the maximum flux that occurs at the center of the core. 

This distribution is called the " chopped" 2/3-cosine distribution. The 
flux is not zero at the core-reflector interface. 


Now, the average heat source in the core equals 


Vt “ 


jf Sit (i)=° s (ir - 1) 6 * 


(E35) 
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Integrating and solving for £^.(1/2) result in 


Snt(!) - i- 21 Srt 

Then 

^ ^(ff 2 - §) 

(E36) 

(E37) 

From a similar development. 


- 1-21 Va C0B (lr - 1) 

(E38) 

Combining equations (E26) and (E37), and (E27) and (E38), 
following equations for X and Y: 

gives the 

x - vih(i? - f) + 4 \ 

(E39) 

Y w H[ sln (lr - 1) + #] 

(E40) 


A constant or uniform axial distribution of fuel, which generally 
results in a heat-source distribution approximately like the one just 
derived, will cause the temperatures of the reactor core to vary in a 
manner that will be quite different from that desired when the allowable 
stress curve is considered. It would be advantageous to design reactors 
so that the temperature resulting would be more nearly like those desired 
from allowable stress considerations. The following is an analysis of 
the heat-source distribution required for any arbitrarily selected tem- 
perature distribution (which could be selected on the basis of the stress 
desired). From the heat-source distribution so determined, the fuel dis- 
tribution would then be determined, if it is assumed that the neutron- 
flux distribution is the same as for the case of uniform fuel 
distribution. 

Suppose the heat-source distribution in fuel element or moderator 
to be represented by an equation of the following type: 

Q(x) = + C^x + C- 2 j (E4l) 

Then 

Q(0) = C 3 (E42) 

\ 2 

«(f) " C 1 T + C 2 \ + 0 3 < E43 > 

Q(Z) = C^ 2 + C 2 Z + C 3 (E44) 
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Also, 


Therefore, 



(E45) 


(E46) 


(E47 ) 


(E48) 


Any three values of the fuel-element or moderator temperature can 
now be set; for example, T(o) , T(— J, and T(l) can be selected at any 

one of the radii a, d, e, m, or f. Suppose is chosen. Then 

three equations based on equation (E32) are available with known temper- 
atures into which equation (E25), for the appropriate axial position, is 
first substituted; and then into the three resulting equations the equa- 
tions (E26), (E27), (E42), (E43>, (E44) , (E46), (E47), and (E48) are sub- 
stituted; again for the appropriate axial position. Three equations re- 
sult with three unknowns , the constants Cq, Cg, and Cj in the heat- 
distribution equation (E4l). These equations can be solved for the con- 
stants; and thus, for the chosen temperature distribution, the required 
heat-source distribution is known. In like manner, Cq, Cg, and C3 can 
be determined for any set of temperatures T a , T^, and so forth. Of 
course only one set can be chosen, and from the resulting Q(x) varia- 
tion the other temperatures will automatically result. By performing 
the algebraic manipulations mentioned, the equations for Cq, Cg, and 
C3 resulting from choosing T m , T a , T^, T e , or T f variations are 


C 


1 


C 


2 



2 

l 


4(A + 2B)D - (A + 4B) E 
(A + 6B)(A + 2Bj - 2(A + 3B)(A + 4B) 

(A + 6B)E - 8(A + 3B)D 
(A + 6B) (A + 2B) - 2(A + 3B)(A + 4B) 


c 3 = 


* - T in 


B 


(E49) 


(E50) 


(E5l) 
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where A, B, D, E, and F have the values given in table XIII depending 
upon which temperature distribution is chosen as the desired one. 

The foregoing analysis depends on the selection of one set of tem- 
peratures at x = 0, x = Z/2, and x = l, but it does not indicate what 
values to choose or limitations on values that can be chosen. From 
equation (E48) it can be shown that 

Q = C ± + C 2 + C 3 (E52) 


and, as will be shown in the following section, and ^ have 

fixed values depending on the power required of the reactor. Therefore, 
temperatures must be chosen so that, when Cj_, C 2 , and C3 are calcu- 
lated and substituted into equation (E52), the average volumetric heat 
source will not be exceeded. Temperatures can be picked at random until 
such a result occurs, but some guidance to the selection of temperatures 
would save computation. A method of doing this would be to substitute 
into equation (E52) for the temperature selected as known, values of cor- 
responding C^, C 2> and C3 obtained from equations (E49), (E50), and 
(E5l). Then one equation results with Q in it, which is known. (The 
method of determining Q is discussed later.) For example, for getting 
(x) for known values, the resulting equation is 


„ T m(!) <W - K 2 T ln ,, 

^°T + K 3 T m (0) ' 4 


(E53) 


where the K are known terms. From the allowable stress curve , some 
idea of the ratio of temperatures can be obtained, and then T^(o) can 
be obtained from the preceding equation. Knowing T m (o), then, and the 
temperature ratios, enough information is available to get the (^.(x) 
variation from the analysis presented. The same type of equation will 
result for the other temperatures listed in table XIII. 


Average volumetric heat source . - According to the assumptions out- 
lined previously, 90 percent of the total reactor power or heat is in 
the fuel elements and 5 percent in the moderator blocks. Also, 5 percent 
is lost to the outside of the core. If the power in megawatts required 
by the powerplant is denoted by then the total heat produced by the 

reactor on the basis of the percentages quoted is 948.1 0.95 Btu per 

second. The heat given up by each fuel-element - moderator-block assem- 
bly Is then 

948.1 <? 

(”l + “2>y T out - T in) * - S PP <E54) 
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where w^ and wg are coolant-flow rates through one assembly. The 
heat given off by a fuel element will then be 

0.9(w^ + w 2^ c p^out ” •'•in) 

__ 


and by a moderator block, 

0 . 05 ( W;L + w 2 )c p (T Qut - T in ) 

0.95 

and the average volumetric heat sources consequently are 

- _ 0.9 ( W 1 + w 2) c p( T out " T in) 

^ 0.95 «z( c 2 _ b 2 ) 

7T 0.05 ( W 1 + w 2) c p("'out “ •'•in) 

W-0.95 ^(fZ.eZ) 


(E55) 


(E56) 


Heat-transfer coefficient . - The convective heat-transfer coeffi- 
cients in the two channels of the fuel-element - moderator-block assembly 
are required in the solution of the equations. Average turbulent-flow 
forced-convection coefficients have been experimentally obtained over a 
wide range of Reynolds numbers, surface temperatures, heat fluxes, en- 
trance conditions, and l/d^ ratios, as reported in reference 30. Eval- 
uating the properties of the fluid and density at the film temperature 
gave good correlations. The empirical equation resulting was 


h = 


0.0216 k F 


hi 




(E57) 


where d h is the hydraulic diameter. For the hole in the fuel element, 
dft = 2a; and for the annulus, d^ = 2(e - d). The subscript B refers to 
the bulk temperature of the fluid, assumed equal to the arithmetic mean 
of T -^ and T QU ^., and F refers to film temperature conditions . The 

film temperature is 


T f = 


_ t B + T w,eff 
2 


(E58) 


where T^ e ff is the effective wall temperature of the entire channel. 
The equation for it is 
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(E59) 


where e in this equation is the base for Napierian logarithms, and S 
is equal to 2jtaZ for the hole in the fuel element and 2n (e + d)Z for 
the annulus around the fuel element. For the coolant used in this proj- 
ect (helium), Cp is 1.25 Btu/(lb) (°R) , and the viscosity and conduc- 
tivity (obtained from ref. 24) are 


The determination of h requires iteration, because it depends on 
'E V)e ff> which in turn is calculated for a known value of h. 


The decrease in pressure of the coolant through the end reflectors 
and the fuel-element - moderator-block assemblies in the core was kept 
below a specified minimum value. The calculation of the coolant pres- 
sure change through two channels of the assemblies (the end reflectors 
having the same geometry as the core but with no uranium in the circular 
elements) was thus required. It was assumed that there was no heat 
pickup by the coolant in the end reflectors and thus no change in total 
temperature in passing through the reflectors. 

The static-pressure drops through the channels were calculated 
using normal momentum and friction pressure-drop formulas for turbulent 
flow through pipes. Or, for momentum pressure drop, 


H p = 2.3072(T F )°' 647 xi0 -7 lb/(ft)(sec) 
k p = 4.085(T p ) 0,650x 10" 7 Btu/ (sec) (ft) (°R) 


Pressure-Drop Analysis 



(E60) 


and for friction pressure drop, 



(E61) 
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where p av is the arithmetic mean of p in and P ou f The friction 
coefficients used were those given by Reynolds' analogy as 


| = 0.023(Re)"°’ 2 


(E62) 


which are in good agreement with Karman-Nikiiradse values for the Reyn- 
olds number range most frequently used. The correlation of friction 
coefficients also depends, as in the case of heat-transfer coefficients, 
on the heat flux and the temperature at which the fluid properties are 
evaluated. Making use of this fact, and combining equations (E6l) and 
(E62), 


Ap 


0.184 


T- 


B l 


G c 


few 

UipTp 


0.2 Tp d^ 2gp, 


av 


(E63) 


For the holes in the fuel elements, the hydraulic diameter d h in equa- 
tion (E63) is equal to 2a; and for the annulus between the fuel element 
and moderator block, it is 2(e - d) . When the end-reflector passages 
are being considered, l in equation (E63) is the end-reflector thick- 
ness t ; and, when the core passages are considered, l is the core 
length . 


Contraction and expansion losses into and out of the end-reflector 
passages were also considered. The change in total pressure from the 
inlet duct to the passage inlets was calculated from 


^in * ^passage inlet “ K en 


G £ 


( n 


(E64) 


where K en was obtained from reference 9 and is given by the equation 


K. 


en 


fRe 

= 0.5 — 


inlet passage 
3000 


3 


0.077 


-0.4a 


where 


a = 


+ (e 2 - d 2 ) 


(E65) 


(E66) 


The change in total pressure from the passage outlets to the outlet duct 
was obtained from 


^passage outlet " *out (l " a ) 


G 


(p ) 2 g 

passage outlet 


(E67 ) 
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Each part of the pressure-drop calculation requires iteration, be- 
cause the densities used to get the Ap or AP values depend on the 
downstream pressures being calculated. 


General Method of Design 

There is some radius e, as brought out at the beginning of this 
analysis, for which will remain constant axially and Tq(x) will 

equal T-^x) at each axial position. By assuming some value of radius 
a and values of the coolant-flow rates (knowing total coolant flow 
through a fuel element and the ratio of w^/wg obtained as explained 
in the following sentence), various values of e can be assumed until 
such a condition occurs. An approximate formula for e to start the 
calculation can be determined by equating equations (E25) and (E28), 
solving for the flow ratio ^q/wg from the resulting equation, then 

equating this ratio to a 2 V&/[( e2 - d 2 )-/ (e~- d)] . This expression 
was obtained by assuming that the only pressure drops in the two chan- 
nels were those due to friction, and of course must be equal. As a con- 
sequence the flow ratio is roughly proportional to the ratio of the 
cross-sectional areas of the channels multiplied by the ratio of the 
square root of the hydraulic diameters of the channels. After determin- 
ing the pressure drops, the temperatures are calculated and iteration is 
required because the pressure drops are dependent on the surface temper- 
atures of the fuel element and moderator block. 

The calculations are made for various assumed values of radii a 
and e and flow rates Wq and Wg until the conditions of (l) the 

same coolant temperature rise in both channels, (2) the same over-all 
pressure drop in both channels with a value below the required limit, 
and (3) the required hardware surface temperature limitations are met. 
Because of the numerous iterative processes involved, computations were 
carried out on an IBM 653 computing machine. 

For the present study the fuel elements and moderator blocks were 
designed by the foregoing method for a given set of conditions. These 
conditions were later changed, but no change was made in the individual 
fuel elements and moderator blocks. More fuel elements and moderator 
blocks were added to the core, thus changing the flow rates. A new IBM 
653 program was made to calculate this "off-design" case, which was also 
briefly discussed at the beginning of this analysis. The problem is to 
find the resulting fuel-element and moderator-block temperatures and, 
for a given total amount of coolant flowing through each assembly, the 
amount flowing through each channel that gives the same pressure drop. 

It can be shown that, because the ratio of a 2 -/it/[( e2 - d 2 )- v /(e - d)] 
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is the same as for the on-design case, then w ]_/ w 2 i s about the same; 
and thus using equations (E25) and (E28) the temperature of the coolant 
out of the annulus and hole will he about the same as for the on-design 
case. Thus, a very good estimate of w-j_ and W2 can be made when the 
sum of the two values for the off-design case is known. 


OPHTJ MI! 


Declassified 1-24-1963 


E-288 


Declassified 1-24-1963 


QJb^XUH-L 


109 


APPENDIX F 


ANALYSIS OF STRESSES IN FUEL ELEMENT AND MODERATOR BLOCK 

A theory for the stresses in a circular fuel element is available 
in reference 31 (pp. 169 to 182) for conditions not quite the same as 
for the fuel element used here. One assumption in the reference is that 
no external or surface forces are exerted on the heated tube. Another 
assumption is that there is no heat flow across the outer boundary of 
the fuel element and therefore the maximum temperature occurs at this ra- 
dius. In the present fuel element the heat is dissipated to the coolant 
at the inner and outer radii of the element, and the temperature reaches 
a maximum (or no heat flow occurs) at some radius within the fuel-element 
"meat." Some initial equations of reference 31 could be used, and the 
reader is referred to the reference for their development. These equa- 
tions are then further developed in the following analysis for the par- 
ticular conditions mentioned in the foregoing discussion. 


A 

A' 

a 

B 


B' 


b 

C' 


Symbols 

The following symbols are used in appendix Fr 
integration constant, (Btu) (in. )/sec 
integration constant (see eq. (F36) ) 
inner radius of fuel element, in. 
integration constant, Btu/ (in. ) (sec ) 
integration constant (see eq. (F37)) 
outer radius of inner cladding of fuel element, in. 
integration constant (see eq. (F42)) 
integration constant (see eq. (F3)), °R 
integration constant (see eq. (F6) ) , °R 
inner radius or outer cladding of fuel element, in. 
specific heat of coolant, Btu/(lb)(°R) 
outer radius of fuel element, in. 
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E 

e 

F z (i) 

F 3 (i) 

f 

k 

l 

m 

P 

Q 

Q 

B 


r 

T 

W 

W 1 

w 2 


x 

a 


modulus of elasticity of material in tube, lb/sq in. 

inner radius of moderator block, in. 

defined by eq. (F29) 

defined by eq. (F30) 

defined by eq. (F3l) 

equivalent outer radius of moderator block, in. (jtf^ = area 
of the block with hexagonal boundary) 

thermal conductivity, Btu/(sec) (°R) (in. ) 

length of core of reactor, in. 

= (R m - a) or (R m - R), in. 

static pressure of coolant, lb/sq in. 

volumetric heat source in material at a radial and an axial 
point in tube, Btu/(sec)(cu in.) 

equivalent average volumetric heat source for solid tube of 
radius equal to outer diameter of tube at an axial position 
x, Btu/(sec)(cu in.), see eq. (F25) and its derivation 

= a, in. (used for convenience in expanding a^( a/d) equation) 

radius at which maximum tube temperature occurs, in. (deter- 
mined from temperature analysis of appendix E) 

any radius in tube, in. 

total temperature, °R 

defined by equation (F32), lb/sq in. 

coolant flow through hole of fuel element, Ib/sec 

coolant-flow rate through annulus of fuel-element - moderator- 
block assembly, lb/sec 

axial position from coolant inlet to reactor core, in. 

coefficient of thermal expansion of material in tube, 
in. /(in.) (°R) 


td 
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3 = (b - a) or (b - R), In. 

r =0/^ 

6 = (d - a) or (d - R) , in. 

e longitudinal strain, in. /in. 

v Poisson 1 s ratio 

t = r/d 

p = (c - a) or (c - R) , in. 

<r stress, psi 

ex equivalent stress, psi (see eq. (P52)) 

Subscripts : 
ad adjusted 

d radius d 

eq equivalent 

in into reactor core 

mod moderator 

mt meat 

o reference 

out out of reactor core 

r radial 

z axial 

6 tangential 

1,2,3 used with r and £ when double or triple integrals are 

used 


srorRwp 
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Fuel Element 

Development of general stress equations . - The fuel element and the 
sketch used for the analysis are shown in figures 25 and 62. The max- 
imum temperature occurs at radius Rj^ in figure 62. Then 

KU. ' 0 (F1) 

From reference 31 (p. 170) for steady-state conditions for a tube, 

ls( r i) + l = ° < F2 > 

where Q is the volumetric heat-source strength, which varies with 
location. 

Integrating equation (F2) once, 


5*|/ r « r i ar i (C i=“ 


(F3) 


Using equation (Fl) in equation (F3) , 


c i - - 1 to i 


Combining equations (F3) and (F4) , 

dT 1 f T _ . _ 

r d? + k 1 ^*1^1=0 

% 

Integrating equation (F5), 

1 f T 4rx /* r 1 

T+ U Qr 

At r = d, 

1 r d drq /* r l 

T d + k 7 — i Qr 2 dr 2 + C 2 = 0 


2 dr 2 + C 2 - 0 


(F4) 


(F5) 


(F6) 


(F7) 
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Subtracting equation (F7) 

from equation (F6), 


T - T d + 

1 /* r drq 

i / — / Qr 2 dr 2 = 0 

k *4 r i x 22 

(F8) 

Let 


*-S 

(F9) 

and 


it 

(F10) 


where Q q is the power or heat-source density at some reference point. 
Then equation (F8) becomes 




(Fll) 


The fuel element is made of three rings, two cladding rings and the 
"meat" ring containing the uranium. For stress analysis, the assumption 
was made that the fuel element consisted of one material only, or the 
cladding and matrix could be represented by one ring. If the element 
had been divided into two cladding rings and one "meat" ring, each would 
have to be treated separately. Because of the complexity of the anal- 
ysis, it was considered beyond the scope of the present study. There- 
fore, the modulus of elasticity E, the coefficient of thermal expansion 
a, and Poisson's ratio v were considered constant in the radial direc- 
tion. For these conditions, reference 31 (p. 177) arrives at the fol- 
lowing equation for radial stress in the tube: 


_d_ / 3 dP r \ Ear 2 dT 
dry dryi-vdr - 


(FI 2) 


Substituting in this equation the value of dT / dr 


from equation (F5), 


d (_3 dCT r\ Ea f T „ 

dr dr ) ~ (l - v)k r J Qr l ^1 ^ F13 ^ 

This is a differential equation for the radial thermal stress in terms 
of properties of the tube and the distributed heat source Q. It is 
to he integrated twice subject to appropriate boundary conditions. 
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Integrating it once, including an arbitrary constant that depends on 
the slope of the stress curve at r = d, 


d- q r _ Ea 
dr ~ k(l - v) 


^ (/ r l dr l I 1 + 

m 


2A 


(F14) 


Integrating again, including a second constant that depends on the 
stress at r = d. 


a. = 


Ea 


* k(l - v) 


l-C [J / 1 f 2 7 2 ^ S] + ■ ?) + B } 


From reference 31 (p. 177), the tangential stress is 

dcr 

°e = r w + J r 

Substituting equations (F14) and (F15) into this equation, 

•r / 


(F15) 


(F16) 


Oa = 


e “ k(i - v) 




Ir^ dr 


1 f 1 ^2 ^2 

X ) 


r 2 ^2 


f 2 ^3 ^3 

Rm 


+ A (? + 7 ) + B } 


(F17) 


Also from reference 31 (p. 176), the axial stress may be written 

0 Z = v(a r + O 0 ) - Ea(T - T 0 ) + e z E (F18) 

where T Q is a reference temperature that is eliminated in the analysis. 
Substituting equations (F15) and (F17) into equation (F18), 

* k(i E - -vT {7 j[ 1 ^2 


+ 2 


■f fr -f ' 


< /* r 2 
r 2 dr 2 J Qr 3 dr 3 

V. Rm 


) + fi + 2B }- E [ ai 


(T - T q ) - e z ] 
(F19) 


S5iK 5 RkT 

Declassified 1-24-1963 


E-288 


Declassified 1-24-1963 




115 


From equation (F8) , 


i r T f * 1 


(F20) 


Putting equation (F20) into (F19) gives 
Eav 


a z " k(l - v) 


{? -C f 1 ^ c ** s 

2 / r [J f 2 dr2 ^ ** **)] 

jT / ri ^2 + ft + ' *[“<% - T °> - «.](«) 


Rm 

Define, as in reference 31, the following: 

Fl(5) - -C [f ( 


1 r«i r J i l h ai2 rfs « 3 | 


P2(s) - 7 -f ( ii d?i £\ rss ii2 ) 
' S(,> ■ ^ ^ ^ 


(F22) 


(F23) 


(F24) 


As "before, let 


and 


I = r/d 

r = Q/Qo 


(F9) 

(F10) 


The average power density Q for the space occupied by the tube 
can be related to the reference value Qq by the following integral 
(the same as in ref. 31) : 
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Q«d 2 = 


r 


Q2jtr dr = 


/ 


r V* 


r dr 


Q 


- / 




or 




Q = -2Qo ./ rli d| x 


(P25) 


Using equations (F22) to (F25) in equations (F15), (F17), and (F2l), the 
latter three equations become the same as those in reference 31 except 
for the limits of integration in equations (F22 ) , (F23) , and (F24) . The 
final equations are 


o r U) = - 


EcoQd 2 

ra/d 

2k(l - v) J ri d| 


[f-lU) + A' (l - + B*j (F26) 


*«(0 - - 


- 2 
EooQd 


/* a 7^ - 

2k(i - v) J n 


as 


[jqd) + F 2 (6) + A» (l + + B'J (F27) 


EoQd 2 v 


°z^) = ~ ~ /*a/d 


2k(l 


-V)/' 


rl d| 


X ^(l) + F 2 (0 + F 3 (0 + 2A' + 2B’ + C'J(F28) 


where 


A' = 




B' = 


B 

Qo d2 


and 
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E [a(l a - T c ) - «J 

^ = 

EoQvd^ 

a/d 

n ± d6 x 

Solution of general stress equations . - The first step in the solu- 
tion of the radial, tangential, and longitudinal thermal stress equations 
is to integrate the F}_(|), FgCl); and F 3 (|) equations. The results are 




d| 


+ -K f ri di + rl 3 d! (F29) 

Kjj/d - 2 * 


4r "i 


>(l) = i f % rl di - f ri di - n 3 d£ (fso) 

2 \/d 2 | 2 \/d 


P 3 (fi) t— 7 - In I f 1 Tl d| - f l rl in | d| 

1 ' " V 1 


(F31) 


Now the boundary conditions are 


Md = " p 


° r (i) - -p 

where p is the coolant static pressure. 
Let 


W = 


EccQd 


/ ojd 

rl 


d& 


(F32) 


Then applying the houndary conditions to equation (F26) , 

-p = -W |Fi It J + A' (l 


- w [ f i(!) + a ’( 1 -f) + B ] 


(F33) 


REfTRM 1 

Declassified 1-24-1963 


Declassified 1-24-1963 


118 


DJ&UKCil' 



-p = -W ^(1) + A'(l - 1) + bJ] 

(F34) 

From equation 

(F29) , 




(F35) 

Using this in equation (F34) and solving the resulting equation and equa' 
tion (F33) for A' and B', the result is 


A ' d 2 - a 2 Pl ( d ) 

(F36) 


ii 

PQ 

(F37 ) 


In the cladding It Is assumed, that there is no heat source; there- 
fore, t - 0. In the meat, the heat source is assumed constant radially 
and equal to Qq (at any axial station), so that y = 1. Then, in the 
stress equations. 



(F38) 


The remaining terms required before the stress equations can be solved 
for a particular radial position are C 1 and F-^a/d), the latter 
being required in equation (F36). 


Using equation (F29) and the variation of 
shown that 


T just noted, it can be 



The third boundary condition is that the total force exerted over 
a cross section normal to the axis is zero (there are no end restraints 
and no acceleration), or 



a z 2jtr dr = 0 
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or 


2*& c 




a/d 


2*d £ 


/ 


,/d 


o z l d| - o 


(F40) 


Using equations (F28) and (F32) in equation (F40) , 


-2«Wd 2 f + F 2 U) + F 3 (e) + 2A' + 2B’ + C']| d| = 0 

*Vd L J 


or 


C' » -2(A' + B') - 


2d £ 


d 2 - a 2 ‘'a/d 


/ [sFqU) + F 2 (0 + F 3 (|)]i d| (F41) 


Using equations (F29) to (F3l) and knowing the variation of X pre- 
viously noted, the expression under the integral sign can he integrated 
so that 


C’ = -2(A' + B') 


2d 2 fl *1 - , i M 2 /* 2 - 

d 2 - a 2 L 8v d 2 0 v W V a 2 , 


1 C 4 - b 4 1_ 

+ 16v d 4 " 8v 


+ 


1_ 

4v 





(F42) 


All terms are now available to solve for the stresses at any radial 
location. The maximum combined or equivalent stress will occur at the 
inner radius. Therefore, equations for the stresses at radii a and 
d were determined. 


Stresses at inner and outer radii . - Examination of the stress equa- 
tions shows that, to get the stresses at radii a and d, F-|_(a/d), 

F-^(l), F 2 (a/d), F 2 (l), FjCa/d), and Fj(l) must be known. The first two 

have been obtained (eqs. (F35) and (F39)). From equations (F30) and 
(F3l) and the variation of y, it can be shown that 
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Equivalent stress . - The radial, tangential, and longitudinal 
stresses are combined into an "equivalent" stress. The shape of the 
element being considered does not enter into the equation for equivalent 
stress, because a point is being considered when the stresses are com- 
bined. The equation is 

% = ^2 V( a r ~ a e )2 + (a r - a z) 2 + C°e - CT z) 2 ( F52 ) 

The equivalent stress is determined at the inner and outer boundaries of 
the fuel element. / 

Adjusted equi valent stress . - The foregoing equivalent stress is 
greater than the actual stress in the material, because some plastic 
flow and stress-relieving occurs. A method based on the short-time 
stress-strain curve of the material as shown in figure 63 is used to get 
the actual stress to compare with allowable stress values for the mate- 
rial. The short-time curve is plotted as shown in figure 63, and the 
elastic part of the curve is extended as shown by the dashed line. If 
the equivalent stress is in the dashed region of the curve, as shown, 
then a vertical line is dropped from this point to a point on the solid 
curve (point a). The stress at this point, the adjusted stress a a q> is 

the one used to compare with the allowable stress . If a g( ^ falls on 
the line Ob, it is not adjusted but is used to compare with the allow- 
able stress. 


Moderator Block 

The foregoing analysis is applicable to the moderator block, if the 
hexagonal outer boundary is considered to be represented by an equivalent 
circle of radius f . The maximum temperature in the moderator block 
occurs at this radius f, and the block is uniformly heated radially. 

In a sense, it is equivalent to the meat of the fuel element (only one 
ring), and f would be 1 throughout from radius e to radius f. Then 
the following changes are required in the equations to make them appli- 
cable to the moderator: Change to f, a to e, b to e, c to f, 

and d to f . The pressure of the coolant around the moderator block 
is again equal to p, but no cooling is assumed, as mentioned in 
appendix E. 


Determination of Q 

In the foregoing analysis, the term Qq is equal to {^(x) or 
Qmod( x ) appendix E, depending on whether the fuel element or moder- 
ator block is considered. Combining equations (E37) and (E55) and 
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equations (E38) and (E56), the following values for and Q^^Cx) 

result t 

W x ) = wh 1 * 21 (Vl - — — — — I 1 - T±n) cosfe - (F53) 

0-95 *z(c 2 - b 2 ) V 3Z V 


( W 1 + w 2 ) c p(^out “ ^in) / 2 itx £ 

31 ' 3 


„ , % 0.05 , V W 1 ^ w 2/ (J p^ x out - x in^ (\ 

W x ) = 0^5 1 - 21 _ e 2 } cos \: 


(F54) 


M 


Then for the fuel-element stress equations, 

•a/d 

'1 

Integrating this equation and substituting equation (F53) result in 


/ a/ a ✓'b/d 

n d| - - 2 ^t(x) / I d| 
c/d 


(F55) 


Q = 0.365 


(w x + w 2 )c p (T out - T ln ) 

( 2 jtx rt 

— ^ COSIttt— - — 


Jd 


/ 2 jtx rt\ 

\3Z " 3/ 


(F56) 


In like manner, for the moderator-block stress equations, 

/*e/f -e/f 

« - - 2 Snod (x > / « « " - 2Q moa( x ) ./ ^ « 


(F57) 


or 


Q = 0-02025 


(vi + w 2 )c (T out - T in ) 


If 1 


( 2 «x 

cos V— 


(F58) 


Accuracy of Calculated Stresses 

Equations (F48) to (F5l) and similar ones for the moderator block 
are of a form such that the W term is very large and the bracketed or 
braced term consists of differences between very small numbers. These 
latter numbers must be carried out to many significant figures for any 
accuracy to result from the calculations. 

In the present calculations the results were carried to eight sig- 
nificant figures on a computing machine ; but even then some uncertainty 
arose as to whether this was accurate enough. As a consequence, one 
formula (that for the fuel-element tangential stress CX 0 ( a/d)) was ex- 
panded by substitution of series expansions for the log terms and use 
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of the inner radius plus small increments for the radii of the fuel ele- 
ment to get the formula in a form more accurate to use than the original 
one. Then a check of the stress using the expanded form was made and 
compared with the original calculations described. The following shows 
the expansion of the Cg(a/d) formula of equation (P48). 


CO 

CO 

CM 

w 


Expansion of Tangential-Stress Formula 

Substituting equations (F32) , (F38) , (F39) , and (F43) into equation 
(F48) , the equation for Oq( a/d), and reducing give the following: 


M _ EaQ d c 

0 \d/ ” 8k(l - v) ( d 2 _ a 2 } (c 2 


T15) [ 4cV ln G) 


• ^ m ) - «> 2a2 + «& 2 m 

■ + 2a 2 R 2 - 2R^d 2 - o 4 + b 4 + 2b 2 d 2 j - 


(F59) 


Letting a=R, b=R+0, R = R + m, c = R + p, d = R + 5, and expand- 
ing the log terms in a series expansion, 

ln (f) - * §) 

" 5 - l(s f + l(s) 5 - l(l) 4 + I® 5 - + • • • ( F6 °) 

which converges for p < R. Similar equations are obtained for In (b/a) 
and In (d/a), the p terms in equation (F60) being replaced by 3 and 
6, respectively. 

Substituting these expansion formulas for the log terms and the new 
nomenclature for the radii in equation (F59), and in turn expanding terms 
such as (R + 0) 4 , substituted for b 4 , and finally reducing give 
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- 8k(l - vK^)^ - 7) H" 82(P ' ^ + 8<P2 ' ?2) 

- 3 (p 3 - P 3 )] + 4|1 | & 3 (p - m) - 5 2 (p 2 - m 2 ) + S 2 (p 2 - p 2 ) 

+ | 5(p 3 - 3 3 ) - i (p 4 - p 4 )] 

+ ~ [ & 4 ( P - m) - 2& 3 (p 2 - m 2 ) + 26 2 (p 3 - p 3 ) - 6(p 4 - p 4 ) 

+ | (p 5 - p 5 )j + | S 5 (p - m) + 6 4 (p 2 - m 2 ) 

- S 2 (p 4 - p 4 ) + | 6(p 5 - p 5 ) - i (p 6 - p 6 >J 

+ — ~ [& 6 (p - m) - 6 5 (p 2 - m 2 ) + S 2 (p 5 - p 5 ) 

15R- 5 L 

- 8(p 6 - P 6 ) + | (p 7 - |J 7 )] + . . . j - p (F6X) 

This type of equation is more amenable to computation and does not 
require the large number of significant figures for accuracy needed by 
the original equation. 


Sources of Properties Used in Stress Calculations 

The properties needed in the calculation of the fuel-element and 
moderator stresses are the modulus of elasticity E, the coefficient of 
thermal expansion a, Poisson's ratio v, and the thermal conductivity 
k, of the materials used for the parts. The fuel elements, as mentioned 
in the body of this report, are considered to be made entirely of molyb- 
denum. The moderator blocks are made of beryllium oxide BeO. The sources 
for the data mentioned for the two materials are given in table XIV. 
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APPENDIX G 


ANALYSIS OP NUCLEAR CRITICAL REQUIREMENTS 

The following analysis gives a method for determining the uranium 
investment required for a given cylindrical reactor to make it critical 
or to have a required reactivity with and without control rod, distrib- 
uted boron, and xenon buildup after shutdown effects. The one- 
dimensional, two-group diffusion method of reference 12 (p. 238) is ex- 
tended to six groups. A fully reflected cylindrical reactor is approx- 
imated by first prescribing values of the axial buckling and solving in 
a radial direction, and then by prescribing values of radial buckling 
based on assumed flux levels, solving in an axial direction, and iter- 
ating. The method is based on dividing the energy spectrum into six 
groups and the reactor into three regions, first in a radial direction 
and then in an axial direction. The core is one region, and the re- 
flector is divided into two regions. It is necessary to obtain the 
neutron cross section and diffusion coefficients for each energy group 
before using the diffusion- theory equations. The averaging procedure 
is similar to that in reference 12 (pp. 227-229). Some details of the 
method are given herein. 


Symbols 

The following symbols are used in appendix Gl 

2 2 
B "buckling, cm 

d reactor-core diameter, in. 

E neutron energy, ev 

E Q energy of source neutrons produced by fission =10 ev 

k ff effective multiplication factor 

1 length of reactor core, in. 

N concentration of reactor constituent, atoms/cc 

& reactor power, Mw 

r reactor-core radius, d/2, in. 

T average temperature of reactor core, °K 
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E , 


u 

V 


X 

\r 

2 


cr 

9 


lethargy of neutrons, In — 
neutron velocity, cm/sec 

any variable distance used in buckling formula, cm 
transport mean free path, cm 
macroscopic neutron cross section, Ncr, cm"^ 
microscopic neutron cross section, cm 2 / atom 
neutron flux, neutrons/ (cm 2 ) (sec) 


i 


Subscripts j 


a 

f 


n 


th 


r 

s 

t 

th 


absorption 
fission 
refers to n w “ neutron group 
radial 
scattering 
total 
thermal 
axial 


1 , 2 , 3 , 

4 , 5,6 


denotes the six energy groups into which neutrons are divided 
(6™ group is thermal group) 


Superscript : 

average when used with 2 and cp 


Detailed Procedure for Determining Required Uranium Investment 

Reactor composition . - The number of nuclei per cubic centimeter 
N of each material in the reactor is first determined. Enriched uranium 
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with 92.5 percent and 7.5 percent was assumed. The methods 

of determining the N values are straightforward and will not be dis- 
cussed. The N of the coolant was determined at an average pressure 
and temperature of the coolant. 

Microscopic neutron cross sections . - The microscopic cross sec- 
tions a & , o s , and cr^ of the materials in the reactor are required as 

well as Cf of the uranium. They were obtained from reference 32, over 
the range of lethargy u from 0 to thermal. The thermal values of the 
l/v absorbers were averaged for the Maxwell-Bolt zmann distribution by 
multiplying them by -/j t/2. For uranium, the thermal values were also 
multiplied by a factor to account for the non-l/v behavior of the ab- 
sorption and fission cross sections. The thermal energy used to deter- 
mine the thermal cross sections was calculated from 

E th = 8.61><10- 5 T, ev (Gl) 

Average macroscopic cross sections and diffusion coefficients . - 
The average macroscopic cross sections and diffusion coefficients for 
each of the six groups are required as mentioned before. These include 
values for both the reflector and the core. The general method of get- 
ting these constants using the N and a values described is reported 
in reference 12. The macroscopic cross sections were flux- weighted. 

The buckling of the cylinder was computed using estimated values of 
the axial and radial reflector savings. This buckling could, for ex- 
ample, represent an equivalent bare sphere for purposes of discussion of 
the results of the constants program, and also leads to an estimate of 
the neutron leakage from the reactor core. Modifications included in 
the program from the one mentioned in the reference are computations of 
the neutron age by Marshak* s formula for mixtures of heavy elements (ref. 
33) and use of Goertzel-Selengut slowing down for hydrogen and Wigner 
slowing down for other elements (ref. 34) rather than the simple Fermi 
approximation . 

Determination of multiplication factor, k e ^. - Knowing the con- 
stants mentioned, the diffusion theory can be used to determine, for a 
given reactor and uranium investment, the effective multiplication factor 
k eff' The details of the method for the six-group, three-region case 
used herein are similar to those fully described for a nine-group four- 
region reactor in reference 35. In the present case, the method was pro- 
grammed for an IBM 653 computer. The restrictions used in the six- group 
theory are as follows : 

(1) All fissions occur in the first region, the core. 

(2) All fissions are assumed to occur in the lowei’ three energy 
groups . 
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(3) Neutrons may be born only in the upper three energy groups. 

The first region, the core, is divided into 20 parts; and the sec- 
ond and third regions, in the reflector, are divided into 10 and 8 parts, 
respectively. The output of the program consists of the flux and the 
power distributions for each group at the 38 points just described and 

k eff • 


Determination of k e ff for fully reflected reactor. - The program 
discussed in the previous section can be applied to a cylindrical reactor 
either radially or axially. As pointed out in the beginning of this 
appendix, in solving in a given direction, the leakage in the other di- 
rection is approximated by using an estimated buckling in that direction. 
The formulas used for the buckling are 



E 

[- 


2 1 

(l + x) + 1.42 AtJ 

2.405 

(r + x) +0.71 


2 


2 


(G2) 


(03) 


2 

where B z is an estimated axial buckling used with a radial solution, 

O 

and B£ is an estimated radial buckling used with an axial solution. 

The term that must be estimated is the reflector savings (not necessar- 
ily the same value in the two equations). This puts some reflector on 
the core in the direction in which the assumed buckling is being deter- 
mined. In addition, the actual thickness of the reflector on the side 
or end of the core is put into the programs discussed when either a 
radial or axial solution, respectively, is being sought. The transport 
mean free path A^ in the equations for buckling is a value determined 
in the constants program for the core. The value of for the third 

group is used as being about an average for all of the groups. 


The method for getting k e f£ for the fully reflected reactor is 
to run the constants program using the B^ for the equivalent spherical 
reactor previously described. Then an estimated B z is calculated using 
equation (G2), and k e ff is found using the six-group program in a 

radial direction. From the results of the program it is possible to set 
up a series of coefficients in a matrix, in a manner similar to that 
described for the nine-group solution of reference 35, and by Crouts 

2 

rule to determine manually the B^ resulting from the calculations. 
Knowing B^ and the estimated Br, B£ is calculated from B^ - B z - 
This is then used and the program is solved axially. Again from these 
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2 2 2 
results B^. Is determined manually and B£ is calculated from B^. 

minus the B^ used. This is compared with the original B^ estimated, 

and if it is different, the procedure is repeated. The fact that the 

B^ used in the constants program may not be the same, as it theoreti- 

O 

cally should be, as the final B£ determined as described is usually of 
little importance, because the constants for most reactors of the type 
involved vary very little with B B . 

The procedure just described will yield k e ff for the chosen re- 
actor; but, for machine calculations, it is awkward in that alter each 

o 9 

machine calculation with an estimated or Br a manual calculation 

2 z r 

of B£ is required instead of leaving the calculation on the machine. 

As a consequence, a method has been worked out that avoids this awkward- 
ness. In addition to running a constants program with B B for the 
equivalent spherical reactor, at least two more are run with widely dif- 
fering values of B^. The multiplication factors are then obtained from 
the programs and are plotted as shown in figure 64(a). The point marked 
A on this plot is the one for which the buckling is that of the equiva- 
lent spherical reactor. Then using the constants from the program in 
which the equivalent spherical reactor buckling is used, multiplication 
factors for a radial solution of the diffusion- theory program are ob- 
tained using three estimated values of axial buckling. This procedure 
is repeated for an axial solution using three estimated values of radial 
buckling. The results are plotted as shown in figure 64(b). 

p 2 

There is only one B^ (= b£ + B^) that results in a given k e ff 

for the reactor. So B+ is the same on either the B„ varied curve or 
2 ^ r~\ ,z 

varied curve for a k e ff such as (l) in figure 64(b). Call this 

B t,® • Then ^ 

B i© = B ? + <d) (°*> 

2 

for the B z varied curve. In like manner for the other curve , 

B t,®= B r,© +B z ( a5 > 

2 

Because the left sides of the two equations are identical, B r in the 

2 2 

first equation equals B r ^ and in the second equation equals 

B z • Then, for a given on the figure, adding the intercepts 

of the constant line with the curve, as shown, will give a point 

on a curve of k e ff against b|. This procedure is repeated for other 
values of k e fp, which results in the curve shown in figure 64(c). 
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The curve just obtained does not give the answer for the reactor 
specified hut gives a curve for many reactors with the chosen core com- 
position. There is some equivalent hare reactor with the chosen core 

o 

composition whose Bf and are the same as that for the specified 

reactor. If the results for such a hare reactor were known , the point 
could he spotted on figure 64(c) , and the k eff for the specified re- 
actor would he known. It is not known, so the approach is to take a 
curve such as in figure 64(a) and plot it on the same sheet with the 
curve in figure 64(c). The results would he like that shown in figure 
64(d) . It will he remembered that the results at the of point A 

in the figure were used to get the results for figure 64(c). Now, if 
point A had coincided with point B in figure 64(d), the ^ e ff cross- 
over point B would he the answer for the specified reactor. The fact 
that the points do not coincide mak.es the B an approximate 

answer. If the fast cross sections are constant over the range of B^ 
used in the constants program, then, even if A does not coincide with 
B in figure 64(d), the k e ff at the crossover point E is satisfactory. 

2 

If the constants mentioned do vary some with B^, the calculations are 
repeated until A and B do coincide, or the constants used in the diffu- 

o 

sion theory are at the B£ that results at the crossover point In fig- 
ure 64(d). 


Average flux in core and percent fissions in each group . - Part of 
the output of the diffusion-theory program, as mentioned previous ly, is 
the distribution of flux for each group through the core and reflector. 
The flux distribution in the reactor is the product of the radial and 
axial distributions. The fluxes obtained are of course relative fluxes. 
The absolute fluxes were determined as follows : For a core diameter of 

d inches, a core length of l inches, and a reactor power of 0* mega- 
watts, the following formula is applicable under the restrictions imposed 
before (fissions occur only in lower three energy groups): 


^xio 6 x5.ixiq 10 

l(2.54) 3 



CPu 


+ 2 


f 5=“ + 
’ ^6 



(G6) 


Therefore, the left side gives the number of fissions per second per 
cubic centimeter of core. The right side gives the same thing. The 
fluxes cp^, cpg, and so forth, determined as described, are relative 

fluxes to some standard of reference. Dividing one by the other as 
shown in the formula gives the average flux of one group relative to the 
other. The cpg term outside the parentheses is the absolute flux for 
the thermal group. Since everything is known in the formula, except this 
value, the equation can be solved for it. The average fission macroscopic 


Hw Tr j lit 1 \ 
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cross sections 2^ for the three groups in the formula are obtained 
from the programs described before. 


Knowing one absolute average flux, in this case cpg, the absolute 


fluxes for the other groups can be obtained by multiplying this value 
by the relative flux ratios cp-^/cpg, cpg/cpg, and so forth, obtained as 
described. The percent fissions in the fourth, fifth, and thermal groups 
are calculated from the following formula: 


(Percent fissions ) n 


Zf,n Pn 


E 2 

n=4 


'f,n 


100 



(G7) 


where for the fourth group n = 4, and so forth. The fluxes are the 
absolute values. 


The details of analyses that use results obtained by methods de- 
scribed herein to determine xenon effects, amount of distributed boron 
required at startup, results of burnout of boron and fuel, and control- 
rod effects are given in the body of this report. 
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APPENDIX H 


ANALYSIS OF SIDE-REFLECTOR , THERMAL-SHIELD, AND 
PRESSURE -SHELL TEMPERATURES 

It Is assumed that heating is caused principally by primary gammas; 
and, although radial attenuation of the gammas occurs in each reactor 
part considered, the attenuation is small enough in each part that a 
constant heat source (an average of the heat source in the part consid- 
ered) is assumed in the analysis of the part. The heat source is assumed 
constant axially in any slab also. The radii of all parts are large 
enough that the equations for heat flow through slabs are considered 
valid. Thermal radiation between parts will be neglected. Attenuation 
of the radiation through the coolant (helium in this project) is also 
neglected. 

The analysis has been made assuming that the coolant flows through 
the passages will be metered so that the coolant temperature in each 
passage at a given axial position from the inlet will be equal. This 
will result in a fixed ratio of coolant weight flows in adjacent pas- 
sages throughout the length of passage, and the maximum temperature in 
a part will remain at the same radial position for all axial positions. 
The particular equations for the shields are developed herein. The fixed 
position of the maximum temperature simplifies the analysis. 


Symbols 

The following symbols are used in appendix H: 

A cross-sectional area of shield normal to gamma radiation, sq in. 

c specific heat of coolant, Btu/(lb)(°R) 

Jr 

d reactor-core diameter, in. 

h convective heat-transfer coefficient, Btu/(sec) (°R) (sq in.) 

I heat flux, Btu/(sec)(sq in.) 

k thermal conductivity, Btu/(in. ) (sec) (°R) 

l reactor-core length, in. 

& reactor power, Mw 
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Q volumetric heat source, Btu/(sec)(cu in.) 

Q average volumetric heat source, Btu/(sec)(cu in.) 

r radius, in. 

T temperature, °R 

t thickness, in. 

w coolant-flow rate, lb/sec 

x distance from inner radius of a shield, in. (see fig. 65) 

Ay incremental axial distance from coolant entrance in passage, in. 

(see fig. 65) 

P a absorption coefficient for gamma radiation, in."’*' 

Subscripts : 

A slab A 

B slab B 

C slab C 

con container around reactor core 
i inner 

in coolant into reactor 

max maximum 

o when used with I, denotes core flux; when used with r, denotes 

outer radius (away from core) 

r. . inner radius of slab A 

A, i 

i’ . inner radius of slab B 

B, i 

t total 

w,eff effective wall 

x any position x in a shield as shown in fig. 65 
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1 passage 1; when used with T only, refers to coolant temperature 

in passage 1 at any axial position; when used with T and 
another subscript , refers to temperature on surface near 
passage 

2 passage 2 with same other connotation as subscript 1 

3 passage 3 


General Equations 

Temperatures . - In figure 65 are shown three slab cross sections 
that are representative of a series of thermal shields. Coolant flows 
in the passages between the slabs , and heat flows from a slab to the 
coolant on either side. Equations for the coolant temperatures and the 
temperatures for the middle slab, a representative one, will be 
developed. 


In a small increment of coolant-passage length Ay (see fig- 65), 
the following equations are applicable with the maximum temperature 
occurring at x^, Xg, and x^ : 


Tl(Ay) = T in + 


T 2 (Ay) = T ln + 


[ 2nr A,o Q A (t A - X A } + 


W l c p 


[^BjO^B " x b) + 2nr C,i^ x Cpy 

_ 


where T(Ay) is the temperature Ay distant from the inlet, and 
hp [t b ^ x (Ay ) - T x (Ay)] (Ay x l) = QgXg(Ay x l) 


(HI) 


(H2) 


where Ay x 1 is the surface area transmitting heat to the coolant for 
a depth of 1 inch. Prom this equation, 




In like manner. 


t B, 2(A7) = T 2 (Ay) + — — 


(H3) 


(H4) 
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Prom reference 15 (p. 657), the general differential equation for 
heat flow through a slab with a constant heat source is 


-k 


d 2 T 

dx 2 


= Q 


which when integrated results in the equation shown in reference 15 (p. 
658) , which in the notation vised herein is 

T b (A 7 ,x) = - x + T BA (Ay) (H5) 

where T(Ay,x) is the temperature of slab B at some axial point Ay 
from the coolant entrance and some radial point x from the left side 
of the slab. 


Differentiating this equation. 


dT B (Ay,x) 2Qg X T B 2 (Ay) - T B i(Ay) Q B ^B 

= -^r-t tjj + ^~ 


dx 


2k- 


B 


2k 


B 


(H6) 


m 


For the position Xg, dT B (Ay,x)/dx = 0, for which equation (H6) results 

(H7) 


_ ^B + T B,2^) ~ T B,l(Ay) k B 
% Or 


X B = — + 


From equations (H5) and (H7), 

Ob 


T B,max^) ' T B,1^) = " 2k 


B 


[% ®B,2(Ay) - TB,i(Ay) 

L 2 + *B 


[ T B,2^ ' T B,l (Ay) . Vb] 


L 2k B J 

L 2 + *B q b J 


(H8) 


If Xg is at the center of the slab, T-g 2 (Ay) - T-g ^(Ay) = 0, 
and equation (H8) becomes 

,2 


QB 


•'\B,max^^ r ) " ^B , 1 (Ay ) gk-, 


B 




(H9) 
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Prom equations (H3) and (H4) , if wj_ and Wg are metered so that 

T-^Ay) = T 2 (Ay), 


- *b) 

Putting this in equation (H7) and solving for Xg, 

_ (tg/2) + (kg/hg) 

Xb ^.(JL 

t B ^h 2 + hj 

If T-^(Ay) = T2 (Ay) , equations (Hi) and (H2) show that 

Zl - ^ + r C,i% x C 

Wl T k,SpS % k - X A) + r B,iVB 


Finally, 


Zl izif' 8 
h i " W) 


(H10) 


(nil) 


(H12) 


(H13) 


Now, equations similar to equation (Hll) can he set up for x^ and 
X C* the values of h, and Q are assumed constant axially, by 

inspection of equations (Hll), (H12), and (H13), it can be seen that Xg 
would be constant axially if x^ and xq were constant. Since x^ 
and xq depend on the same type of equations as Xg, they would be 
constant. 


_Since the maximum shield temperatures occur at the passage exits 
if Q is assumed constant axially, and since x values are constant, 

Ay in the equations can be replaced by the axial length of the passages 
to find the outlet coolant temperatures from equations (Hi) and (H2) and 
then the shield temperatures from equations (H3) , (H4), and (H8) . 

Volumetric heat sources . - Five percent of the energy generated in 
the core is assumed to be dissipated as heat outside the core, as men- 
tioned previously in appendix E. Then the heat flux or intensity coming 
from the core to the reflectors and shields per unit area of core sur- 
face is 


I 


o 


0.05^x948.1 



(hi 4 ) 
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where 9 is the power of the reactor in megawatts. Because the con- 
tainer is so very thin, the gamma flux can be considered constant, and 
the value of the heat-source strength in it due to the absorption of 
gammas is 


^con “ ^o^^a^con 


(H15) 


Reverting to figure 65, if slab A is considered to be the one next 
to the core container with a coolant passage between, and attenuation 
through the coolant is neglected, then the flux at the entrance to slab 


A will be 



(H16) 


The factor, ratio of radii squared, is the allowance for the geometric 
attenuation. 


The volumetric heat source at this face will be 


(Qa)i A, 1 (lA)r A,i^ a)A 


(HI 7 ) 


Now at any point x in slab A, 


(^A^x ^A^r Aj . 


■(^a) A W 


in which it is assumed that omitting the geometric attenuation 
[ r A,i/( r A,i + x )] 2 compensates for omission of an energy buildup factor 


within the slab: 


Q A At A 


■l 


'A 


(Q a )xA to 


where A is the cross-sectional area of the slab normal to the gamma 
ratiation. Integrating the latter equation after inserting the preced- 
ing expression for (Q^) x results in the following equation for the aver- 
age volumetric heat source: 


Q A - ( Q A^r A ^ 


"^a)A t A 
1 - e 


(H18) 
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In similar manner, but with attenuation being assumed in slab A, formu- 
las for slab B would be as follows : 


(r \ - (t \ l^ a)A V r A,o\ 

(lB)r B,i = (l A } %i e y r B ,i J 


'Wa 


(H19) 


2 

where e ~~ is the attenuation for the material., and ( r A^o/ r B,i^ 

is the geometrical attenuation of the flux at the radius r^. Q . As be- 
fore, omission of (r^ ±/ T ^ Q )^ is assumed to compensate for the lack of 
an energy buildup factor: 


4 " ^B^r-g^^a^B 


%i 


(H20) 


and 


Qb - (%) %1 - 

Equations similar to the foregoing can be set up for all the slabs. 


(H21) 


Convective heat-transfer coefficients . - Equation (E57) for the 
fuel-element heat-transfer coefficients was also used to determine the 
coefficients in the coolant passages between the shields. For the hy- 
draulic diameter, the value for an annulus was used which, for the pres- 
ent case, is twice the space between the shields or twice the coolant- 
passage thickness. As in the case of the fuel elements, the coefficients 
depend on the shield temperatures, which in turn depend on the coeffi- 
cients, so that iteration is required. 


General design method . - The steps in the design are as follows: 

(1) Values of x where the maximum temperatures in each shield 
occur are assumed. 

(2) Coolant-flow-rate ratios in adjacent passages are calculated 
using equations of the type of equation (H12) . 

(3) Since the total coolant-flow rate through the reactor is known, 
wq can be determined from 



and the ratios calculated from step (2). 
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(4) When w-]_ and the coolant-flow ratios are known, Wg, W3, and 
so forth, can be calculated. 

(5) The heat-transfer coefficient hj_ Is calculated using an 
assumed effective wall temperature and w^. 

(6) The heat-transfer-coefficient ratios are calculated from the 
approximate formula (eq. (H13)). 

(7) The coefficients hg, hj, and so forth are calculated using 
steps (5) and (6). 

(8) The coolant outlet temperature is calculated using equation 
(Hi). (This is the same for all passages because of the type of anal- 
ysis made, as brought out previously. ) 

(9) Values of x are calculated for the slabs where maximum tem- 
peratures occur with equations like equation (Hll). 

(10) The procedure is repeated until x in step (9) nearly equals 
values used in step (l) . On second trial, step (6) is changed to cal- 
culating h for each passage using equation (E57) and iterating on 
T w , e ff as in appendix E for the fuel- element coefficients. 

(11) The shield temperatures are calculated using equations like 
equations (H3) , (H4) , and (H8) . 
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APPENDIX I 

REACTOR COMPONENTS AND ASSEMBLY PROCEDURE 
By D. W. Drier 

A discussion of the various parts of the reactor is given in the 
"body of this report. This appendix summarizes the various reactor parts 
and outlines the suggested assembly procedure. 

Components 

Core . - The core consists of the following members : 

(a) Two outer support plates (0SP1 and 2, fig. 5) 

(b) Five inner support plates (iSPl, 2 , 3, 4, and 5, fig. 5) 

(c) 2632 Moderator blocks (figs. 5 and 25) 

(d) 1316 Reflector blocks (figs. 5 and 25) 

(e) 640 Fuel-element tubes (figs. 5 and 25) 

(f) 18 Tie bolts (fig. 5) 

(g) 3872 Support bushings (figs. 5 and 25) 

(h) 3912 Spring spacers (figs. 5 and 25) 

(i) Flow-dividing container (fig. 5) 

( ,j ) 120 l/8-Inch-diameter radial support pins (figs. 5 and 47) 

(k) 48 l/ 4-Inch- diameter radial support pins (fig. 5) 

(l) Circular side reflector (ARl, fig. 5) 

Side reflector . - In addition to the side -reflector piece included 
in the core components, there are three circular side-reflector pieces, 
AR2, AR3, and AR4 (fig. 5). 

Inner shield . - The components of what may be called the inner bio- 
logical shield, some components of which form the thermal shield, too, 
are : 

(a) Pressure vessel (fig. 5) 
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(b) Five frontal shields (FS1, 2, 3, 4, and 5, fig. 5) 

(c) Four radial centering vanes 

(d) Four 1^- inch-diameter bolts 

(e) Four annular flow dividers (AFD1, 2, 3, and 4, fig. 5) 

(f) Six annular thermal shields (AS1, 2, 3, 4, 5, and 6, fig. 5) 
Rear reactor shield . - The components of the rear reactor shield 

are : 

(a) Pressure-vessel head (fig. 5) 

(b) Four rear shields (RSI, 2, 3, and 4, fig. 5) 

(c) Four radial vanes (fig. 49) 

(d) Four clamping bars (fig. 5) 

(e) Annular support (AS, fig. 5) 

(f) Five rear annular shields, RAS1, 2, 3, 4, and 5, fig. 5) 

(g) 72 l/2-Inch-diameter radial support pins (fig. 5) 

(h) 96 3/8-Inch-diameter spacer plugs (fig. 5) 

(i) Flow container and ducts (fig. 5) 

(j) Flow-distributing shell (fig. 5) 

Control rod and actuators . - The components of the control- rod and 
actuator assemblies are: 

(a) Three pairs of shim rods (fig. 7) 

(b) Three pairs of housing tubes for shim rods (fig. 7) 

(c) Three sets of electrically controlled servo actuators used for 
shim rods 

(d) One pair of central regulating rods (fig. 7) 

(e) One pair of housing tubes for central regulating rods (fig. 7) 

(f) One pair of electrically controlled servo actuators for regu- 
lating rods 
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(g) Three single shim-scram rods (fig. 7) 

(h) Three housing tubes for single shim-scram rods (fig. 7) 

(i) Three sets of hydraulically controlled servo actuators for 
single regulating rods 

(j) Two scram rods (fig. 7) 

(k) Two housing tubes for scram rods (fig. 7) 

(Z) Two sets of pneumatically controlled servo actuators for scram 
rods 

(m) Nine split positioning rings (fig. 5) 


Assembly Procedure 

Core . - The procedure for assembling the reactor core follows: 

(1) Insert the bolts through 0SP1. 

(2) Turn tie bolts vertically with bolt heads supporting the 
assembly. 

(3) Slide flow-dividing container and AR1 over 0SP1. 

( 4 ) Insert 24 l/ 4- inch-diameter support pins radially through AR1 
and flow-dividing container and into 0SP1. 

(5) Insert 18 reflector blocks over the 18 tie bolts, fitting the 
circular ends in 0SP1. 

(6) Insert spring spacers and support bushings into 640 front re- 
flector blocks used for fuel elements. 

(7) Insert the 640 front reflector blocks Into 0SP1. 

(8) Insert odd-shaped reflector blocks surrounding control-rod open- 
ings and filling spaces on outside diameter. These blocks are held in 
plane with end pins. 

(9) Place ISP1 over reflector blocks. 

(10) Insert 18 support bushings with spring spacers over tie bolts 
and through ISP1 into reflector blocks . 
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(11) Insert 640 support bushings with spring spacers through ISP1 
into reflector blocks. 

(12 ) Insert 24 l/ 8- inch-diameter support pins radially through AR1 
and flow-dividing container and into ISP1. 

(13) Assemble first layer of moderator blocks over tie bolts and 
support bushings. 

(14) Fill in spaces surrounding control-rod openings and spaces on 
outside diameter with odd-shaped moderator blocks, and insert end pins. 

(15) Place ISP2 over reflector blocks and continue assembly of core 
in a repetitive manner, as described in the previous steps, to 0SP2. 

(16) Place 0SP2 over the rear layer of assembled reflector blocks 
and insert 24 l/4- inch-diameter support pins radially. 

(17) Tighten tie-bolt nuts in place. 

Side reflector and inner shield . - The five frontal shields (FS1, 

2, 3, 4, and 5) and four radial centering vanes are assembled externally 
with the four annular flow dividers acting as spacers. The assembly of 
the shields, reflectors, and core is accomplished by lowering these mem- 
bers into the pressure vessel. The pressure vessel is supported in a 
fixture attached to its central support pin with its axis of symmetry 
vertical. The assembly procedure follows: 

(1) The four radial centering vanes are positioned vertically in an 
assembly fixture with the rear straight edge down. 

(2) The five frontal shields are then slid over the radial vanes 
beginning with FS1 and separated by the four annular flow dividers. 

(3) After FS5 is lowered into place, it is bolted to each of the 
four radial centering vanes. This completes the external assembly of 
the five frontal shields, clamping them positively together. 

( 4 ) The externally assembled frontal shields are lowered into the 
pressure vessel. They are positioned by the centering pin and the four 
pads with the torque pins fitting in corresponding slots in the radial 
vanes. A cross section of this construction is shown in figure 66, 
section F-F of figure 5. 

(5) The annular shields are then lowered into the pressure vessel 
beginning with AS5. The forward end of AS5 is positioned by the four 
radial vanes over which it slides and the rearward end by radial splines 
mating with the pressure vessel. AS4, 3, 2, and 1 are positioned by 
the four radial vanes. 
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(6) Three annular reflectors (AR4, 3, and 2) are then lowered into 
place , also being positioned by the four frontal radial vanes. 

(7) The assembled core, with several control rods already inserted, 
is then lifted vertically by attaching to the threaded tie-bolt ends and 
lowered into position with the extending radial pins sliding in the in- 
ternal grooves in AR2. 

(8) The 640 fuel-element tubes can then be vertically lowered into 
the core, being guided by the spring spacers. 

Rear reactor shield . - The rear-shield assembly procedure follows: 

(1) The four radial vanes and the flow-distributing shell are bolted 
to RS4 with four bolts. 

(2) RS4 Is supported in a fixture with its axis vertical and its 
rear face down. 

(3) RSI Is laid flat with its front face down. 

( 4 ) The flow container with ducts and the annular support are low- 
ered over RSI. 

(5) 24 l/2-Inch-diameter radial support pins are then inserted 
through AS and the flow container into RSI. 

(6) RS2 And RS3 are then lowered into position, and the remaining 
48 l/2-inch-diameter radial support pins are inserted. 

(7) This assembly is then lowered into RS4. 

(8) RAS5, 4, 3, 2, and 1 are then lowered into position, being 
guided by the four radial vanes and each being separated from the pre- 
ceding by 24 3/8-ineh-diameter spacer plugs. 

(9) The entire rear reactor shield assembly is then clamped tight 
by bolting the four clamping bars to the four radial vanes. 

(10) The rear reactor shield assembly is lowered on the main assem- 
bly containing the core. The four radial slots in the protruding annu- 
lar shields guide the rear reactor shield assembly into position by mat- 
ing with its four radial vanes. 

(11) The pressure-vessel head is lowered over the rear assembly and 
rotated to secure positively the clamping bars that position the rear 
assembly. 
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Control rods and, actuators . - The control- rod and actuator assembly 
is accomplished as follows: 

(1) The remaining control rods with their respective housings are 
externally assembled and inserted into the reactor through the pressure- 
vessel head. 

(2) The split positioning rings are inserted around their respective 
housings in the pressure- vessel head. 

(3) The servo actuators are mechanically attached to their respec- 
tive control rods, and the housing flanges for these actuators are se- 
cured to the pressure- vessel head. 

(4) The reactor is now ready to be installed in the fuselage of the 
airplane . 
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APPENDIX J 


REACTOR BIOLOGICAL SHIELD 

In order to calculate the biological-shield weight, it is necessary 
to know the reactor power, reactor-core size, reflector thicknesses, 
thermal-shield thicknesses, pressure-shell dimensions, and the materials 
to be used. These items are assumed to be known. The biological-shield 
problem then reduces to the determination of the gamma- and neutron- 
shield thicknesses that result in a shaped unit shield of minimum weight 
for a prescribed total radiation dose rate at a specified distance from 
the reactor center. Such calculations currently depend upon the use of 
various assumptions and approximations, and at best only an approximate 
shield weight can be obtained at this time. A brief description of the 
method employed is given herein. 

The leakages of the various gamma energy groups and of the neutrons 
out of the reactor are evaluated. The reactor is then replaced by an 
isotropic point source having a source strength equal to this leakage. 

The reactor and shield assembly is divided into n conical shell sec- 
tors (see fig. 54). The radiation that emerges from each sector contrib- 
utes to the prescribed total dose rate at the crew compartment as di- 
rect and single air-scattered radiation. The gamma- and neutron-shield 
thicknesses external to the pressure shell are determined for each sec- 
tor by the use of Lagrange multipliers so that a minimum total shield 
weight will result. 


Symbols 

The following symbols are used in appendix J : 
a distance from reactor center to exterior of reflector, cm 

B buildup factor 

b distance from reactor center to exterior of pressure shell, 

cm 

D dose rate, rem/hr 

d distance from reactor center to crew compartment, cm 

E q average energy of gamma rays, Mev/photon 

H height of reactor core, cm 
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factor for conversion from neutrons/ (cm 2 ) (sec) to rem/hr 
factor for conversion from photons/ (cm 2 ) (sec) to rem/hr 
total leakage from reactor core, neutrons/sec or photons/sec 
Lagrange multiplier 

leakage from reactor core per cm 2 , neutrons/ (cm 2 ) (sec) or 
photons/ (cm 2 ) (sec) 

number of sectors 

reactor power, Mw 

scattered gamma dose rate per unit energy emitted from iso- 
tropic source for various shadow-shield angles, 

(rems/hr)/ (Mev)/sec) 

radius of reactor, cm 

radial variable, cm 

source strength, neutrons/ (cc) (sec) or photons/ (cc) (sec) 

leakage from shield through solid angle dfl, neutrons/sec or 
photons/ sec 

thickness , cm 

actual metal thickness of thermal shields and pressure shell, 
cm 

volume, cc 

void fraction, (b - a - "^ths^/fa ” a ) 
weight, g 

angle of distribution of radiation from shield with normal 
to surface, radians 

angle from which neutrons can emerge from shield and enter 
detector directly (see fig. 67) 

absorption coefficient, cm"^ 

shadow- shield angle, radians 
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p 

density, g/cc 

2 

macroscopic cross section, cm - '*' 

9 

neutron flux, neutrons/ (cm^) (sec) 

t 

sector angle, radians 

A cos k \|r 

cos - cos ^2 

dfl 

solid angle, radians 


Subscripts : 


c 

reactor core 

in 

inside of pressure shell 

k 

sector 

N 

neutron 

out 

outside of pressure shell 

R 

reflector 

S 

shield 

s 

refers to thermal shield modified to account for void 

T 

total for all sectors 

t 

total for any one sector 

ths 

thermal shield 

V 

volume 

r 

gamma 

1,2 

sector angle limits 


Superscripts : 


A 

attenuated through "biological shield 

D 

direct 
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m exponent of neutron distribution from shield 

n exponent of gamma distribution from shield 

s scattered 


Sources and Leakages 

The fast-neutron and gamma leakages from the reactor were evaluated 
and the reactor was then replaced by an isotropic point source giving 
these leakages. Assuming a uniform power generation, the number of fast 
neutrons produced per unit volume becomes 


S V,N " 


7.75(10 16 )^ 

V 


) 


neutrons 
(cc) (sec) 


(Jl) 


Assuming that five 3-Mev photons are emitted per fission, the number of 
photons produced per unit volume becomes 


_ _ 1.55(10 17 )^ photons 

v,T - v ’ (cc) (sec) 


(J2) 


The average energy of the neutrons and gammas was assumed to he 3 Mev. 


The uncollided leakages out of the reactor core were assumed to he 
the same as those out of the surface of a semi-infinite slah. These 
leakages are good approximations along most of a cylindrical surface if 
the radius of the cylinder is sufficiently large. Near the ends these 
leakages exceed those out of a cylinder. 


These leakages are 


and 


l - IM neutrons 
^ 4 ^c ’ (cm 2 ) (sec) 


i = fL-x photons 
^ 4 ^c (cm 2 ) (sec) 


and the total leakages are 

L N = V 2 **" + 2ltRH )j> 


2 . o neutrons 


sec 


(J3) 


(J4> 


(J5) 
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Iy = Z r (2«R 2 + 2«RH) , £ h °g° nS (J6) 


Dose Rates 


The radiation leaking out of the surface of the shield was taken to 
be that coming from the reactor core only. Capture gammas , inelastic 
scattering gammas, and other gamma radiation originating outside the 
core were not considered. Of the radiation leaking out of the shield, 
the amount reaching the detector (or crew compartment) directly and the 
amount reaching the detector after single air scattering were determined. 
These radiation amounts were obtained for a number of sectors into which 
the reactor and shield were divided. 


Direct neutron dose rate . - The neutrons entering a sector of solid 
angle dfl were attenuated through the known thicknesses of the reflector, 
the thermal shields, and the pressure shell. Allowance was made for the 
amount of void within the pressure shell. The neutrons were further 
attenuated through gamma and neutron shields of unknown thicknesses ex- 
ternal to the pressure shell. The neutron leakage out of the k^ 1 sec- 
tor of the shield was thus expressed as 


Sr 


k,N 


= L N 


dft 

4ji 



e 


(J7) 


where the only unknowns are t^y* ^k,N a ^ er "the shielding mate- 

*%^k,R 


rials are selected. The factor e 


accounts for the neutron 


,th 


sector of the reflector. The factor 


attenuation through the 

e accounts for the neutron attenuation through the k^ 1 

sector thermal shields and pressure shell. If the thermal shields and 
the pressure shell are made of different materials, this factor must be 
replaced by the product of two exponentials, one for attenuation through 
the thermal shields and the other through the pressure shell. The void 
within the pressure shell is accounted for as follows: 

tk,s = ^thsA 1 " v k) 311(1 2 s = 2 ths^ " v k) > where v k ls the void 
fraction and -^is is the actual metal thickness of the thermal 

- ( 2 r t k,r +2 N t k ; N^ 

shields and pressure shell. The factor e accounts 

for the neutron attenuation through the biological shield exterior to 
the pressure shell. 


The neutron leakage out of the shield surface is assumed to have 
an angular distribution proportional to cos m 0, where 6 is the angle 
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the emergent neutron makes with the normal to the surface. If 
neutrons leaking out of the surface have this distribution, there will 
be |(m + l)/2«J cos m 0^ neutrons per unit solid angle in direction 

0. From figure 67 it can be seen that only those neutrons that emerge 
at an angle 0' will enter the detector directly. For large separation 
distances as encountered here, 0' is approximately equal to \|r. 


The total flux into a detector at a distance d from the reactor 

■f - Vi 

center, for the k^ sector, was then found by integration over the 
angular width of the sector (\|r^ to tyg) • Introduction of a factor K 

for converting neutrons/ (cm^) (sec) to dose-rate units rem/hr finally 
produced the equation for the direct neutron dose rate into the crew 

compartment from the k bn reactor sector: 


1 % / _2 R^k,R\ / “ z ths^k,ths ] 

wl e A e ) 

x ^ ^ _ cosm+1 ^ (js) 


Scattered neutron dose rate . - Consider an isotropic source emit- 
ting Ljj neutrons per second, and consider only those neutrons emitted 

into the solid angle dfl contained in the conical shell bounded by 
angles and tyg* The single scattered flux into a detector at dis- 

tance d from the source due to the neutrons in dft is 


q> = 



8jtd 



(J9) 


g 

where Z is the total macroscopic scattering cross section for air 
evaluated at the appropriate altitude (see ref. 36). 


Now s k ^ N represents the neutrons entering the solid angle dft. 

To facilitate scattering calculations, assume that these neutrons are 
emitted radially from the shield surface (instead of with a cos m 0 
distribution) . The corresponding isotropic source becomes 

= s k,N ( J1 °) 
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Insertion of equations (J7) and (J10) into (J9) and introduction of the 
conversion factor K finally yields the scattered neutron dose rate at 
the crew compartment as 



Direct gamma dose rate . - The direct gamma dose rate at the detec- 
tor at a distance d from the reactor center is calculated in a manner 
analogous to that used for the direct neutron dose rate. In the case of 
gammas, however, the gamma leakage replaces the neutron leakage 

Ljj, a new conversion factor K' to convert photons/ (cm^) (sec) to rems 

per hour is used, reflector and shield buildup factors are incorporated, 
\i replaces 2, and a new distribution at the shield surface (cos n 6) is 
assumed. The direct gamma dose rate can then be written as 



L r 

4rtK'd 2 



-n 


ths kjths 


)i 






] 


x 


/ x / H+l 

(ErEs) ( c °s if 


n+l . 
cos t 2 ) 


(J12) 


Scattered gamma dose rate . - The total isotropic source, 
as for neutrons, is 


obtained 


T A 4k 

L r - s k,r an 


= L r 




(J13) 


Reference 37 presents data for evaluating p(^), the scattered gamma dose 
rate per unit energy emitted from an isotropic source for various shadow- 
shield angles £. An approximate analytical expression for p(|) for 
3-Mev photons was found to be 

/ x Paltltude 1500 2 [5. 898(l0 -8 ) | 2 -9. 4l(l0" 4 )] , x 

PU; _ 0.0765 d 10 15 S W } 
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The scattered dose rate from any sector bounded by angles and 

will be p(i|f-j_) - p(i 2 ) (rems/hr)/(Mev/sec) . Hence, the dose rate from 
the k^h sector is 


E k,r ‘ [ p( *i> ’ 


(J15) 


where E Iy is the energy emission from the source. Substitution for 
Lp leads to the equation 


,S .. k-^M/'/ p thB^,ths\|‘y(Vk,r +p N%,N)] BRBgp ^ i _ (J16) 


E k,r - *<*■ le 


Total dose rate . - The total radiation dose rate reaching the crew 
compartment is the sum of the individual dose rates for all sectors, or 


n 


% - £ <i£,n + E k,N + <r + D k, r ^ 


(J17) 


k=l 


Weight 

An element of volume of one of the conical-shell sectors of the 
shield is 


Therefore 


dV = 2jtr sin r di|r dr 


2 

dW = p2rtr sin i|r di|f dr 


(J18) 


where p is the density of the shield material considered. The weight 
of the sector can now he found hy double integration between the 

appropriate limits : 


W = 2itp f 2 f 


* 


2 ? 

r dr sin \| r d+ 


i % 

/ r 2 2 

r dr 


A cos i k (r^ - r^) 


(J19) 


Declassified 1-24-1963 


UJH^XUUX 


Declassified 1-24-1963 


154 


D-Ciunjux 


where rj_ and r 2 are the radial limits for the particular piece of 
shield under consideration, and A cos i|r k represents the difference in 
cos evaluated at the angular limits of the sector under consideration. 


Gamma-shield weight inside pressure shell and including pressure- 
shell weight . - In this case, the lower radial limit is the distance 
from the reactor center to the outer edge of the reflector a^. This 
distance increased by the thermal- shield thickness adjusted to account 
for the void is used for the upper radial limit. Measurements are made 
along the centerline of the sector under consideration. To account for 
the void in this region, the density y. is replaced by y(l - v^), 

and the thermal-shield thickness t^ s is replaced by t^-t^g/ (l - v k ) . 
Hence, 


W k,r, in 


2 * 

3 


py(l - v k )A cos i|/ k 



3 

a k 


(J20) 


Gamma shield exterior to pressure shell . - The lower limit in this 
case is the distance from reactor center to the outer edge of the pres- 
sure shell denoted "by b^, and the upper limit is this distance increased 
by the gamma-shield thickness yet to be determined. Again, measurements 
are made along the centerline of the sector. The result is 


w k,r,out - A cos + t k<r ) 3 - b 3 ] (J21) 

Weight of neutron shield . - The neutron shield is located exterior 
to the outer gamma shield. The inner radial limit is therefore equal to 
the outer radial limit for W k ,^., ou t> the outer radial limit is this 
value increased by the yet to be determined neutron-shield thickness. 

The neutron- shield weight is given by 


W k,N = A cos ♦kjfak + + t k,N) 3 “ ( b k + ^/r) 3 ] ( J22 ) 

Total shield weight . - The total shield weight is the sum of the 
over-all sectors of the three weights just described, or 

n 

W T “ ^ ( w k,r,in + w k,r,out + w k,N^ (J23) 

k=l 
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Gamma- and Neutron-Shield Thicknesses 

For a given reactor power and size, the number and thicknesses of 
the thermal shields and the thickness and size of the pressure shell are 
known. Hence, Wy j_ n also remains fixed. The gamma- and neutron-shield 
thicknesses for those parts of the shield exterior to the pressure shell 
must be determined. Use of Lagrange multipliers permits the calculation 
of these thicknesses that will give a minimum-weight shield for a pre- 
scribed total dose rate at a specified location. The known necessary 
conditions are 


3d, 


T 


cWj _ 

^k,y ■*- ,a ^k,y 


and 


3w 


T 


SD, 


T 


^tk,N + La dt k,N 


= 0 


(J24) 


(J25) 


In conjunction with the subsidiary condition imposed by the prescribed 
total dose rate Drj\, these equations serve to determine the gamma- and 

neutron-shield thicknesses and the constant of proportionality La, 


From equations (J24) and (J25), the value of La can be written 


as 


L a — _ 


SDip 

fek,y 

C5W rji 


3d, 


T 


dt k,N 


(J26) 


"5t 


k ,r 


*5t 




as shown in reference 34. Differentiation of equations (J8), (Jll), 
(J12), and (J16) with respect to tj^y* 


t^ but assuming that the 


buildup factors remain constant , results in the following: 
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^Vr 

^k/r 


= -v£. 


bn. 


k,r 


2JL. = 


^k,r 

3 Vr 


N 


“V^N 


■V°k,r 

“V4,r 


^°£,n 

otv TlT 


k,N 


¥k,N 


" _2 N D k,N 

^ * ~ M N E 5[>r 


If designated the total dose rate for any one sector, then 


SDiJI 

*5t 


Bd 4 


k ,r dt k,r 


" "¥( D k,N + D k ; N) ' ^r( D k,r + D k/r) 


and 


= SD t 

dt k,N ^ t k ; N 


- “%( D k,N + ^n) " + D k,r) 
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Also, if designates the total weight for any one sector, then 

dw T 2w t 


dt, dt, 
k,r k ,r 


= {2«fy(t k , r + 

+ 2*Pg[(tfc, N + ‘k,r + \) Z - (**,>■ + \) 2 ]} A cos (J30a) 


and 


cWip c)W+ 


at k,N ^k,N 

= 2jt PN( t k,N + \,r + b k) 2 A cos \ (J30b) 

Therefore, from equations like (J26) for each individual sector, 

* D k,iir) + ^(^,r + D k,r) 

^fyCt^+bj^) 2 + Snpjj^^ + t^ + b^ 2 - (t k , r + b k ) 2 ]}A cos ^ 

(J31) 


L a — 


and 


La = 


%( D k,N + D k,w) + + P k,r) 

+ \, r + b k ) 2 A cos *k 


(J32) 


If equations (J3l) and (J32) are equated for any sector, the only 
unknowns are the gamma- and neutron-shield thicknesses y. ^k,N’ 

For a series of assumed values for t^ corresponding values of t^. jj 
can then be found, and a relation between t^. ^ and t^. ^ can be 
determined. ^ ’ 

From either equation (J3l) or (J32), values of La can be found 
for the corresponding pairs of values of t^y an( ^ ^k N P rev iously 

determined. A relation between La and either t^ ^ or t^. ^ can 
then be established for each sector. 


By use of these values of La and the corresponding values of 
t-, _ and t v in-, equations (J3l) and (J32) can be solved simultaneously 

.for (D? - t + Df J\ and (DT + D? | for each sector. The sum of these 

\ k,N k,N/ \ k,r k,r/ 


SECRET 

Declassified 


Declassified 1-24-1963 


158 


Bj&ujtusr 


quantities then yields for each sector, and the sum of the 

gives the total dose rate Drp. When such calculations are completed 
for a series of values of La, a relation between La and Up can be 

established. From this relation, the value of La corresponding to 
the prescribed allowable total dose rate at the crew compartment can be 
found. The values of t^ y. and tj^jj corresponding to this La can 
then be found, and the shield weights can then be calculated for each 
sector. Finally, the sum of the individual sector shield weights gives 
the desired total weight Wrp. 

Tables XV and XVI give some material properties and constants used 
in the shielding calculations . The values of the final shield thick- 
nesses and weights are given in tables XI and XII. 
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APPENDIX K 


ADDITIONAL SHIELD WEIGHTS TO COMPENSATE FOR DUCTS 

The previous discussion of the "biological-shield calculations did 
not account for the presence of ducts in the shield. Allowances must 
"be made for additional shielding in order to compensate for radiation 
streaming through the ducts and for the voids in the shielding caused 
by the presence of ducts. A simplified method for incorporating these 
effects is described herein. 

Only fast-neutron shielding must be considered for streaming and 
void effects; such neutron shielding is generally believed sufficient 
to shield against gamma rays, too. A method of calculating these addi- 
tional shield weights will be described. The additional shield weights 
will be determined on the assumption that hemispherical bumps will be 
added to the previously calculated biological shield. These bumps are 
to be tanks of decalin, which will also provide additional chemical fuel 
for emergency use. 


Symbols 

The following symbols are used in appendix K: 
d duct diameter, cm 

l length of duct leg, cm 

a reflection coefficient (assumed =1) 

0 angle variation from elbow corner to elbow corner of duct-leg 

exit, radians 

cp^ fast-neutron flux at duct-leg inlet, neutrons/ (cm^) (sec) 
cp Q fast-neutron flux leaving duct-leg exit, neutrons/(cm )(sec) 
Subscript : 
m duct leg 

Calculation Method 

A method was derived that would yield the optimum lengths of ducts 
and an accompanying minimum additional shield weight. However, the 
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geometry of the reactor discussed herein prevented the use of optimum- 
length ducts. As a consequence , the method of optimization could not 
he used, and calculations were made for the particular duct configura- 
tion of the present reactor. Using the hiological-shield procedure 
(appendix j), it is possible to determine the flux at the inlet of each 
duct leg; these calculated fluxes are the fluxes that would exist at 
the specified locations if there were no ducts (e.g., for a solid shield) 
An additional flux due to streaming through each duct must also he con- 
sidered. The streaming fluxes were determined by use of the equation on 
page 278 of reference 38: 

9o _ 1 1 (ii\ 2 _o_ 1 feV _a_ 1 K) 2 a 

CP! 8 yi x J 8 \ l 2 J sec ^ 8 \ l 5 ) sec 0 2 ' ’ ' 8 \l m J sec 0 m 

The values of d^, 2 m , <3^/2^ and 0 m for the reactor considered herein 
are as follows : 


Leg 

Duct 

diameter, 

4, 

cm 

Duct 

length, 

l, 

cm 

d/2 

V 

deg 

1 

1.27 

12.70 

l/lO 

90 

2 

20.32 

30.48 

2/3 

60 

3 

27.94 

50.80 

11/20 

90 

4 

27.94 

121.92 

11/48 

45 

5 

27.94 

137.16 

11/54 

— 


The value of cfn used in the determination of the streaming fluxes of 
each leg was taken as the sum of the solid-shield flux and the streaming 
flux of the preceding leg. The total flux for each leg was the sum of 
this calculated streaming flux plus the so-called direct or solid-shield 
flux at the leg outlet. 

The sum of the two fluxes (that for the solid shield plus that for 
the streaming) was then attenuated through the remaining hiological- 
shield lengths. The thickness of the additional shielding required to 
attenuate this flux to the allowable flux was then determinable; the 
allowable flux was found from the biological-shield calculations of the 
neutron dose rates in the appropriate sector of the shield. The weight 
of the additional shield required because of streaming was then deter- 
mined under the assumption that it would be added as a hemispherical bump 
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These bumps would then be placed on the exterior of the biological 
shield in direct line with the respective duct under consideration. 

Additional shielding due to the voids in the biological shield 
caused by the presence of ducts is also required. The volume of each 
duct was calculated, and its weight was determined. This weight was 
then projected to the outside of the biological shield, and the addi- 
tional shielding required because of the presence of ducts was thereby 
determined. 
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TABLE I. - PRINCIPAL ASSUMPTIONS, DESIGN, AND 
OPERATING CONDITIONS FOR TYPICAL AIRPLANE 


Airplane : 


Gross weight, lb 

400,000 

Flight Mach number 

0.72 

Altitude, ft 

30,000 

Wing : 


Area, sq ft 

6250 

Span, ft 

250 

Sweep, leading edge, deg 

18.5 

Design lift coefficient 

0.28 

Tail area, sq ft 

2500 

Fuselage : 


Length, ft 

175 

Diameter (nominal), ft 

15 

Lift-drag ratio (including nacelles) 

19.7 

Engines : 


Number 

8 

Airflow, each engine at flight conditions, lb/sec 

57.0 

Nacelle diameter, ft 

4.0 

Compressor pressure ratio 

11.0 

Turbine-inlet temperature, °R 

1800 

Propeller shaft horsepower, each engine at flight 


conditions 

3410 

Ratio of turbine-inlet to compressor-exit pressure 

0.85 

Reactor : 


Power to coolant, Mw 

98.5 

Helium: 


Flow, lb/sec 

74 

Inlet pressure, lb/sq ft 

180,000 

Inlet temperature, °R 

1250 

Outlet pressure, lb/sq ft 

172,800 

Outlet temperature, °R 

2250 

Heat exchanger: 


Number 

8 

Helium: 


Inlet pressure, lb/sq ft 

170,200 

Inlet temperature, °R 

2250 

Inlet Mach number 

0.024 

Outlet pressure, lb/sq ft 

168,500 

Outlet temperature , °R 

1200 

Air : 


Inlet pressure, lb/sq ft 

9070 

Inlet temperature, °R 

985 

Inlet Mach number 

0.063 

Outlet pressure, lb/sq ft 

8330 

Outlet temperature, °R 

1200 

Helium pump: 


Pressure ratio 

1.085 

Efficiency 

0.80 

Temperature rise, °R 

50 

Horsepower per engine 

872 
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TABLE II. - WEIGHTS OF TYPICAL AIRPLANE COMPONENTS 


Aircraft structure, lb 


108,000 

Fuselage, lb 

28,000 


Wing, lb 

41,000 


Tail, lb 

9,000 


Landing gear, lb 

30,000 


Equipment, lb 


22,000 

Engines, pumps, and nacelles, lb 


58,000 

Heat exchangers and lines, lb 


25,000 

Reactor and shield, lb 


120,000 

Core, lb ) 

4,500 


Reflector, lb/ 



Thermal shield, lb) 

36,000 


Pressure shell, lbj 
Gamma shield, lb 



30,500 


Neutron shield: 



Lithium hydride, lb 

18,500 


Decalin (chemical fuel), lb 

30,500 


Decalin for takeoff and landing, lb 


4,000 

Payload, lb 


63,000 

Gross weight, lb 


400,000 


TABLE III. - DESIGN CONDITIONS HELD CONSTANT 


Efficiencies : 


Compressor polytropic 

0.88 

Turbine poly tropic 

.85 

Reduction gear (propeller and helium pump) 

.95 

Propeller 

.80 

Inlet total- pressure ratio 

.95 

Exhaust-nozzle velocity coefficient 

.96 

Compressor rotor inlet relative Mach number at tip 

.95 

Turbine rotor inlet Mach number at hub 

.60 

Turbine rotor exit axial Mach number at hub 

.47 

Hub -tip radius ratio: 


Compressor inlet 

.50 

Turbine exit 

.60 

Turbine blade centrifugal stress at rotor exit hub, psi 

30,000 

Turbine blade material density, Ib/cu ft 

500 

Blade taper factor, ^ 

.70 

Compressor equivalent tip speed, ft/sec 

1000 

Ratio of turbine-exit to compressor-inlet pressure 

1.038 
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TABLE IV. - FINAL REACTOR-CORE VOLUMES 


Power = 100 Mw 

Number of fuel elements = 640 
Core length = 24 in. 

Side reflector thickness = 3.75 in. 

End reflector thickness = 4 in. 

Uranium = 140.5 lb of 92.5 percent an a 7.5 percent U®®® 

Fuel element "meat": UC-Mo mixture 

Core radius = 17.36 in. 


Item 

Volume, 
cu in. 

Volume, 

percent 

UC in fuel elements meat 

312 


1.4 

Mo in fuel elements tubes 

320 



Mo in fuel elements spacers 

21.6 


1.72 

Mo in fuel elements sleeves 

48 



Mo in moderator-block support plates 

265 

> 

1.53 

Mo in core container 

81.8 


Structure 

Mo in control rods 

635.8 


2.80 

Inconel in through -bolts (70 percent Ni) 

47.8 


.21 

He in core (void) 

3094 


13.61 




Void 

BeO in moderator blocks 

17,598 


77.30 

BeO in control rods (rods out) 

326 


1.43 

Total 

22,750 

100.00 


Flow area, 118 sq in. 
Frontal area, 947 sq in. 


Flow area 
Frontal area 


0.125 
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TABLE V. - REACTOR DATA FOR TWO OPERATING CONDITIONS 



Loaded core for ini- 
tial startup 

Partially depleted core 
37.4 days at 100 Mw 

Wu, lb 

140. 

5 

130 

.0 

k eff 

(cold clean) 

1.18 



k eff 

(hot clean) 

1.12 

1. 

11 

(hot dirty) 



1. 

08 

Group 

Energy 

Neutron 

Fissions, 

Neutron 

Fissions, 



flux. 

percent 

flux. 

percent 



neutrons 


neutrons 




(sec) (cm 2 ) 


(sec) (cm 2 ) 


1 

10 Mev — 2 *46 Mev 

1.57 x 10^ 


1,48x10^*3 


2 

X-V X I'w* v U ft J. v 1 V 

2 46 Mev — 0 9 OR Mev 

1 6 40 


16 1 R 


3 

tj • * v l iw v v « w vu riw v 

0*908 Mev - 25 Kev 

65 86 


J-U * -LCJ 

65 52 


4 

V l v w vj X lv^ V Li vJ Ilw Y 

25 Kev - 50.6 ev 

62.19 

31.13 

63.58 

29.58 

5 

50.6 ev - 0.411 ev 

18.02 

34.50 

19.31 

34.28 

6 

0.102 ev 

3.22 

34.37 

3.66 

36.13 
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TABLE VI. - WEIGHTS AND DENSITIES 
OF CONSTITUENTS IN CORE 


Constituent 

Weight, 

lb 

Weight, 

percent 

UC 

153.1 

6.1 

He 

.5 

Negligible 

Mo 

506 

20.1 

BeO 

1845 

73.2 

Ni 

15.4 

.6 

Weight 
of core 

2520 

100.0 


Constituent 

% 

atoms/ cc 
core 

He 

8.6*10 19 

Mo 

387.5 

Be 

5401 

0 

5401 

U235 

40.39 

y 238 

3.28 

C 

43.67 

Ni 

19.19 
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TABLE VII. - VOLUMES AVAILABLE AND RE- 
QUIRED IN FUEL ELEMENTS FOR FUEL 
[Uniform radial power distribution.] 


Reactor 

sector 

Volume 
available, 
cu in. 

Volume 
required 
for UOg, 
cu in. 

Volume 
required 
for UC, 
cu in. 

1 

0 

a 2 .47 

a 1.85 

2 

13.09 

7.48 

5.58 

3 

20.37 

12.82 

9.56 

4 

40.72 

18.82 

14.04 

5 

34.18 

25.89 

19.32 

6 

40.72 

34.69 

25.88 

7 

61.8 

46.16 

34.44 

8 

87.3 

61.79 

46.09 

9 

77.1 

a 83.57 

62.34 

10 

90.2 

a 108 .87 

81.22 


a Available volume insufficient. 


TABLE VIII. - PROPOSED CONTROL-ROD 
SYSTEM COMPONENTS AND WORTH 


Control rods 

Worth, 

8k/rod 

Actuating 

mechanism 

Number 

Type 


2 

3 

(Pairs ) 
3 
2 

Regulating 

Shim 

Shim-scram 

Scram 

0.005 

.04 

.04 

.02 

Linear re- 
luctance 
motor 

Synchronous 

reluctance 

motor 

Hydraulic or 
pneumatic 

Pneumatic 

Total 

10 


Total 

0.285 + 
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TABLE IX. - REACTOR STRESSES 
fin psi. Loading, 3 g's.] 


Pressure vessel 

Core 

Rear shield 
assembly 

Control rod 
housing tubes 

Front cen- 
tering pin 

1. Wall (1.75 in.) 

Hoop membrane 22,400 

Average longitudinal 11,200 

1. Supporting bushings 
Longitudinal 

bending Negligible 
Bearing 57 
Shear 58 
Hoop bending 3780 

1. Fin posi- 
tioning 
Shear 5133 

Max. bending 3, 
1970 

Max . def lec - 
tion, in. 

0.015 

Shear 2100 

Bending 2700 

2. Flange^ (without gussets) 
Longitudinal hub 31,200 

Radial flange 17,220 

Tangential flange 10,700 

Maximum equivalent 26,435 

2. Spring spacers 

Tensile 176 

2 . Fin bolt 
Tensile 

13,410 



3. Bolt 

Tensile 19,400 

3. Pins for fuel-element 
support plates c (angu- 
lar acceleration = 

6 radians/sec^) 

Max. shear 8900 

Max . bearing 3850 

Max . bending 5030 





a In most unfavorable position. 

^Method of computation given in reference 39. 

c It is assumed that support plates are at uniform temperature throughout and no thermal stresses 
are generated. 
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TABLE X. - SHIELD WEIGHTS FOR VARIOUS 

COMBINATIONS OF MATERIALS 

[lOO Mw, 100 ft, 0.025 rem/hr, 0.72 Mach number, 
30,000 ft.] 


Inner 

gamma 

shield 

material 

Outer 
gamma 
shield 
mater ial 

Six 15° 
neutron 
shield 
material 

1° 

Four 22 -g 

neutron 

shield 

material 

Weight, 

lb 

« 

U 

*U 

LiH 

LiH 

73,000 

F77 

**F77 

Decalin 

I 

89,000 

W 

W 

Decalin 

\ 

91,000 

Fe 

W 

Decalin 

LiH 

96,000 



F77 

| 


96,000 



Fb 

1 

\ 

97,000 



Fe 

LiH 

LiH 

98,000 



1 

Decalin 

LiH 

103,000 






(Design) 



T 

LiH 

Decalin 

107,000 



Pb 

h 2 o 

h 2 o 

111,000 



Fe 

Decalin 

Decalin 

112,000 



Fe 

H20 

H 2 O 

121,000 


Depleted. 

**Fansteel 77 metal (89 percent W) . 
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TABLE XI. - EXTERIOR SHIELDING IN EACH SECTOR 

[lOO Mw, 100 ft, 0.025 rem/hr, 0.72 
Mach number, 30,000 ft *3 
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Sector 

Angle, 

deg 

Gamma-shield 

thickness 

(Fe), 

in. 

Neutron-shield 

thickness, 

in. 

Decalin 

LiH 

1 

15 

4.8 

55.8 

- 

2 



5.7 

35.8 

— 

3 



7.0 

36.7 

— 

4 



10.7 

38.4 

— 

5 



7.2 

40.3 

— 

6 



6.0 

39.7 

— 

7 

22 

:.5 

5.0 

_ 

36.8 

8 

1 

1 

1.3 

— 

37.8 

9 



.7 

— 

38.2 

10 



0 

— 

37.1 


TABLE XII. - SHIELD WEIGHTS FOR VARIOUS REACTOR POWERS 
AND SEPARATION DISTANCES 

[0.025 rem/hr, 0.72 Mach number, 30, 000 -ft altitude, 
Fe -Fe-decalin-LiH 


Reactor 

power, 

Mw 

Distance, 

ft 

lb 

V 

Wr,o, 

lb 

lb 

lb 

50 

50 

35,000 

28,000 

25,000 

88,000 

21,000 


100 

31,000 

28,000 

18,000 

77,000 

18,000 


150 

29,000 

28,000 

15,000 

72,000 

16,000 

100 

50 

45,000 

35,000 

38,000 

118,000 

26,000 


100 

40,000 

35,000 

28,000 

103,000 

23,000 

(Design) 


150 

37,000 

35,000 

24,000 

96,000 

21,000 

150 

50 

50,000 

38,000 

45,000 

133,000 

30,000 


100 

45,000 

38,000 

36,000 

119,000 

26,000 


150 

42,000 

38,000 

30,000 

110,000 

24,000 

200 

50 

55,000 

41,000 

52,000 

148,000 

32,000 


100 

49,000 

41,000 

42,000 

132,000 

29-, 000 


150 

46,000 

41,000 

38,000 

125,000 

27,000 
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TABLE XIII. - FACTORS REQUIRED IN SOLUTION OF FUEL-ELEMENT AND MODERATOR TEMPERATURES 



882 -H 


Declassified 1-24-1963 


E-288 


Declassified 1-24-1963 




175 


TABLE XIV. - SOURCES OF PROPERTIES USED IN STRESS CALCULATIONS 


Material 

E 

a 

V 

k 

Mo 

Ref. 16 

(pp. 195-198) 
Ref. 11 

(fig., p. 17) 

Ref. 16 
(p. 193) 

Ref. 11 

Ref. 16 
(p. 193) 
Ref. 40 

BeO 

Ref. 15 

(table A-3) 
Ref. 41 

(table IX) 

Ref. 15 

(table A-3) 
Ref. 41 
(p. 50) 

Ref. 15 

(table A-3) 

Ref. 41 
(p. 51) 


TABLE XV. - MATERIAL PROPERTIES 
USED IN SHIELD ANALYSIS 


Material 

Density, 

/' 

g/cc 

I ; 1 
cm L 

^1 
cm ^ 

Core BeO 

2.20 

0.1111 

0.0714 

Reflector BeO 

2.20 

.1111 

.0714 

Thermal- and gamma - 
shield Fe 

7.8 

.1667 

.2703 

Pressure-shell Fe 

7.8 

.1667 

.2703 

Neutron-shield LiH 

.79 

.1160 

.0280 

Neutron-shield decalin 

.895 

.1019 

.0349 
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TABLE XVI. - CONSTANTS USED IN SHIELD ANALYSIS 

[Reactor power, 100 Mwj reactor-core height, H, 60.96 
cmj reactor-core radius, R, 44.02 cmj reflector 
thickness, tp, 10.16 cmj distance from reactor center 
to crew compartment, d, 3048 cm.] 


Sector 

Range of 
sector, 
deg 

a k> 

cm 

hk, 

cm 

Je 

fy^ths* 

cm 

b k " a k “ hkjths 
v k t»k - a k 

1 

0-15 

40.44 

108.80 

55.36 

0.190 

2 

15-30 

43.30 

109.80 

50.54 

.240 

3 

30-45 

50.34 

106.70 

42.00 

.255 

4 

45-60 

68.10 

99.16 

23.78 

.234 

5 

60-75 

61.60 

88.62 

21.74 

.195 

6 

75-90 

57.38 

82.56 

20.30 

.194 

7 

90-112^ 

57.80 

83.24 

20.54 

.193 

8 

112'|-135 

68.00 

100.50 

27.40 

.157 

9 

135-157| 

47.80 

110.00 

27.64 

.556 

10 

157^-180 

40.92 

109.00 

42.78 

.372 
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Figure 2. - Picture of airplane model showing location and arrangement of equipment. 
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Figure 3. - Layout of engine, heat exchanger, pump, and headers. 
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Figure 4. - Schematic cross section, of inner pipe of the 
concentric helium lines. 
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Figure 5, - Plan view of reactor showing control and mechanism (section A-A of fig. 52). 
he obtained by using the request card bound in the back of the report.) 


(A large copy of this figure may 
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Figure 6. - Enlarged plan view of reactor (section A-A of fig. 52). 
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Figure 7. - Section C-C of figure 
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Figure 8. - Closeup plan view of reactor pressure shell in airplane model 
with shield indicated. 
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Figure 9. - Details of mounting reactor and shield in airplane fuselage. 
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Figure 10. - Effect of design flight Mach number on airplane component weights. Alti- 
tude, 30,000 feet. 
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- Effect of design flight altitude on airplane component 
Flight Mach number, 0.72. 
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Figure 12. - Effect of design landing speed on air- 
plane component weights. Flight Mach number, 
0.72; altitude, 30,000 feet. 
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Figure 13. - Effect of design gross weight on airplane component 
weights. Flight Mach number, 0.72; altitude, 30,000 feet. 
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Figure 14. - Effect of aspect ratio on airplane component weights. Flight 
Mach number, 0.72; altitude, 30,000 feet. 
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Net jet thrust per unit airflow, F/w a , lb/(lb/sec) 
Airflow per unit turbine frontal area, w a /A T , 
(lb/sec )/sq. ft 
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Figure 16. - Variation in F/w a , w a /Ay, and HP/w a 
with compressor pressure ratio. Flight Mach num- 
ber, 0.72; altitude, 30,000 feet; turbine-inlet 
temperature , 1800° R . 
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Figure 17. - Variation in F/w a , w a /Arp, and HP/^ 


w. 


a 


with turbine-inlet temperature. Flight Mach num- 
ber, 0.72; altitude, 30,000 feet; compressor pres- 
sure ratio, 11. 
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Figure 18. Effect of turbine-inlet temperature on component 
weights. Flight Mach number, 0.72; altitude, 30,000 feet; 
compressor pressure ratio, 11. 
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Figure 18. - Concluded. Effect of turbine -inlet temperature on 
component weights. Flight Mach number, 0.72} altitude, 30,000 
feet; compressor pressure ratio, 11. 
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Figure 19. - Effect of compressor pressure ratio on component weights. 
Flight Mach number, 0.72} altitude, 30,000 feet} turbine- inlet tem- 
perature, 1800° R. 
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Figure 19. - Concluded. Effect of compressor pressure ratio on component 
weights. Flight Mach number, 0.72; altitude, 30,000 feet; turbine- 
inlet temperature, 1800° R. 
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Figure 20. - Variation of weight (shield + ducts + cooling) with reactor power. 
Flight Mach number, 0.72; altitude, 30,000 feet; dose rate, 0.025 rem per hour. 
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Figure 21. - Stress data employed in heat -exchanger studies. 
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and return line 



Figure 22. - Schematic sketch of end view of header arrangement for heat exchanger. 
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Figure 23. - Effect of compressor total -pres sure ratio (inlet pressure) on 
heat -exchanger dimensions and weight for reference airflow rate of 100 
pounds per second. 
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Figure 24. - Proposed design of pump bearings and seals. 
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Figure 25. - Section D-D of figure 5, through support bushing for 
moderator block. 
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Figure 26. - Temperatures of two-coolant-passage fuel element 
and of helium in axial direction for 100-megawatt reactor 
(for "chopped" 2/3-cosine axial distribution of heat). 

( w h = 0-0459 lh/secj w an = 0.0684 lh/sec; ^passages * 2.6$ 
of initial pressure . ) 
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Figure 27. - Stress -rupture strength of molybdenum + 0.45 percent titanium 
for 1000-hour rupture life (extrapolated data). 
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Figure 28. - Allowable and actual fuel-element stresses at outer 
radius (for "chopped" 2/3-cosine axial distribution of heat). 
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Figure 29. - Temperatures of two- coolant -pas sage fuel element and 
of helium in axial direction for 100-megawatt reactor (for axial 
distribution of heat resulting in satisfactory fuel-element 
stresses). (w h = 0.0459 lb/sec; w an = 0.0684 lb/sec; APp assages » 

2.6^ of initial pressure.) 
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Figure 30. - Variations of volumetric heat-source distribution 
in fuel element. 
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Figure 31. - Allowable and actual fuel-element stresses at outer 
radius (for conditions resulting in satisfactory fuel-element 
stresses) . 
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Figure 32. - Effect of reactor diameter (core diam. plus twice reflector 
thickness) on amount of required for constant effective multipli- 

cation factor. 
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Figure 33. - Effective multiplication factor for reactor for various operating 
conditions . 
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Figure 34. - Variation of flux radially along reactor center- 
line for six groups. (Startup, 140.5 lb of uranium.) 
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- Variation of flux axially along reactor centerline for six groups. 
140.5 lb of uranium.) 
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Figure 36. - Required and available volume distribution of uranium carbide 
in reactor core. (Shaded regions indicate that required volume exceeds 
available volume . ) 
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Figure 37. - Temperatures of moderator block and helium in axial 
direction for 100-megawatt reactor (for axial distribution of 
heat resulting in satisfactory fuel-element stresses), (w^ = 
0.0459 lb/ sec; w an = 0.0684 lb/sec; ^passages * 
initial pressure.) 
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Figure 39. - Variation of allowable and actual moderator ~block 
stresses at outer radius with axial length of block (for con- 
ditions resulting in satisfactory fuel-element stresses). 
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Figure 40. - Side -reflector temperatures and flow rates in passages between reflector 
shells . 
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Figure 41. - Thermal- shield temperatures and flow rates through passages. 
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Figure 42. - Block diagram showing power or neutron-flux variations that affect 
reactivity. 
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Figure 43. - Xenon buildup in response to step decrease in neutron flux from equilib- 
rium xenon. 
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Figure 44. - Reactivity variations with time of operation at full power. 
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Figure 45. - Section E-E of figure 5 
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Figure 46. - Fuel tube in reactor core. 
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Figure 47. - Inner moderator block and fuel-element support plate; five required. 
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Figure 48. - Fuel tube and positioning "block. 
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Figure 49. - Hear support shield ’RS 4 . 
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Figure 50. - First front shield plate FS1. 
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Figure 51. - First rear annular shield plate RAS1. 
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Figure 52. - Section B-B of figure 5. 
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Figure 53. - Control-rod details. 
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Figure 54. - Schematic elevation showing shaping of biological shield. 
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Figure 55. - Variation in shield weight with altitude. Reactor power, 100 megawatts; dose rate, 0.025 rem per hour at 
100- foot distance. (Additional weights for streaming, voids, and cooling not included.) 
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Figure 56. - Variation of shield weight with dose rate at 100-foot distance. 
Altitude, 30,000 feet; reactor power, 100 megawatts . (Additional weights for 
streaming, voids, and cooling not included.) 
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Figure 57. - Block diagram of basic nuclear-powered turboprop control loops. 
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Figure 58. - Turboprop control scheme utilizing turbine bypass power modulation. 
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Figure 60. - Curve used to determine sizes of concentric lines (eq. ( C2 7 ) ) - 
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Figure 61. - Estimated performance map of axial- flow 
helium pump assuming operation at constant actual 
blade speed. 
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Figure 62. - Axial cross section of fuel element. 
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Figure 63. - Illustration of method of determining ad- 
justed equivalent stress cr a ^. 
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(b) Results from diffusion- theory program. 
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(d) Final curves. 

Figure 64. - Method of determining effective multiplication factor k eff . 
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Figure 66. - Section F-F of figure 5. 
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Figure 67. - Sketch used in. determining direct neutron dose rate at detector. 
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